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1.0  OBJECTIVES 


This  final  report  documents  work  conducted  as  part  of  a  five-year  Multi-disciplinary  University 
Research  Initiative  (MURI)  in  the  area  of  Radiation  Effects  on  Emerging  Electronic  Materials  and 
Devices.  The  objectives  of  this  program  are  to  examine  and  understand  the  underlying  physical 
phenomena  that  control  the  radiation  response  of  semiconductor  technologies  incorporating 
emerging  materials  and  devices.  The  radiation  response  and  electrical  properties  of  technologies 
that  exhibit  exceptional  promise  for  application  in  DoD  systems  are  investigated  experimentally 
and  through  application  of  advanced  theory  and  simulation.  The  overall  purpose  is  to  develop 
knowledge  and  tools  that  will  guide  development  of  future  radiation-hardened  electronics.  This 
MURI  program  includes  researchers  from  Vanderbilt  University,  the  University  of  Florida,  Georgia 
Institute  of  Technology,  North  Carolina  State  University,  Rutgers  University,  and  Arizona  State 
University  and  strong  collaboration  with  leading  industrial  and  government  labs. 

The  combined  effects  of  advances  in  microelectronic  materials  and  device  structures  have  resulted 
in  more  changes  to  underlying  integrated-circuit  technologies  over  the  past  decade  than  in  the 
previous  forty  years.  As  integrated  circuits  become  smaller  and  more  complex,  it  is  increasingly 
difficult  to  predict  their  responses  to  radiation  and  to  guarantee  reliable  operation.  These  problems 
are  compounded  by  the  incorporation  of  new  materials  and  device  structures.  Before  the  1990s, 
only  seven  elements  were  in  common  use  in  integrated  circuits;  currently  more  than  half  the 
elements  in  the  periodic  table  are  in  use  or  under  evaluation.  Simultaneously,  device  dimensions 
have  reached  approximately  1  nm  vertically  and  several  10s  of  nm  horizontally.  The  interaction  of 
radiation  with  these  new  materials  in  such  small  volumes  can  cause  unexpected  and  challenging 
reliability  problems.  For  some  advanced  technologies,  the  error  rates  observed  in  satellite  systems 
have  been  orders  of  magnitude  greater  than  those  predicted  by  conventional  ground-based  testing 
and  simulation  methods.  Understanding  and  mitigating  these  radiation-induced  errors  is  the  key 
reliability  challenge  associated  with  deploying  future  generations  of  electronics  in  space  and 
defense  systems.  Energy  absorption,  carrier  generation,  carrier  transport,  charge  trapping  and 
defect  formation  and  dynamics  depend  on  the  specific  materials  used  in  the  ICs.  Sensitivity  to  the 
electrostatic  effects  of  radiation-induced  trapped  charge,  lifetime  degradation,  and  device-edge  and 
inter-device  leakage  depend  on  the  detailed  device  geometries  and  doping  profiles.  Moreover,  high¬ 
speed  circuits  exhibit  increased  vulnerabilities  to  single-event  effects,  including  multiple-bit  upsets 
that  result  from  aggressive  scaling,  and  large  enhancements  relative  to  ion-strike  angle  that  are 
much  greater  than  for  earlier  generations  of  technology.  While  it  was  possible  to  study  total  dose 
and  single  event  effects  separately  in  larger  devices,  the  boundary  is  now  increasingly  less  clear  as 
a  single  event  may  induce  charging  or  damage  in  the  entire  device. 

The  research  conducted  as  part  of  this  program  combines  experiments  on  emerging  materials  and 
devices  with  atomic-scale  materials  theory  and  a  high  performance  computing-based  radiation- 
effects  simulation  approach  to  develop  guidelines  for  design  and  application  of  advanced  electronic 
technologies  in  radiation  environments .  The  emphasis  is  on  technologies  that  show  the  most  promise 
for  defense  and  space  systems.  Experimental  samples  were  fabricated  in  the  labs  of  MURI  team 
members  and  acquired  through  collaborations  with  industrial  and  government  laboratories.  The 
technologies  examined  include  alternative  gate  dielectrics,  strained-layer  transistors,  silicon-on- 
insulator  (SOI)  devices,  FinFETs,  SiGe  devices,  and  ultra-small  CMOS  transistors.  Samples  were 
irradiated  using  particle  and  photon  sources  and  the  results  were  interpreted  using  an  approach  that 
combines  first-principles  quantum  mechanical  analysis  of  defects  with  engineering-level  modeling. 
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A  new  simulation  approach,  based  on  determining  the  device-level  response  to  an  ensemble  of 
realistic  particle-generated  track  structures,  was  used  to  analyze  single-event  effects  in  ultra- small 
devices.  This  approach,  which  employs  a  Monte  Carlo  tool  developed  at  Vanderbilt,  allows  the 
most  accurate  analysis  to  date  of  single-event  upsets  and  transients,  as  well  as  emerging  issues  like 
localized  displacement  damage  and  microdose  effects.  The  simulation  tool,  Monte  Carlo  Radiative 
Energy  Deposition  (MRED),  includes  the  GEANT4  libraries  developed  by  the  high  energy  physics 
community  as  well  as  other  models  for  the  interaction  of  radiation  with  matter.  The  approach  is 
complemented  by  atomic-scale  calculations  of  radiation-induced  atomic  displacements,  defect 
generation,  and  subsequent  charging. 

2.0  STATUS  OF  EFFORT 

The  program  ended  in  August  of  2010,  following  a  three-month  no-cost  extension.  Some  of  the 
major  accomplishments  of  this  program  include: 

•  The  groups  at  North  Carolina  State ,  Rutgers ,  and  Vanderbilt  collaborated  to  develop  the  most 
extensive  understanding  to  date  of  the  radiation  response  of  alternative  gate  dielectrics.  This 
is  a  very  significant  issue  for  application  of  future  CMOS  technology  generations  in  DoD 
systems.  These  new  gate  dielectrics,  particularly  those  based  on  hafnium,  were  introduced 
into  production  in  the  45 -nm  CMOS  technology  node.  The  higher  dielectric  constant  of 
hafnium-based  films  compared  to  silicon  dioxide  allows  a  greater  physical  thickness  for  the 
film  while  still  providing  sufficiently  high  capacitance  per  unit  area.  This  greater  thickness, 
however,  can  increase  the  sensitivity  to  total  ionizing  dose.  Among  the  most  important 
findings  is  that  HfSiON  films  have  superior  total-ionizing-dose  radiation  response  compared 
to  other  Hf-based  alternative  dielectrics.  In  addition,  the  combined  effects  of  radiation  and 
electrical  stress  were  found  to  be  worse  than  a  linear  combination  of  the  individual  effects  in 
some  films. 

•  North  Carolina  State  produced  Hf02  films  on  Ge  substrates  that  have  interfacial  properties 
equivalent  to  those  of  thermal  oxide  on  Si.  The  key  step  is  the  formation  of  the  interfacial 
layer.  This  material  system  is  important  because  of  the  potential  to  build  devices  with  better 
drive  current,  particularly  for  p-channel  transistors,  but  the  carrier- transport  properties  have 
been  limited  by  poor  interfacial  quality.  NCSU  provided  devices  to  Vanderbilt  for  radiation 
characterization.  Vanderbilt  conducted  detailed  radiation  studies  of  Hf02  on  Ge  devices  and 
found  that  the  total-dose  radiation  response  is  promising. 

•  Vanderbilt  performed  systematic  studies  of  defect  stability,  hydrogen  transport,  and 
interfacial  reactions  in  irradiated  Si02  on  Si  using  density  functional  theory.  DFT  is  the 
state-of-the-art  method  for  analyzing  defects  in  solids.  The  resulting  understanding  of  the 
properties  of  hydrogen  in  Si02/Si  is  the  most  thorough  examination  of  the  fundamental 
processes  responsible  for  radiation-induced  interface  traps  performed  to  date. 

•  Vanderbilt  developed  a  physical  model  of  bias-temperature  instability  (BTI)  in  Si  devices 
and  showed  that  the  processes  responsible  for  BTI  and  radiation-induced  degradation  have 
a  common  origin.  In  particular,  it  was  found  that  hydrogen  originating  in  the  Si  substrate 
plays  a  key  role  in  both  negative  bias  temperature  instability  and  enhanced  low-dose-rate 
sensitivity. 

•  The  University  of  Florida  and  Vanderbilt  collaborated  to  perform  the  first  controlled 
measurements  of  radiation  response  in  strained-Si  technologies.  Strained  Si  is  the  key 
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3.0  ACCOMPLISHMENTS/FINDINGS 


3.1  TID  Effects  on  1-T  DRAMs: 

We  have  characterized  the  total  ionizing  dose  response  of  two  SOI  FinFET  technologies  operated 
in  1T-DRAM  mode,  one  with  poly-crystalline  Si  gates  and  a  Si02  gate  dielectric,  and  the  other  with 
metal  gates  and  a  high-K  gate  dielectric.  These  devices  were  programmed  using  both  back-gate 
pulse  and  gate  induced  drain  leakage  (GIDL)  programming  methods.  1T-DRAM  cells  programmed 
with  back-gate  pulses  are  quite  sensitive  to  total  ionizing  dose  radiation,  with  the  memory  sensing 
margin  and  retention  time  decreasing  significantly  with  increasing  dose.  In  contrast,  the  sensing 
margin  and  retention  time  show  high  tolerance  to  total  ionizing  dose  (TID)  irradiation  when  the 
1T-DRAM  cells  are  programmed  via  gate  induced  drain  leakage  (GIDL);  see  Fig.  1.  The  sensing 
margin  increases  significantly  with  decreasing  gate  length  for  1T-DRAM  cells  programmed  via 
GIDL.  We  conclude  that  capacitor- less  lT-DRAMs  programmed  via  GIDL  are  strong  candidates  for 
embedded  memory  applications  for  ultimately  scaled  CMOS  devices  in  high-TID  applications.  [E. 
X.  Zhang,  D.  M.  Fleetwood,  F.  El  Mamouni,  R.  D.  Schrimpf,  M.  L.  Alles,  W.  Xiong,  K.  Akarvardar, 
and  S.  Cristoloveanu,  “Total  ionizing  dose  effects  on  FinFET-based  capacitor-less  lT-DRAMs,” 
accepted  for  publication  in  IEEE  Trans.  Nucl.  5c/.,  57,  no.  6,  Dec.  2010.] 

3.2  Low  Frequency  Noise,  Moisture  Effects,  and  Aging 

We  have  found  that  the  room  temperature  1/f  noise  gate-voltage  and  frequency  dependences  of 
pMOS  transistors  are  affected  significantly  by  moisture  exposure  and  total  dose  irradiation.  The 
voltage  noise  power  spectral  density  Syd  ~  (V  -  Vt)"p,  where  Vt  is  the  threshold  voltage,  Vg  is  the 
gate  voltage,  and  P  is  a  measure  of  the  gate- voltage  dependence.  For  pMOS  devices  (Fig.  2),  pre¬ 
irradiation,  P  ranges  from  0.4  to  0.9,  and  the  frequency  exponent,  a  =  -d  lnSyd/dln  f,  is  greater  than 
unity.  Post-irradiation,  the  gate- voltage  and  frequency  dependences  change  significantly,  with  p  » 
1,  and  a  much  closer  to  unity.  For  nMOS  devices  (Fig.  3),  pre-irradiation,  P  >  1.6  and  a  ~  1,  with 
little  change  after  irradiation.  We  attribute  these  observed  changes  in  pMOS  noise  to  changes  in  the 
trap  density  and  energy  distribution,  Dt(Ef),  of  these  devices.  Before  irradiation,  Dt(Ef)  increases 


Fig.  1.  Sensing  margin  as  a  function  of  total  ionizing  dose  for  nMOS  1T-DRAM  cells 
programmed  by  back-gate  pulse  and  GIDL  methods. 
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toward  the  valence  band  edge,  but  after  irradiation,  the  distribution  typically  is  more  uniform. 
Moreover,  for  some  moisture-exposed  devices,  Svd  p  ~(Vg-Vt)  3  after  irradiation,  indicating  a  Dt(Ef) 
that  increases  toward  midgap.  We  conclude  that  irradiation  and/or  moisture  exposure  can  greatly 
affect  the  defect  energy  distributions  for  these  devices,  and  that  the  observed  nMOS  and  pMOS 
noise  can  be  described  by  a  simple  trapping  model  with  an  energy-dependent  trap  distribution.  [S. 
A.  Francis,  A.  Dasgupta,  and  D.  M.  Fleetwood,  “Effects  of  total  dose  irradiation  on  the  gate-voltage 
dependence  of  the  1/f  noise  of  nMOS  and  pMOS  transistors,”  IEEE  Trans.  Electron  Dev.  57,  503- 
510(2010).] 


Frequency  (Hz) 


Fig.  2.  SVd  as  a  function  of  frequency  for  3  jam  x  16  jam  pMOS  transistors  from  control  and 
moisture-exposed  parts,  before  and  after  500  krad(Si02)  total  dose  irradiation. 


3.3  Defects  and  Hydrogen 

We  have  performed  an  extensive  investigation  of  the  energetics  of  the  interactions  between  molecular 
hydrogen  and  common  defects  in  Si02  that  are  typically  associated  with  O  deficiency.  We  find  that 
H2  does  not  easily  crack  at  neutral  vacancies,  but  it  will  crack  efficiently  at  O  vacancy  sites  that 
have  captured  a  hole  and  relaxed  into  the  puckered  configuration  of  an  E^  defect  (Fig.  4),  releasing 
a  proton  into  the  oxide.  This  leads  to  a  reduction  in  oxide-trap  charge  and  an  increase  in  interface- 
trap  charge  with  time.  Isolated  Si  dangling  bonds  also  can  play  a  role  in  cracking  H2  depending  on 
their  concentration.  These  results  provide  significant  insight  into  the  underlying  causes  of  latent 
interface  trap  buildup  in  MOS  devices  and  enhanced  low-dose-rate  sensitivity  in  linear  bipolar 
devices.  (Fig.  5.)  These  are  longstanding  problems  in  the  radiation  effects  community.  [B.  R.  Tuttle, 
D.  R.  Hughart,  R.  D.  Schrimpf,  D.  M.  Fleetwood,  and  S.  T.  Pantelides,  “Defect  interactions  of  H2  in 
Si02:  Implications  for  ELDRS  and  latent  interface  trap  buildup,”  accepted  for  publication  in  IEEE 
Trans.  Nucl.  Sci.,  57,  no.  6,  Dec.  2010.] 

We  have  also  found  that  enhanced  low-dose-rate  sensitivity  (ELDRS)  in  lateral  and  substrate  pnp 
bipolar  devices  can  occur  because  of  the  much  lower  probability  for  electron  capture  by  protons  in 
Si02,  as  compared  to  mobile  or  trapped  holes.  New  experimental  results  and  a  critical  evaluation 
of  previous  work  demonstrate  that,  at  high  dose  rates  and/or  in  oxides  with  low  concentrations  of 
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N 


Vg-Vt(V) 

Fig.  3.  SVd  at  ~10  Hz  as  a  function  of  Vg— Vt  for  3  jam  x  16  jam  nMOS  transistors  from  the 
control  and  moisture-exposed  parts,  before  and  after  500  krad(Si02)  total  dose  irradiation. 


Reaction  Coordinate  Number 


Fig.  4.  The  reaction  energy  for  H2  dissociating  at  a  positively  charged,  puckered  vacancy. 


hydrogen,  electrons  can  more  easily  neutralize  slowly  diffusing  or  metastably  trapped  holes  via 
annihilation  (recombination)  or  compensation  (offsetting  trapping)  before  the  holes  can  release 
H+.  In  contrast,  at  low  dose  rates  and/or  in  oxides  with  higher  concentrations  of  hydrogen,  which 
can  react  with  and  modify  the  structure  of  O- vacancy-related  defects,  it  is  more  likely  that  holes 
can  release  H+  during  transport.  This  is  because  the  cross  section  for  electron  capture  by  H+  is 
several  orders  of  magnitude  smaller  than  the  electron  capture  cross  section  for  a  slowly  moving  or 
metastably  trapped  hole  (Fig.  6).  This  enhanced  proton  release  at  low  dose  rates  or  in  oxides  with 
high  hydrogen  concentrations  can  lead  to  increased  interface  trap  formation,  which  is  the  most 
common  source  of  enhanced  gain  degradation  in  lateral  and  substrate  pnp  bipolar  transistors.  This 
provides  significant  insight  into  the  longstanding  problem  of  ELDRS  in  linear  bipolar  transistors. 
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Fig.  5.  Threshold  voltage  shifts  due  to  oxide-  and  interface-trap  charge  AVot  and 
AVit,  as  well  as  the  ratio  of  post-irradiation  to  pre-irradiation  mobility  p/ppre,  for 
pMOS  transistors  from  Oki  Semiconductor,  irradiated  to  75  krad(Si02)  at  a  dose 
rate  of  ~  240  rad(Si02)  / s  and  annealed  at  1 00  °C.  The  irradiation  and  annealing  bias 
was  +6  V.  AVot  and  AVit  were  estimated  via  the  midgap  charge  separation  technique. 
(After  J.  R.  Schwank,  et  al.,  IEEE  Trans.  Nucl.  Sci.,  vol.  39,  no.  6,  pp.  1953-1963,  Dec. 

1992.) 
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Fig.  6.  Relative  charge  neutralization  rates  for  trapped  holes  and  protons  in  ~  1  pm 
oxides  grown  on  zone-melt  recrystallized  Si02  during  5  eV  ultraviolet  illumination 
for  the  same  exposure  conditions  at  295  K  and  77  K.  This  energy  is  high  enough 
to  inject  electrons  over  the  3.1  eV  Si-to-Si02  barrier,  but  too  low  to  create  electron- 
hole  pairs  in  Si02.  (After  K.  Vanheusden,  et  al.  Appl.  Phjjs.  Lett.,  72,  no.  1,  pp.  28-30, 

1998.) 
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[D.  M.  Fleetwood,  R.  D.  Schrimpf,  S.  T.  Pantelides,  R.  L.  Pease,  and  G.  W.  Dunham,  “Electron 
capture,  hydrogen  release,  and  enhanced  gain  degradation  in  linear  bipolar  devices,”  IEEE  Trans. 
Nucl.  Sci.  55,  2986-2991  (2008).] 

3.4  Defects  in  High-K  Gate  Oxides 

We  have  characterized  radiation  induced  charge  trapping  in  ultrathin  Hf02-based  n-channel 
MOSFETs  as  a  function  of  dielectric  thickness  and  irradiation  bias  following  exposure  to  10  keV 
X-rays  and/or  constant  voltage  stress.  Positive  and  negative  oxide-trap  charges  are  observed, 
depending  on  irradiation  and  bias  stress  conditions.  No  significant  interface-trap  buildup  is  found 
in  these  devices  under  these  irradiation  and  stress  conditions.  Enhanced  oxide-charge  trapping 
occurs  in  some  cases  for  simultaneous  application  of  constant  voltage  stress  and  irradiation, 
relative  to  either  type  of  stress  applied  separately  (Fig.  7).  Room  temperature  annealing  at  positive 
bias  after  irradiation  of  transistors  with  thicker  gate  dielectric  films  leads  to  positive  oxide-trapped 
charge  annihilation  and/or  neutralization  in  these  devices,  and  net  electron  trapping.  The  oxide 
thickness  dependence  of  the  radiation  response  confirms  the  extreme  radiation  tolerance  of  thin 
Hf02  dielectric  layers  of  relevance  to  device  applications,  and  suggests  that  hole  traps  in  the  thicker 
layers  are  located  in  the  bulk  of  the  dielectric.  A  revised  methodology  is  developed  to  estimate  the 
net  effective  charge  trapping  efficiency,  f  t,  for  high-K  dielectric  films.  As  a  result,  estimates  of  f  t 
for  Hf  silicate  in  previous  work  by  Felix  et  al.  (2002)  are  reduced  by  up  to  18  %.  [S.  K.  Dixit,  X.  J. 
Zhou,  R.  D.  Schrimpf,  D.  M.  Fleetwood,  S.  T.  Pantelides,  R.  Choi,  G.  Bersuker,  and  L.  C.  Feldman, 
“Radiation  induced  charge  trapping  in  ultrathin  Hf02-based  MOSFETs,”  IEEE  Trans.  Nucl.  Sci. 
54,1883-1890  (2007).] 

Charge  trapping  characteristics  were  investigated  for  MOS  capacitors  with  6.8  nm  Hf02  layers  and 
1.0  nm  interfacial  silicon  oxynitrides;  the  effective  oxide  thickness  of  the  high-K  gate  dielectric 
layers  is  2.1  nm.  These  devices  were  irradiated  with  10-keV  X-rays  or  subjected  to  constant 
voltage  stress,  and  then  annealed  for  ten  minute  intervals  of  alternating  positive  and  negative  gate 
bias  at  temperatures  between  25  and  150  °C.  The  resulting  oxide-trap  (Fig.  8)  and  interface-trap 
charge  (Fig.  9)  densities  exhibit  reversible  buildup  and  annealing  that  depend  strongly  on  bias  and 
temperature.  Additional  defect-density  growth  with  time  was  observed  as  a  result  of  charge  injection 
into  the  gate  stack  during  the  annealing  process  (Fig.  10).  This  defect-density  growth  increases  with 
increasing  annealing  time  and  temperature.  After  irradiation,  the  most  of  the  reversibility  in  the 
charge  trapping  is  due  to  metastable  electron  traps  in  the  near-interfacial  dielectric  layers.  After 
constant  voltage  stress,  the  motion,  reactions,  and  trapping  of  protons  at  or  near  the  Si/oxynitride 
interface  are  more  important  to  the  observed  device  response  than  are  metastable  electron  traps. 
This  is  a  result  of  reduced  electron-hole  pair  creation  during  low-energy  constant- voltage  stress,  as 
compared  to  high-energy  X-ray  irradiation.  These  results  illustrate  the  importance  of  both  electron 
traps  and  protons  to  the  ionizing  radiation  response  and  long-term  reliability  of  MOS  devices  with 
high-K  gate  dielectrics.  [X.  J.  Zhou,  D.  M.  Fleetwood,  L.  Tsetseris,  R.  D.  Schrimpf,  and  S.  T. 
Pantelides,  “Effects  of  switched-bias  annealing  on  charge  trapping  in  Hf02  gate  dielectrics,”  IEEE 
Trans.  Nucl.  Sci.  53, 3636-3643  (2006).] 
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Fig.  7.  Threshold  voltage  shifts  (AVT)  due  to  constant  voltage  stress,  CVS  (-2  V),  and  CVS  (-2  V) 
+  irradiation.  The  bottom  x-axis  shows  the  dose  in  Mrad(Si02).  The  top  x-axis  shows  the  time 

during  CVS. 


Fig.  8.  AVt  (a)  for  Al/Hf02+SiOxNy/Si  pMOS  capacitors  irradiated  to  1.0  Mrad(Si02)  with 
10-keV  X-rays,  followed  by  a  series  of  switched  bias  anneals  at  50  to  150  °C.  The  gate  bias  for 
irradiation  is  0.3  V.  The  switched  bias  anneals  are  ±  0.3  V  (PBTS,  NBTS),  and  the  stress  time  was 

600  s  each. 
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Switched  Bias  Experiments 


Fig.  9.  AV.t  for  Al/Hf02+SiOxNy/Si  pMOS  capacitors  irradiated  to  1.0  Mrad(Si02)  with  10- 
keV  X-rays,  followed  by  a  series  of  switched  bias  anneals  at  50  to  150  °C.  The  gate  bias  for 
irradiation  is  0.3  V.  The  switched  bias  anneals  are  ±  0.3  V  (PBTS,  NBTS),  and  the  stress  time 

was  600  s  each. 


-rad 


Switched-bias  experiments 


Fig.  10.  Cumulative  injected  charge  densities  for  switched-bias  annealing  at  elevated 
temperatures  after  irradiation.  The  capacitor  area  A  =  5.6  '  10-4  cm2.  The  experimental 
conditions  are  the  same  as  Fig.  1 .  Cumulative  charge  densities  during  PBTS  are  shown 
in  the  upper  part  of  the  figure,  and  cumulative  charge  densities  during  NBTS  are  shown 
in  the  lower  half.  These  densities  do  not  include  the  portions  of  the  anneal  when  the 
device  was  in  the  opposite  bias  state,  so  the  total  charge  is  obtained  by  adding  the  charge 
densities  in  the  upper  and  lower  halves  of  the  figure. 
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3.5  A  Comprehensive  Understanding  of  the  Efficacy  of  N-Ring  SEE 
Hardening  Methodologies  in  SiGe  HBTs 

Georgia  Tech’s  contributions  to  this  five  year  MURI  on  radiation  effects  centered  on  the  radiation 
response  and  damage  mechanisms  in  device  and  circuits  built  from  silicon-based  heterostructures, 
including  primarily  silicon-germanium  (SiGe)  heterojunction  bipolar  transistors  (HBTs),  SiGe 
modulation  doped  FETs  (MODFETs),  and  strained- silicon  CMOS. 


We  investigate  the  efficacy  of  mitigating  radiation-based  single  event  effects  (SEE)  within  circuits 
incorporating  SiGe  heterojunction  bipolar  transistors  (HBTs)  built  with  an  N-Ring,  a  transistor-level 
layout-based  radiation  hardened  by  design  (RHBD)  technique.  Previous  work  of  single-device  ion- 
beam  induced  charge  collection  (IB ICC)  studies  has  demonstrated  significant  reductions  in  peak 
collector-collected  charge  and  sensitive  area  for  charge  collection;  however,  few  circuit  studies 
using  this  technique  have  been  performed.  Transient  studies  performed  with  Sandia  National 
Laboratory’s  (SNL)  36  MeV  160  microbeam  on  voltage  references  built  with  SiGe  N-ring  SiGe 
HBTs  have  shown  mixed  results,  with  reductions  in  the  number  of  large  voltage  disruptions  in 
addition  to  new  sensitive  areas  of  low-level  output  voltage  disturbances.  Similar  discrepancies 
between  device-level  IBICC  results  and  circuit  measurements  are  found  for  the  case  of  digital 
shift  registers  implemented  with  N-ring  SiGe  HBTs  radiated  in  a  broadbeam  environment  at 
Texas  A&M’s  Cyclotron  Institute.  The  over-all  error  cross-section  of  the  N-ring-based  register 
is  found  to  be  larger  than  the  standard  SiGe  register,  which  is  clearly  counter-intuitive.  We  have 
worked  to  resolve  the  discrepancy  between  the  measured  circuit  results  and  the  device-level 
IBICC  measurements,  by  re-measuring  single-device  N-ring  SiGe  HBTs  using  a  time-resolved 
ion  beam  induced  charge  (TRIBIC)  set-up  that  allows  direct  capture  of  nodal  transients.  Coupling 
these  measurements  with  full  3-D  TCAD  simulations  provides  complete  insight  into  the  origin  of 
transient  currents  in  an  N-ring  SiGe  HBT.  The  detailed  structure  of  these  transients  and  their  bias 
dependencies  are  discussed,  together  with  the  ramifications  for  the  design  of  space-borne  analog 
and  digital  circuits  using  SiGe  HBTs. 
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Fig.  11.  RHBD  approach  using  n-rings  for  SiGe  HBTs. 
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Fig  12.  Simulated  transient  waveforms  for  an  ion  strike  through  the  N-ring  of  an  N-ring  SiGe  HBT  biased  with  3 
V  on  the  N-ring  and  collector  and  all  other  terminals  ground. 


3.6.  Single  Event  Transient  Hardness  of  a  New  Complementary 
(npn  +  pnp)  SiGe  HBT  Technology  on  Thick-film  SOI 

We  report  heavy-ion  microbeam  and  total  dose  data  for  a  new  complementary  (npn  +  pnp)  SiGe 
on  thick-film  SOI  technology.  Measured  transient  waveforms  from  heavy-ion  strikes  indicate 
a  significantly  shortened  transient  upset  current,  while  maintaining  the  total  dose  robustness 
associated  with  SiGe  devices.  Heavy-ion  broad-beam  data  confirms  a  reduced  single  event  upset 
(SEU)  cross-section  in  a  high-speed  shift  register  circuit. 


Fig.  13.  SEM  cross-section  of  a  SiGe  HBT  on  thick  film  SOI. 
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Fig.  14.  Measured  improvement  in  SEU  cross-section  and  LET  for  SiGe  on  SOI. 


3.7  Single  Event  Transient  Response  of  SiGe  Voltage  References 
and  Its  Impact  on  the  Performance  of  Analog  and  Mixed-Signal 
Circuits 

We  investigated  the  single-event  transient  (SET)  response  of  bandgap  voltage  references  (BGRs) 
implemented  in  SiGe  BiCMOS  technology  through  heavy  ion  microbeam  experiments.  The  SiGe 
BGR  circuit  was  used  to  provide  the  input  reference  voltage  to  a  voltage  regulator.  SiGe  HBTs 
in  the  BGR  circuit  were  struck  with  36  MeV  oxygen  ions  and  the  subsequent  transient  responses 
were  captured  at  the  output  of  the  regulator.  Sensitive  devices  responsible  for  generating  transients 
with  large  peak  magnitudes  were  identified.  To  determine  the  effectiveness  of  a  transistor-layout- 
based  radiation  hardened  by  design  (RHBD)  technique  with  respect  to  immunity  to  SETs  at  the 
circuit  level,  the  BGR  circuit  implemented  with  HBTs  surrounded  by  an  alternate  reverse-biased 
pn  junction  (n-ring  RHBD),  were  also  bombarded  with  oxygen  ions  and  subsequent  SETs  captured. 
Experimental  results  indicated  that  the  number  of  events  with  large  peak  magnitude  have  been 
reduced  in  the  RHBD  version;  however,  with  the  inclusion  of  the  n-ring  RHDB,  new  locations  for 
the  occurrence  of  transients  (albeit  with  smaller  peak  magnitude)  were  created.  Transients  at  the 
transistor  level  were  also  independently  captured.  It  was  demonstrated  that  while  the  transients  are 
short  at  the  transistor  level  (ns  duration),  relatively  long  transients  are  obtained  at  the  circuit  level 
(100’s  of  ns).  In  addition,  the  impact  of  the  SET  response  of  the  BGR  on  the  performance  of  an 
ultra-high  speed  3-bit  SiGe  analog-to-digital  converter  (ADC)  was  investigated  through  simulation. 
It  was  shown  that  ion-induced  transients  in  the  reference  voltage  could  eventually  lead  to  data 
corruption  at  the  output  of  the  ADC. 


3.8  Re-examining  TID  Hardness  Assurance  Test  Protocols  for  SiGe 
HBTs 

We  investigated  the  applicability  of  current  TID  test  protocols  in  the  context  of  advanced  transistor 
technologies  such  as  SiGe  HBTs.  In  SiGe  HBTs,  an  unexpected  shift  in  collector  current  was 
observed  during  total  dose  irradiation.  Using  both  device  and  circuit  measurements,  we  investigated 
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this  phenomenon  and  assessed  its  potential  importance  in  hardness  assurance  of  SiGe  components. 
TCAD  simulations  were  performed  to  explain  the  observed  current  shifts. 


Fig.  15.  Collector  and  base  transient  currents  for  strikes  over 
and  outside  the  emitter  area  for  standard  HBTs. 


Fig.  16.  Single-event  transient  response  of  the  regulator  for  three 
different  strike  locations. 


Fig.  17.  Collector  and  base  currents  of  a  SiGe  HBT  under 
irradiation  of  1  Mrad.  Radiation  stopped  at  1,275  s,  but  biasing 
continued  until  a  total  of  1,600  s  was  reached. 
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Fig.  18.  Normalized  collector  current  during  gamma  TID 
testing.  Each  device  was  forward-biased  at  the  same  VBE  and 
VCE,  and  the  dose  rates  were  17.6,  11.8,  2.08,  1.17,  and  0.22 
rad(Si)  /  s. 

3.9  A  Novel  Device  Architecture  for  SEU  Mitigation:  The  Inverse- 
Mode  Cascode  SiGe  HBT 

We  investigated,  for  the  first  time,  the  potential  for  SEE  mitigation  of  a  newly-developed  device 
architecture  in  a  3rd  generation  high-speed  SiGe  platform.  This  new  device  architecture  is  termed 
the  “inverse-mode  cascode  SiGe  HBT”  and  is  composed  of  two  standard  devices  sharing  a  buried 
subcollector  and  operated  in  a  cascode  configuration.  Verification  of  the  TID  immunity  was 
demonstrated  using  10  keV  x-rays,  while  an  investigation  of  the  SEE  susceptibility  was  performed 
using  a  36  MeV  160  ion.  IBICC  results  show  strong  sensitivities  to  device  bias  with  only  marginal 
improvement  when  compared  to  a  standard  device;  however,  by  providing  a  conductive  path  from 
the  buried  subcollector  (C-Tap)  to  a  voltage  potential,  almost  all  “collected”  charge  is  induced 
on  the  C-Tap  terminal  instead  of  the  collector  terminal.  These  results  were  confirmed  using  full 
3D  TCAD  simulations  which  also  provides  insight  into  the  physics  of  this  new  RHBD  device 
architecture.  The  implications  of  biasing  the  C-Tap  terminal  in  a  circuit  context  are  also  addressed. 
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Fig.  19.  (a)  Cross-section  of  a  representative  inverse-mode  cascode  SiGe  HBT  device  showing  active  layers  as  well 
as  the  modification  needed  to  form  the  second  variant  ofthe  device  with  C-Tap  in  addition  to  the  (b)  schematic 

of  the  inverse-mode  cascode  SiGe. 
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Fig.  20.  Measured  2-D  cut  of  the  electrical  collector-collected  charge  data  for  the  inverse 
cascode  with  C-Tap  showing  nearly  complete  mitigation  of  charge  collection  on  the  electrical 

collector  terminal. 

3.10  Novel  Total  Dose  and  Heavy-Ion  Charge  Collection  Phenomena 
in  a  New  SiGe  HBT  on  Thin-Film  SOI  Technology 
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We  investigated  the  impact  of  63  MeV  proton  radiation  (up  to  a  total  dose  of  2  Mrad(Si02))  on 
the  AC  and  DC  performance  of  a  new  high-performance  SiGe  HBT-on-SOI  technology  from 
STMicroelectronics.  For  the  first  time,  we  compare  the  radiation  response  of  this  SiGe  HBT-on- 
SOI  with  that  of  a  bulk  SiGe  HBT  fabricated  with  an  identical  emitter-base  structure.  That  is,  the 
only  differences  between  the  devices  are  the  substrate  (bulk  vs.  SOI)  and  the  collector  doping. 
More  importantly,  the  devices  under  study  feature  an  innovative  CBEBC  layout  (with  off-plane 
base  contacts),  which  significantly  improves  the  AC  performance.  This  work  analyzes  for  the  first 
time  the  impact  of  the  novel  CBEBC  layout  on  both  Total  Ionizing  Dose  (TID)  and  SEU  response. 
Experimental  data  and  calibrated  3D  TCAD  simulations  demonstrate  that,  in  the  inverse  mode, 


Fig.  21.  Cross-sectional  TEM  micrograph  of  the  SiGe  HBT- 
on-SOI  with  CrEbC  layout. 
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the  current  flow  along  the  large  oxide  surface  between  the  collector  and  the  base  can  be  altered  by 
substrate  bias,  reducing  the  radiation-induced  leakage.  In  addition,  calibrated  3D  TCAD  simulations 
of  heavy  ion  strikes  in  the  center  of  the  emitter  indicate  that  the  CBEBC  layout  is  characterized  by 
novel  charge  collection  phenomena.  Finally,  for  the  first  time,  the  thermal  resistances  (R^)  of  the 
bulk  and  SOI  SiGe  HBTs  have  been  compared  before  and  after  irradiation  over  the  temperature 
range  of  300  K  to  390  K,  demonstrating  that  radiation  exposure  increases  RTH  in  SOI  devices. 


Fig.  22.  ID  plots  of  power  density  in  a  SiGe  HBT  on  SOI  for  V§  =  0 
V  and  V§  =  20  V,  along  the  line  ^  indicated  in  the  inset. 


3.11  On  the  Radiation  Tolerance  of  SiGe  HBT  and  CMOS-based 
Phase  Shifters  for  Space-based,  Phased-Array  Antenna 
Systems 

We  report  the  first  irradiation  results  on  high-frequency  SiGe  HBT  and  CMOS  phase  shifters  for 
space-  based,  phased-array  antennas  used  in  radar  or  wireless  communications  systems.  Both  phase 
shifter  circuits  remain  functional  with  acceptable  dc  and  RF  performance  up  to  multi-Mrad  proton 
exposure,  and  are  thus  suitable  for  many  orbital  applications.  In  addition,  simulation  results  probing 
the  limits  of  phase  shifter  performance  in  a  radiation  environment  are  presented.  The  CMOS  phase 
shifter  simulations  show  a  large  immunity  to  changes  caused  by  increased  leakage  currents,  while 
the  SiGe  HBT  phase  shifter  exhibits  only  marginal  performance  degradation  caused  by  shifts  in  the 
bias  currents  due  to  leakage  in  the  bias  networks. 

3.12  Dose  Rate  Effects  In  SiGe  HBTs 

As  an  ongoing  study  of  radiation  dose  rate  effects  in  advanced  SiGe  HBTs,  we  have  refined  our 
previous  results  with  better  statistics  with  a  large  amount  of  data  from  gamma,  proton,  and  x-ray 
radiation  sources,  with  dose  rates  ranging  from  0.1  rad(Si)/sec  to  1  krad(Si)/sec.  Shown  below 
are  excess  base  current  density  characteristics  from  3  different  radiation  sources  extracted  from 
forward  and  inverse  mode  operation,  respectively. 

From  the  collapsed  data,  we  found  that  the  DC  characteristics  showed  noticeable  (and  unexpected) 
dependencies  on  both  dose  rates  and  radiation  source  type.  Moreover,  the  emitter-base  (E-B)  spacer 
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Fig.  23.  Micrograph  of  the  SiGe  HBT  (left)  and  CMOS  (right)  phase 
shifters. 


Fig.  24.  Change  in  insertion  loss  between  pre-  and  post-irradiation  of  500  krad,  1 
Mrad,  and  3  Mrad  for  both  MOS  and  HBT  phase  shifters. 

and  shallow  trench  isolation  regions  seem  to  react  differently  and  do  not  provide  a  good  correlation 
with  macroscopic  dose  and/or  dose  rate  from  different  sources.  Therefore,  as  suggested  in  our 
previous  report,  microdosimetry  effects  in  SiGe  HBTs  are  being  examined,  with  the  possibility 
of  non-equilibrium  conditions  in  the  secondary  electrons  resulting  from  interactions  with  created 
photon  field.  We  have  found  that  the  x-ray  dose-rate  trends  are  qualitatively  consistent  with  the 
observed  inverse-mode  A JB;  that  is,  increasing  the  x-ray  dose  rate  yields  a  considerably  larger 
inverse-mode  noise  power  spectral  density.  We  attribute  the  dose  enhancement  effects  observed 
from  x-ray  irradiated  devices  to  photon  interaction  with  Cu/W  metallization  and  these  results 
show  a  good  qualitative  agreement  with  GEANT4  -MRED  simulations  on  SRAMs  (conducted 
by  Vanderbilt).  It  should  be  noted  that  the  E-B  spacer  is  a  complicated  oxide  nitride  composite 
stack,  whereas  the  STI  is  a  CVD  oxide  deposited  in  an  RIE  etched  trench  -  ie,  fundamentally 
different  in  composition  and  hence  potential  radiation  response.  These  types  of  films  can  trap  both 
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electrons  and  holes,  depending  on  the  processing  conditions,  the  radiation  dose  and  source,  and  the 
radiation  bias  condition.  We  conclude  that  fundamental  differences  exist  in  a  manner  by  which  the 
local  recombination  rates  in  the  E-B  and  STI  interface  regions  demonstrate  a  differing  response  to 
the  radiation  type  and  dose  rate.  For  instance,  charge  trapping  in  the  STI  may  effectively  “push” 
the  spatial  recombination  rate  peak  away  from  the  STI  region  and  into  the  bulk  Si  of  the  CB 
junction,  where  x-rays  do  not  change  recombination  rates,  but  protons  (and  gamma  rays)  may. 
In  this  scenario,  the  x-ray  response  of  the  device  would  “appear”  to  be  decreased  with  respect  to 
the  proton/gamma  response.  We  are  continuing  to  probe  these  mechanisms  via  calibrated  TCAD 
simulations. 
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Fig.  25.  Forward  mode  AJB  for  proton,  gamma,  and  X-ray  as  a  function  of  equivalent  gamma 

dose. 


Fig.  26.  Inverse  mode  AJB  for  proton,  gamma,  and  X-ray  as  a  function  of  equivalent  gamma 

dose. 
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3.13  Investigations  of  Cryogenic  Radiation  Effects  in  SiGe  HBTs 


Cryogenically-operated  devices  simultaneously  exposed  to  irradiation  is  a  newly  emerging  topic 
special  interest  to  NASA’s  Lunar  and  Martian  robotics  and  human  exploration  missions  as  well 
as  DoD  for  orbital  space  systems.  To  investigate  total  dose  effects  in  SiGe  HBTs  at  cryogenic 
temperatures,  we  have  irradiated  SiGe  HBTs  from  both  IBM  1st  generation  (5AM)  and  3rd  generation 
(8HP)  technologies  at  77K  with  63  MeV  protons,  and  measured  the  DC  characteristics  at  both  77K 
and  room  temperature.  Shown  below  are  the  forward  Gummel  characteristics  from  both  8HP  and 
5  AM  HBTs,  respectively,  irradiated  at  300K  and  77K. 


Fig.  27.  Forward  Gummel  characteristics  of  5AM  and  8HP  SiGe  HBTs, 
respectively,  from  proton  radiation  performed  at  room  temperature  and  77K. 
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3.14  Advanced  Simulation  Tools  for  Radiation  Effects 


The  Florida  Object  Oriented  Device  Simulator  (FLOODS)  is  a  general  purpose  device  simulation 
tool  built  that  is  capable  of  1,  2,  and  3  dimensional  simulation.  In  the  past  it  has  been  used  to 
simulate  SiGe  RF  bipolar  devices,  noise  in  advanced  CMOS  structures,  GaN  chemical  sensors, 
and  infrared  detectors.  It  can  handle  a  wide  variety  of  device  types  because  the  governing  physics 
is  not  hard-coded.  The  user  can  define  scripts  in  a  simple  language,  alagator,  that  describes  the 
governing  equations  for  the  system.  This  makes  it  relatively  straightforward  to  handle  different 
physical  situations  out  of  the  mainstream  of  CMOS  logic  devices. 

We  have  been  enhancing  FLOODS  to  handle  SEU  simulations.  This  has  included  activity  in  both 
the  numerical  and  physical  approximations.  The  numerical  approximation  work  has  centered 
around  trying  to  create  simpler  three-dimensional  adaptation  during  simulation.  SEU  problems 
deposit  a  large  charge  that  requires  fine  grid  spacing  and  then  diffuses  away,  which  reduces  the  grid 
requirement.  In  addition,  automatic  simulation  of  thousands  of  events  requires  the  grid  to  adapt  to 
the  charge  distribution  automatically.  On  the  physical  modeling  side,  we  have  focused  on  mobility 
updates  and  work  on  hydrogen  modeling. 

3.15  Numerical  Enhancements 

The  major  immediate  issue  was  enabling  transient  simulation.  FLOODS  had  been  developed 
primarily  for  DC  and  AC  small-signal  simulation,  both  of  which  are  inappropriate  for  SEU 
simulation  needs.  FLOODS  was  enhanced  to  include  transient  simulation  using  the  TRBDF2 
method.  TRBDF2  has  been  used  in  device  and  process  simulators  for  over  twenty  years.  The 
method  is  A-stable  and  includes  techniques  to  automatically  adapt  the  time  step  size  based  on  local 
truncation  error.  Algorithms  already  existed  inside  alagator  to  do  TRBDF2,  but  these  needed  to  be 
connected  to  the  device  command  and  contact  specifications. 

We  also  explored  parallel  processing,  although  the  results  did  not  appear  to  be  promising.  We 
completed  a  port  to  an  IBM  cluster  under  the  AIX  operating  system.  One  of  the  linear  algebra 
packages  in  FLOODS  does  have  parallel  compute  capability  and  we  enabled  that  and  create  a 
threaded  version  of  FLOODS.  The  alagator  scripting  language  is  coded  to  take  full  advantage  of 
modern  pipelined  processors.  Operators  are  applied  to  vectors  of  results  from  multiple  elements  at 
once.  Precomputation  and  storage  of  results  is  already  enabled  to  prevent  redundant  calculations. 
For  small-scale  parallelism,  we  can  subdivide  the  problem  physically.  One  processor  can  assemble 
Poisson’s  equation,  another  can  assemble  electron  continuity,  and  so  on.  Depending  on  the  modeling 
complexity,  this  efficient  for  up  eight  processors.  Breaking  the  assembly  this  way  is  attractive 
because  there  are  very  few  memory  overlap  issues.  The  code  is  already  structured  in  some  ways 
to  take  advantage  of  this  type  of  parallelism.  However,  our  work  here  was  only  able  to  provide  a 
factor  of  2  improvement  on  a  four  processor  system. 

The  set  of  coupled,  time-dependent  partial  differential  equations  that  govern  semiconductor  device 
behavior  are  usually  expressed  as  Poisson’s  equation  and  drift-diffusion  continuity  equations  for 
electrons  and  holes.  These  expressions  are  normally  evaluated  in  a  numerical  device  simulator  using 
Scharfetter-Gummel  (SG).  The  SG  approach  is  very  stable,  but  it  has  the  limitation  the  current  is 
not  defined  continuously.  For  most  simulations,  this  is  not  an  issue  as  a  grid  can  be  constructed 
that  has  edges  aligned  with  the  current  flow.  This  is  not  possible  in  the  case  SEU,  because  particle 
current  can  be  generated  anywhere  in  the  device  from  a  radiation  strike.  A  change  of  variables  can 
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allow  the  particle  currents  to  be  expressed  in  terms  of  the  quasi-fermi  levels  and  then  this  can  be 
discretized  directly.  This  approach  allows  the  current  to  be  defined  continuously  over  space. 

Our  work  on  this  problem  shows  that  the  quasi-fermi  level  approach  works  well  and  often  offers 
a  computational  advantage  over  the  SG  approach.  For  some  cases,  a  30%  improvement  in 
performance  could  be  noted  for  the  quasi-fermi  level  approach  and  only  rarely  did  this  ever  exceed 
the  SG  computation  time.  Solution  accuracy  was  equivalent  for  equivalent  grid  spacings  in  the 
simulation. 

We  have  also  developed  a  nested  grid  scheme  that  is  specific  to  SEU  simulation.  After  a  radiation 
strike,  we  know  that  the  peak  carrier  concentration  exists  at  that  time.  From  there  on  time,  the 
carrier  concentrations  only  reduce.  As  the  carrier  concentration  reduces,  we  can  also  reduce  to  a 
coarser  grid.  Figure  28  shows  our  scheme,  which  requires  successive  grid  refinement  along  the 
particle  strike  at  time  zero.  Multiple  passes  at  grid  refinement  are  used  and  each  pass  is  stored  for 
future  reference.  As  the  carrier  concentration  spreads  and  the  peak  reduces,  each  prior  refinement 
is  used  in  turn.  So  the  grid  eventually  coarsens  back  to  the  original  structure.  Figure  29  shows  what 
the  grid  looks  like  after  three  levels  of  refinement. 

Figure  30  shows  the  computation  result  of  this  process.  The  vertical  scale  represents  error  as 
measured  by  the  change  in  the  total  collected  charge.  The  horizontal  scale  represents  the  simulation 
CPU  time  normalized  to  the  fastest  simulations.  As  expected,  the  blue  line  for  the  uniform  grid 
shows  that  error  increases  and  CPU  time  decreases  as  the  grid  is  reduced.  Uniform  grids  are 
rarely  used,  so  the  green  line  represents  our  best  attempt  at  generating  a  custom  grid.  The  red  line 
represents  results  from  the  adaptive  grid  scheme.  The  same  level  of  collected  charge  accuracy  can 
be  obtained  with  an  order  of  magnitude  less  CPU  time. 


Fig.  28.  Proposed  grid  refinement  algorithm. 
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Fig.  29.  Adaptive  grid  at  peak  refinement  level  (X=3)  for  the  N+/P  diode  simulations. 
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Fig.  30.  nMOSFET  SEU  results.  The  number  of  m  grid  points  is  given  next  to  each  data  point.  The  results 
were  normalized  and  a  value  of  1 .0  on  both  scales  represents  the  optimal  result. 
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3.16  Physical  Model  Enhancements: 


Most  of  our  work  on  the  physical  side  focused  on  improvements  in  mobility  modeling  for  SEU 
simulation.  There  have  been  three  main  areas  -  improved  handling  of  piezoresistance  changes, 
carrier-carrier  scattering,  and  a  unified  approach  for  mobility. 

Advanced  devices  have  built-in  mechanical  strain  to  enhance  channel  mobility.  This  is  easy  to 
model  for  a  MOSFET,  because  the  primary  current  path  is  fixed  along  the  silicon  /  oxide  interface. 
The  strain  can  be  computed  here  and  a  mobility  enhancement  computed  for  the  channel  based  on 
the  strain.  This  approach  will  not  work  for  SEU,  because  current  flow  can  be  in  arbitrary  directions. 
In  addition,  the  strain  changes  over  the  device  structure  significantly  so  current  flow  generated  well 
below  the  surface  by  the  strike  cannot  use  the  channel  strain  to  estimate  mobility. 

We  have  implemented  the  entire  piezoresistance  matrix  and  related  the  changes  in  mobility  to  both 
the  local  strain  and  current  flow  direction.  This  couples  more  easily  with  the  quasi-fermi  level 
approach  with  finite  elements,  because  the  basis  function  for  computing  current  and  mechanical 
strain  are  identical  and  easy  to  map  on  to  each  other.  In  the  general  case,  we  use  the  gradient 
of  the  quasi-fermi  level  to  get  the  local  current  direction.  This  is  then  mapped  to  appropriate 
crystallographic  directions  and  transformed  geometrically  to  a  piezoresistance  change  which  is 
implemented  as  a  change  in  the  carrier  effective  mass. 

Figure  3 1  shows  a  comparison  between  experiment  and  simulated  collected  charge  as  a  function 
of  applied  strain  to  the  device  In  this  experiment,  p-n  diodes  were  used  as  the  collector  and  charge 
was  generated  from  a  laser  strike.  The  wafers  were  bent  externally  using  a  calibrated  four-point 
bending  jig.  Therefore  the  strain  was  able  to  be  varied  independent  of  changes  in  the  device 
structure.  Collected  charge  varies  with  the  applied  strain  nearly  linearly.  Simulation  of  the  strain 
influence  on  mobility  correlated  well  with  the  observed  experimental  results. 
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Fig.  31.  Collected  Charge  from  experiment  and  simulation  as  function  of  externally 
applied  strain. 
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One  advantage  of  simulation  is  that  it  can  predict  effects  beyond  those  that  can  be  reached 
experimentally.  Figure  32  shows  the  simulated  current  spikes  at  a  wide  range  of  strains,  beyond 
what  can  be  reached  with  the  four  point  bending  jig  without  breaking  the  samples.  The  change 
between  1  GPa  of  compressive  stress  to  1  GPa  of  tensile  stress  results  in  about  a  70%  change  in 
the  peak  current  value.  lGPa  is  typical  of  stress  being  used  in  advanced  strain  engineered  device 
structures. 


Time  (ns) 


Fig.  32.  Simulated  laser-induced  current  transients  as  a  function  of  <11 0>  uniaxial 
mechanical  stress. 


On  the  physical  side,  we  have  implemented  Klaassen’s  mobility  model  that  includes  terms  for 
minority  carrier  scattering.  These  are  important  and  can  be  controlling  for  SEU  event  simulations. 
During  our  work,  we  discovered  that  the  Klaassen’s  mobility  model  does  not  account  for  carrier- 
carrier  scattering  at  high  carrier  concentrations  accurately.  This  is,  of  course,  the  exact  condition 
that  a  SEU  simulation  faces. 

We  added  a  carrier-carrier  scattering  term  to  the  Klaassen  model  for  simulation  of  SEU.  Figure  33 
shows  the  effect  of  carrier-carrier  scattering  and  compares  models  to  experimental.  As  previously 
discussed,  the  Klaassen  /  Philips  model  highly  overestimates  mobility  at  electron-hole  levels  over 
1017  cnr3.  In  contrast,  the  Dorkel-Leturcq  model  uses  a  similar  approach  to  carrier-carrier  scattering 
as  the  proposed  model.  The  Dorkel-Leturcq  model  fits  well  for  lower  carrier  concentrations  but  at 
higher  level  begins  to  under  predict  mobility.  Another  issue  is  that  at  high-injection  levels  of  more 
than  5xl019cm~3,  the  Dorkel-Leturcq  model  predicts  a  negative  mobility.  Our  new  model  fits  this 
region  quite  well. 
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Fig.  33.  Sum  of  electron  and  hole  mobility  as  a  function  of  carrier  concentration  at  300  K  versus  experimental  data. 

Symbols  represent  experimental  data. 


This  model,  however,  does  not  handle  bipolar  high  injection  conditions  accurately.  It  should  only 
be  used  for  SEU  cases.  Recently,  we  have  developed  a  new,  unified  mobility  approach  to  address 
this  short-coming.  This  paper  is  in  preparation  and  should  appear  shortly. 

3.17  Advanced  Dielectrics 

X-ray  absorption  and  photoemission  spectroscopy,  XAS  and  XPES,  measurements  performed  at 
SSRL  have  been  used  to  obtain  the  conduction  band  edge  and  O-vacancy  defect  states  for  the 
second  and  third  row  TM  oxides,  (i)  Zr02  and  Y203  and  (ii)  Hf02,  respectively.  For  the  first  time 
many-electron  charge-transfer  multiplet  (CTM)  theory  has  been  extended  to  interpretation  of  O  K 
edge  spectra.  A  related  multiplet  approach,  based  on  Tanabe-Sugano  energy  level  diagrams,  has 
also  been  applied,  again  for  the  first  time,  to  a  d2  model  for  O-atom  vacancy  defect  features  in  the 
O  K  pre-edge  regime. 

O  K  spectral  assignments  have  generally  been  based  on  density  of  state  comparisons  that  neglect 
final  state  effects  due  to  O  atom  Is  state  core  holes.  Similarly,  theoretical  studies  for  O-vacancy 
defects  in  Hf02  and  Zr02,  also  based  on  density  functional  theory  (DFT)  have  not  taken  proper 
account  of  the  strongly  correlated  nature  of  TM-atom  d-states.  Resulting  defect  state  energies 
are  inherently  flawed  with  serious  errors  in  interpreting  defect  state  energies  in  spectroscopic 
ellipsometry  (SE).  A  self-consistent  interpretation,  based  on  the  new  approaches  for  O  K  edge  and 
defect  states  demonstrates  how  these  two  studies  can  be  combined  to  correctly  assign  SE  defect 
state  energy  relative  to  conduction  band  edges.  CTM  theory  is  applied  to  M23  and  O  K  edge  XAS 
for  nano-grain  thin  film  Zr02,  and  the  O  K  edge  XAS  of  Hf02.  Tanabe-Sugano  diagrams  have  been 
used  to  assign  spectral  features  in  O  K  edge  spectra  to  O  vacancy  defects  using  the  d2  occupancy 
model. 
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There  has  been  considerable  interest  in  Ge  substrate  channel  regions  for  applications  in  scaled 
complementary  MOS  (CMOS)  devices,  particularly  for  p-MOSFETs.  Many  studies  have  attempted 
to  use  Ge  native  dielectrics  such  as  Ge02,  GeON  and  Ge3N4  as  interfacial  transition  regions  (IRTs) 
between  a  Ge  substrate  and  a  high-K  dielectric.  In  general,  p-  channel  metal  oxide  semiconductor 
field  effect  transistors  (p-MOSFETs)  fabricated  on  Ge  with  Ge-based  ITRs  and  high-K  dielectrics 
such  as  Hf02  and  Zr02  have  yielded  acceptable  levels  of  performance  for  integration  into  scaled 
CMOS  with  interfacial  trap  densities  in  the  mid  1011  cm'2  range. 

In  contrast,  our  previously  reported  study  of  n-MOS  capacitors  (n-MOSCAPS)  with  a  n-Ge/Ge02/ 
Si02  gate  stack  structure  displayed  interfacial  charge  densities  >  mid  1012  cm'2.  More  recently, 
n-MOSFETs  with  Ge3N4  and  GeON  ITRs,  and  high-K  dielectrics  have  been  essentially  written 
off  technology  roadmaps  because  of  high  defect  densities,  >5xl012  cm'2  at  the  Ge/gate  dielectric 
interface.  None  of  these  studies  has  identified  the  source  of  trapping,  e.g.,  a  dangling  bond  or  other 
intrinsic  bonding  defect. 

This  work  focused  on  an  alternative  approach  based  on  spectroscopic  evaluation  of  Ge  surface 
preparation  and  remote  plasma  enhanced  chemical  vapor  deposition  (RPECVD)  of  Ge02.  This 
process  is  qualitatively  similar  to  previous  studies  by  our  group  which  first  identified  a  pathway 
to  producing  device  quality  deposited  Si02.  This  approach  is  based  on  a  processing  protocol 
that  ensures  separate  and  independent  control  of  interface  and  bulk  properties.  The  Ge02  films 
prepared  in  this  way  represent  a  break-through  in  technology  that  will  have  a  significant  impact  on 
radiation  hard  devices.  Continuing  support  for  this  research  by  DTRA  and  the  NSF  will  complete 
the  evaluation  process,  with  test  devices  and  radiation  stress  studies. 

3.18  XAS  studies  of  Gate  Dielectrics 

A.  Background  and  Research  Issues 

Soft  X-ray  (or  vacuum  ultra  violet,  VUV)  absorption  and  photoemission  spectroscopies  (XAS  and 
XPES)  have  been  used  to  study  conduction  band  edge,  and  O-vacancy  defect  states  in  Zr02  and 
Hf02.  XAS  studies  were  performed  at  the  Stanford  Synchrotron  Research  Lightsource  (SSRL). 
Hund’s  rules  apply  for  strongly  correlated  d- state  occupancy  in  O-vacancy  defects  so  that  Tanabe- 
Sugano  diagrams  for  d  to  d’  transitions  can  be  used  to  describe  the  excited  states  of  these  defects. 
This  includes  negative  ion  states  that  act  as  electron  traps.  Density  functional  theory  (DFT)  studies 
of  O-vacancy  defects  in  Hf02  and  Zr02  have  previously  been  applied  to  the  electronic  states  of  these 
defects.  These  calculations  assumed  that  two  electrons  associated  with  the  vacancy  are  distributed 
on  the  transition  metal,  Hf  and  Zr,  d-states  bordering  the  vacancy  in  a  way  that  was  not  consistent 
with  Hund’s  rule  occupancy,  yielding  doubly  occupied  ground  and  negative  ion  states.  This  has 
lead  to  quantitative  errors  in  interpreting  defect  state  energies  in  spectroscopic  ellipsometry  (SE) 
studies,  and  in  applying  results  to  electrical  measurements,  as  well.  Tanabe-Sugano  energy  level 
diagrams  in  this  report  in  Fig  34  are  applied  to  an  equivalent  d2  state  model  for  O-vacancy  defects 
in  the  O  K  pre-edge  regime  providing  a  correct  assignment  for  SE  features  and  their  relationship  to 
electronically  active  traps. 

B.  Conduction  Band  Edge  States 

O-vacancy  defect  states  are  introduced  into  the  forbidden  energy  gap,  and  as  such  their  energies 
relative  to  conduction  and  valence  band  edges  are  important  for  identifying  the  energies  of 
electronically  active  hole  and  electron  traps  relative  to  these  band  edges.  Doubly  occupied  defect 
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conduction  band  edge 
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singlet  states  triplet  states 


valence  band  edge 

Fig.  34.  Energy  level  diagram  for  O-vacancy  defects,  constructed  from  Tanabe-Sugano  diagram  for  a  d2 

occupancy. 


states  have  previously  been  identified  at  the  valence  band  edge  for  Ti02,  Zr02  and  Hf02  using 
XPES  spectroscopy.  The  excited  states  of  these  defects  have  been  detected  by  second  derivative 
analysis  of  pre-edge  of  O  K  edge  XAS  spectra,  but  have  not  as  yet  been  addressed  quantitatively 
using  a  many-electron  theory  approach. 

Charge  transfer  multiplet  (CTM)  approach  is  applied  for  first  time  applied  to  O  K  edge  spectra 
for  nano-grain  thin  film  Zr02  and  Hf02,  the  most  important  aspect  of  this  illustrated  in  Fig.  35  for 
Zr02.  This  approach  takes  into  account  a  localization  of  the  O  K  edge  core  hole,  and  a  coherent 
process  in  which  this  is  filed  from  O-atom  2p  jt  states  in  the  valence  band.  For  Zr02,  this  coherent 
excitation  process  yields  a  final  O  2p54d1  state  that  has  the  same  symmetry  as  the  M23  final  states. 
Consider  first  M23  feature  spin-orbit  splitting.  For  the  Zr  M23  spectrum  this  is  the  Zr  3p  spin  orbit 
splitting  of -13  eV.  For  the  O  K  edge  states,  the  corresponding  spitting  of -10  eV  is  between  O  2p 
a-bonding  and  Jt-bonding/non-bonding  valence  band  states  as  indicated  in  Fig.  33.  For  purposes 
of  this  section  of  the  report  that  focuses  on  O-atom  vacancies,  this  approach  to  conduction  band 
edge  states  provides  a  way  to  correlate  the  energies  of  excited  O-atom  vacancy  and  negative  ion 
states  with  the  band  edge  Zr  4d  and  Hf  5d  Eg  states,  and  as  such  with  effective  band  edge  energy 
for  alignment  of  states  Zr(Hf)02-Si  interfaces.  This  aspect  of  CTMs  for  O  K  edge  features  and  its 
relationship  to  SE  and  electrical  measurements  will  be  addressed  in  more  detail  elsewhere. 

C.  Spectroscopic  Studies  of  O-vacancy  States 

Figure  36  includes  O  K  edge  -  pre-edge  spectra  for  5  nm  thick  films  of  Hf02,  and  cubic  Hafnia  with 
Y  atoms  substituted  on~  15%  of  the  Hf  sites.  The  +3  charge  of  each  Y  atom  is  compensated  by  an 
O-atom  vacancy  that  is  part  of  the  host  crystal  electronic  structure  and  does  not  contribute  localized 
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Fig.  35.  M  and  O  K  traces  and  2nd  derivative  for  nano-grain  Zr09.  Symmetries,  and  spin- 
orbit  and  O  2p  splittings  are  shown. 


defect  states  with  d-d’  transitions,  including  negative  ion  states.  The  ground  state  symmetry  for 
the  d2  configuration  in  8 -fold  cubic  symmetry  is  the  same  as  the  d8  symmetry  for  6-fold  octahedral 
symmetry.  This  ground  state  is  an  “s-like”  or  one-dimensional  triplet  state,  3A1ct.  Based  on  the 
appropriate  Tanabe-Sugano  diagram  energy  levels  in  Fig.  34,  the  allowed  optical  transitions  are  to 
“p-like”,  3Tlg  and  3T2g  states  as  indicated  in  the  diagram.  The  negative  ion  states  are  *A  and  !Tlg, 
and  these  can  be  occupied  by  electron  injection  from  the  conduction  band  and  contribute  to  trap- 
assisted  tunneling  (TAT)  observed  in  both  Hf02  and  Zr02  films  that  are  thicker  than  3-5  to  4  nm  [1] . 
TAT  is  suppressed  in  films  thinner  than  ~  2  nm. 

The  unit  cell,  and  extended  unit-cell  symmetry  in  the  tetragonal  Hf02  is  higher  than  that  in  the 
compositionally  averaged  alloy  structure  of  cubic  Hafnia,  and  this  is  reflected  in  the  differences 
in  symmetry  of  the  respective  3Tlg(P)  structures  in  Fig.  36.  The  three  features  in  each  symmetry 
designated  d-d’  transition  is  a  manifestation  of  the  3-fold  degeneracy. 

The  symmetry  of  the  same  three  features  in  monoclinic  Hf02  with  a  seven-fold  coordinated  Hf  is 
significantly  reduced.  There  is  an  increase  in  the  number  of  features  for  d-d’  transitions  to  14  from  9 
in  tetragonal  Hf02.  This  is  due  to  singlet-triplet  mixing,  and  is  present  for  negative  ion  states  as  well 

D.  Summary  of  Significant  Results 

(a)  The  conduction  band  states  observed  in  O  K  edge  X-ray  spectroscopy  have  the  same  relative 
energies  as  those  observed  in  SE,  and  cannot  be  interpreted  as  a  joint  density  of  states  as 
originally  proposed. 
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Fig.  36.  O  K  pre-edge  spectra  for  d-d’  dipole  allowed  transitions  and  dipole  forbidden  negative  ion  states:  (right)  in 

tetragonal  Hf02,  and  (left)  cubic  hafnia. 


(b)  The  Zr  M2  3  spectrum,  and  O  K  edge  spectrum  of  both  Zr02,  as  well  as  the  L2  3  spectrum  and 
O  K  edge  of  Ti02  must  be  analyzed  in  terms  CTM  theory,  contrary  to  the  previous  model 
proposed. 

(c)  In  this  report  O-vacancy  defects  are  described  by  a  theory  that  takes  proper  account  of  Hund’s 
rule  occupancy,  and  therefore  DFT  theories  used  previously  are  not  valid.  The  energy  levels 
of  the  d-  dD  transitions  from  the  d2  ground  state,  and  negative  ion  states  as  well  are  obtained 
by  application  of  Tanabe-Sugano  energy  level  diagrams. 

(d)  Differences  in  the  band  edge  symmetry  between  non-crystalline  Si02  and  Hf02  account  for 
absence  of  trap-assisted  tunneling. 

(e)  Gd203  and  Lu203  display  similar  pre-edge  O-vacancy  defects  in  their  O  K  edge  spectra.  The 
as-deposited  Lu203  spectra  is  included  in  Fig.  37 

3.19  Conduction  Band  Edge  Excitonic  States  and  Electron  Trapping 
in  Dielectrics 

A.  Introduction 

The  performance  and  reliability  of  a  gate  dielectric,  (i)  non-crystalline  Si02,  Si3N4,  a  Si  Oxynitride 
alloy,  or  (ii)  a  high-x  transition  metal  oxide  are  determined  by  band  edge  intrinsic  bonding 
states,  intrinsic  bonding  defects  and  macroscopic  strain.  This  final  report  identifies  significant 
differences  between  band  edge  excitonic  states  in  (i)  Si02,  and  (ii)  Si3N4,  and  Si  oxynitride  alloys, 
(Si3N4,X(Si02)l-X.  These  provide  an  explanation  for  differences  in  trapping  and  trap-assisted 
(TAT).  Differences  between  symmetries  of  band  edge  states  in  Si02  and  transition  metal  oxides 
also  explain  differences  in  TAT  injection  from  negatively  biased  n-type  Si  substrates. 

A  connection  between  (i)  strain-reducing  medium  range  order  (MRO),  and  (ii)  nano-scale 
separation  into  hard-soft  mixtures  in  non-crystalline  Si02  is  explained  by  many-electron  theory 
applied  to  Si-atom  “d-like  states”.  These  states  participate  in  and  O  pjt  to  Si  djt  back-donation  of 
electrons,  yielding  shortened  Si-O  bond-lengths  and  MRO  atom-pair  correlations.  These  bonding 
interactions  specifically  identify  a  coherence  length  of  ~1  nm  associated  with  hard  6-member  rings/ 
ring  segments  encapsulated  by  compliant  or  5-  and  7-member  rings/ring  segments  into  a  nano-grain 
hard-soft  nano- structure,  reducing  macroscopic  strain  and  giving  Si02  unique  reliability  properties. 
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Fig.  37.  O-vancancy  states  in  the  pre-edge  spectra  of  nano-grain  Lu203. 

B.  Experimental  Procedures 

Five  nm  thick  films  of  Si02,  Si3N4,  and  Si  oxynitride  alloys  were  remote  plasma  deposited  on 
nitrided  superficially  oxidized  Si(001).  These  films  were  annealed  in  Ar  at  a  temperature  of  ~950°C. 
Si  L23,  and  O  and  N  K  edge  spectra  were  obtained  by  X-ray  absorption  spectroscopy  (XAS)  at  the 
Stanford  Synchrotron  Research  Lightsource  (SSRL). 

C.  Experimental  Results 

Figure  38  compares  the  Si  L23  spectra  of  non-crystalline  Si02  and  crystalline  Si.  The  band  edge 
states  in  these  two  materials  that  comprise  the  Si-Si02  interface  are  “s-like”  non-degenerate  2Alg 
states.  Figure  39  is  the  2nd  derivative  spectrum  of  the  band  edge  states  in  non-crystalline,  thin  film 
plasma  deposited  Si02  that  has  been  annealed  at  ~950°C.  There  is  a  one-to-one  correspondence 
between  the  energy  differences  in  eV  units  of  band  edge  states  as  determined  by  visible  and  vacuum 
UV  spectroscopies,  and  in  the  O  K  edge  spectra.  An  X-ray  energy  of  529.25±0.1  eV  corresponds 
to  an  energy  of  8.9  eV  for  the  band-gap  of  non-crystalline  Si02.  In  order  of  decreasing  X-ray 
energy,  and  as  marked  in  Fig.  39,  the  features  in  the  O  pre-edge  derivative  spectra  correspond 
to:  i)  the  band-gap,  Eg,  ii)  two  bound  excitons,  E2  and  Ex,  iii)  negative  ion  states  between  529.1 
and  527.4  eV,  and  iv)  four  intra-d  state  transitions  for  the  O-atom  vacancy  defect  represented 
by  a  high-spin  d2  state.  The  symmetries  of  these  states  are  identified  using  the  Tanabe-Sugano 
diagrams.  The  combination  of  the  localized  2Alg  excitonic  states  at  the  Si02  band  edge,  and  the 
xAlg  and  1T2g  symmetries  of  unoccupied  negative  ion  states  of  O-vacancy  defects  does  not  support 
dipole  allowed  fast  radiative  decay  into  negative  ion  states.  This  explains  why  TAT  processes  have 
not  been  reported  for  negatively  biased  n-type  Si  in  Si-Si02  gate  stack  structures.  In  contrast,  the 
symmetries  of  band  edge  Zr02  and  Hf02  states  allows  tunneling  injection  to  the  negative  ion  states 
of  O-vacancy  defects.  This  explains  the  TAT  processes  initiated  by  electron  injection  from  n-type 
Si  substrates  that  have  been  reported  for  Zr02  and  Hf02  MOS  structures. 
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Fig.  38.  L23  spectra  for  c-Si  and  nc-Si02. 
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Fig.  39.  Pre-edge  O  K  spectra  for  nc  thin  film  Si02. 


Figures  40  and  4 1  indicate  respectively,  (i)  O  K  band  edge ,  and  pre-band  edge  states  of  non-crystalline 
Si3N4„  and  expanded  scale  plots  of  the  negative  ion  states  of  (ii)  Si3N4,and  (ii)  a  (Si3N4,0.5(SiO2)0.5 
Si  oxy-nitride  alloy  that  has  been  used  as  alternative  gate  dielectric  in  the  first  eight  to  ten  years  of 
the  21st  century.  There  is  a  significant  qualitative  difference  between  Fig. 39  for  Si02  and  Fig.  40 
for  Si3N4,  that  is  related  to  the  difference  in  the  number  of  p  electrons  in  the  ground  states  of  O-  and 
N-atoms,  four  for  O,  and  five  for  N.  This  correlates  with  the  singly  occupied  2pjt  state  that  gives 
rise  to  sharp  spectral  feature  at  400  eV  in  Si3N4„  and  Si  oxynitride  alloys  including  compositions 
both  Si02-rich  and  Si3N4,-rich  as  well.  It  is  significant  to  note  that  the  final  state  for  the  N-atom  “pjt 
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to  pjt*”  transition  is  at  an  energy  that  is  between  the  band  edge  excitonic  states  of  Si3N4,  and  Si 
oxynitride  alloys,  and  the  negative  ion  states  associated  with  N-vacancy  defects  in  these  dielectrics. 

This  difference  in  the  ordering  of  electronic  states  manifests  itself  in  Si3N4  dielectric  thin  films  by 
promoting  TAT  and/or  Poole-Frenkel  transport.  The  final  state  symmetry  of  the  N-atom  pjt  to  pjt* 


Fig.  40.  N  K  edge  and  pre-edge  spectra  -  Si3N4. 


X-ray  photon  energy  (eV) 


Fig.  41.  Expaded  x-axis  for  N  K  pre-edge  spectra  for  Si3N4  and  Si  oxynitride  alloys. 


MURI  FINAL  REPORT  •  2005  -  2010  •  VANDERBILT  UNIVERSITY 

PAGE  32 


transition  provides  a  transport  channel  between  the  “s-like”  symmetry  of  the  band  edge  states,  and 
the  even  symmetry  of  the  negative  ion  states.  In  addition,  it  is  likely  that  this  state  plays  a  role  in 
the  Negative  Bias  Temperature  Instability,  NBTI,  that  is  much  stronger  in  N-containing  dielectrics 
than  Si02. 

Finally,  Fig.  42  presents  the  O  K  pre-edge  spectra  of  (SiO2)03(HfO2)03(Si3N4,)04,  or  HfSiON334. 
The  three  relative  sharp  features  between  530.2  and  531 .5  eV  and  associated  with  band  edge  states 
of  SiON  alloys  as  indicated  in  Figs.  39  and  40.  The  lower  eV  features  in  this  figure,  are  associated 
with  O- vacancies  with  Hf  neighbors,  as  represented  by  a  d2  state.  The  symmetry  designations,  and 
spectral  width  of  these  states  are  very  nearly  the  same  as  the  defect  states  for  O-vacancies  in  Hf02 
as  displayed  in  Fig.  37.  The  SiON  band  edge  states  overlap  the  negative  ion  states  of  the  O- vacancy 
defects.  They  are  on  different  atoms,  and  the  density  of  the  Hf02  negative  ion  states  is  at  least 
three  orders  of  magnitude  smaller  than  the  SiON  states.  This  assignment  is  supported  by  the  results 
in  Fig.  43  which  indicate  that  the  conduction  band  edge  states  in  the  HfSiON334  alloy  are  those 
displayed  in  Fig.  42  and  designated  there  as  Si  band-edge  “s-like”,  one  dimensional  Alg  states.  This 
overlap  explains  why  this  current- voltage  studies  of  gate  stacks  with  this  alloy  do  display  TAT,  and 
why  the  response  to  X-ray  stress  is  essentially  the  same  as  that  of  Si02,  in  particular  there  is  no 
build  up  of  negative  space  charge  since  the  negative  ion  states  are  not  accessible  due  to  overlap  in 
energy  that  has  been  used  to  explain  the  spectral  data  in  Fig.  42. 

D.  Summary  of  Significant  Results 

The  different  symmetries  of  2Alg  excitonic  states  at  the  Si02  band  edge,  and  O-vacancy  2A1ct  and 
1T2g  unoccupied  negative  ion  states  effectively  prevents  injection  into  these  negative  ion  states 
suppressing  trap-assisted  tunneling  (TAT).  The  combination  of  symmetries  of  Zr02and  Hf02  band 
edge  states  allows  TAT  for  electron  injection  from  n-type  Si  substrates.  The  differences  in  O  and 
N-atom  p-states  (4  compared  with  3)  results  in  a  N-atom  pjt  to  pjt*  transition  sandwiched  between 
the  “s-like”  symmetry  of  band  edge  states,  and  the  even  symmetry  of  N-vacancy  negative  ion 
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Fig.  42.  O  K  pre-edge  spectra  -  SiOn334. 
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Fig.  43.  O  K  edge  and  pre-edge  SiOn 


states.  This  provides  a  “pathway”  for  TAT  processes,  well  known  in  Si3N4,  and  extending  to  Si 
oxynitride  alloys,  e.g.,  HfSiON334,  as  well. 

3.20  Low  Defect  Density  Remote  PECVD  Ge02 

A.  Introduction 

Contrary  to  the  conventional  wisdom,  there  are  not  a  larger  number  of  different  defect  states  within 
the  band-gap  of  Si02.  Instead  the  energies  of  absorption  and  luminescent  transitions  have  been 
assigned  to  transition  energies  of  a  d2  high  spin  state  localized  on  one  the  two  Si  atoms  associated 
with  the  formation  of  an  O-atom  vacancy.  When  one  of  these  electrons  is  removed  by  exposure 
to  UV  radiation,  y-rays,  x-rays,  and/or  high-energy  electrons,  the  singly  occupied  state  is  the 
well-known  E’  center.  The  O-atoms  generated  in  the  defect  activation  process,  3/20-Si-0-Si-03/2  = 
3/2O-Si[0e]+  +  -[2e]Si-03/2  +  O0,  contribute  to  hole-related  states,  O0  +  3/2Q-Si-0-Si-03/2  =3/20-Si-(0- 
[lh]-0)°  +  0[le]Si-O3/2,  identified  by  electron  spin  resonance  (ESR). 

A  similar  decomposition  reaction  to  the  one  quoted  for  Si02,  occurs  for  Ge02  in  the  quartz  (4-fold 
Ge  and  2-fold  O)  bonding  arrangement,  that  present  in  non-crystalline  bulk  Ge02  glass.  Based  on 
bond-enthalpies,  the  barrier  for  the  Ge02  decomposition  reaction  is  reduced,  and  hence  the  defect 
density  of  bulk  glasses  is  higher,  and  the  surface  reaction  rate  for  oxidation  of  bulk  Ge  is  higher. 
I  have  performed  experiments  on  this  glass  in  1976,  and  it  did  indeed  have  a  yellow  color  in 
transmission.  All  studies  have  to  date  have  yield  a  band-gap  of  -5.5  eV,  consistent  with  this  yellow 
color.  The  same  band  gap  has  been  reported  for  Ge02  formed  by  oxidation  of  bulk  Ge.  The  issue 
is  how  to  produce  Ge02  with  a  lower  defect  density  than  that  of  the  bulk  glass,  and  thermally- 
oxidized  Ge  surface.  The  answer  is  simple  and  straight-forward:  our  group  uses  a  non-equilibrium 
low  temperature  plasma-assisted  process,  and  then  determine  how  the  defects  created  by  the 
decomposition  reaction  change  with  post-deposition  annealing  temperature  in  inert  ambients.  This 
is  exactly  what  we  have  done,  and  it  works  very  well! !  It  must  be  pointed  out  that  the  rutile  6-fold 
coordinated  Ge  form  of  non-crystalline  (and  crystalline  as  well)  Ge02  is  highly  soluble  in  water 
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(hygroscopic),  while  the  quartz  4-fold  coordinated  form  is  not.  This  is  important  for  technological 
applications,  and  is  validated  by  our  experience  with  measurements  made  on  films  of  the  quartz¬ 
like  non-crystalline  form  at  two  SSRL  runs  spaced  by  more  than  4  weeks. 

3.21  Summary  of  Preliminary  Results  for  Plasma  Deposited  Ge02 
Thin  Films 

A.  Introduction 

Most  of  the  results  below  come  from  XAS  measurements  made  at  SSRL.  These  are  combined  with 
depth  resolved  CLS  by  the  Brillson  group  at  The  Ohio  State  University.  Finally,  by  a  model  for 
the  Ge-Ge02  interface  is  proposed  based  on  our  measurements  in  combination  with  comparison 
between  spectroscopic  ellipsometry  measurements  of  Aspnes  and  co-workers  on  crystalline  Si  and 
Ge,  and  theoretical/experimental  values  of  band  offset  energies  between  Si  and  Si02,  and  Ge-Ge02 
interfacial  regions. 

B.  RPECVD  Growth  of  Ge02  Thin  Films 

The  thin  films  of  Ge02  were  deposited  by  remote  plasma  enhanced  chemical  vapor  deposition  onto 

(i)  remote  plasma  nitrided  (RPAN)  surfaces  of  crystalline  Si  and  Ge  to  prevent  oxidation  during 
film  growth,  and  (ii)  wet-chemical  cleaned  Ge  surfaces,  using  a  process  that  gives  minimal  Ge02 
native  oxide  growth.  Films  were  subsequently  annealed  (1  minute  RTA’s)  in  Ar  at  temperatures  of 
400°C,  600°C  and  700°C. 

C.  Summary  of  XAS  and  DRCLS  Results 

(i)  empty  conduction  band  states  (nitrided  Si  substrates)  -  Fig.  44  -  (a)  separation  in  “s-like”  and 
“d-like”  states  a  higher  energies;  (b)  similar  AE  with  respect  to  band  gaps  of  9.5  eV  for  Ge02 
and  8.9  eV  for  Si02;  (c)  2nd  derivative  spectra  (not  shown)  indicate  two  bound  excitons  - 
spacing  ~  0.2  eV;  (d)  “d-like”  sharper  in  Ge02. 

(ii)  empty  conduction  band  states  (cleaned  Ge  substrates)  -  Fig.  45  -  (a)  as  deposited  Ge02  films, 
and  films  annealed  at  700°C  in  Ar  display  spectra  qualitatively  different  than  those  in  Fig.  44 
with  features  that  are  broader,  and  less  well  defined,  whereas  (b)  films  annealed  at  400°C  and 
600°C  display  sharper  features  similar  to  those  in  Fig.  44;  these  spectra  are  higher  energy 
because  the  2p  states  of  the  Si  conduction  band  are  at  lower  energy  than  the  3p  states  of  Ge 
because  of  repulsion  effects  associated  with  the  occupied  Ge  3d  shallow  core  states. 

(iii)  Figs.  46  and  47  compare  respectively  (a)  the  spectral  dependence  of  the  DRCLS  spectrum, 
and  (b)  the  integrated  response.  The  spectral  dependence  and  relative  integrated  response  are 
different  for  the  lease  defect  processing  —  deposition  on  clean  Ge  followed  by  annealing  at 
400  C  and  600C.  The  spectral  response  is  peaked  at  lower  eV  significantly  higher  emission 
in  the  three  more  defective  samples,  and  the  integrated  response  is  a  factor  of  10  or  higher  in 
the  defective  samples  as  well. 

D.  Band  Off-sets  at  “ldeal”(No  Interfacial  Transition  Regions)”  Ge-Ge02  and 
Si-Ge02  Interfaces 

Based  on  the  valence  and  offset  energies  for  Si-Si02,  and  Ge-Ge02,  conduction  and  valence  band 
energies  relative  to  Si  and  Ge  have  be  estimated.  The  conduction  band  energy  is  about  1.5  eV 
higher  for  Ge-Ge02,  and  the  valence  band  energy  is  about  1  eV  smaller.  This  is  indicated  in  Fig.  48. 
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Fig.  44.  O  K  edge  spectra:  Si02  and  GeOr 
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substrate  temperature  (degrees  C) 


Fig.  47.  Integrated  DRCLS  spectra  of  Ge02  films  on  Ge  substrate. 
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Fig.  48.  c-Si  and  c-Ge  band  edge  alighments  with  Si02  and  GeOz,  respectively. 

3.22  Total  Ionizing  Dose  and  Single  Event  Effects  in  Strained  Si 
Technologies 

Uniaxial  strained- silicon  (Si)  and  high-k  gate  dielectrics  are  key  technologies  used  to  enhance 
transistor  performance  for  sub  100-nm  logic  technology  nodes.  Uniaxial  mechanical  stress 
improves  device  characteristics  such  as  mobility  and  gate  tunneling  current,  with  minimal  stress- 
induced  threshold- voltage  shifts.  Strained- silicon  technology  increases  drive  current  using 
additional  process  features  such  as  tensile  SiN  capping  layers  to  create  tensile  stress  in  nMOSFETs 
and  compressive  SiN  capping  layers  and  SiGe  source/drain  regions  to  create  compressive  stress 
in  pMOSFETs.  Channel  stresses  of  up  to  1.5  GPa  can  be  produced  using  these  stressors.  These 
strained-Si  processes  also  induce  additional  stress  (<  1  GPa)  such  as  uniaxial  or  biaxial  stress  in  the 
source/drain  region.  The  combined  effects  of  the  stressors  cause  changes  in  carrier  mobility  in  the 
source/drain  and  bulk  (substrate)  regions,  as  well  as  in  the  channel. 

Radiation-hardened  devices  technology  for  space  and  military  electronics  market  has  strongly 
influenced  by  MOS  devices  technology  for  commercial  electronics  market.  Since  the  cost  of  making 
radiation  hardened  devices  becomes  more  expensive,  people  who  are  working  for  electronics  in  the 
radiation  environment  make  their  effort  to  reduce  the  cost  using  commercial  technology.  However, 
no  systematic  study  how  strained  Si  technology  responds  to  radiation  environment  has  been  done 
yet.  This  work  investigates  the  effect  of  strained  mechanical  stress  on  total  ionizing  dose  effects  on 
strained  nMOSFETs  and  laser  induced  current  transients  in  strained  Si  diodes  experimentally  and 
theoretically. 

A.  Total  Ionizing  Dose  Effects  on  Strained  Hf02-based  nMOSFETs 

High-k  gate  dielectrics  have  been  implemented  to  reduce  transistor  gate  leakage  current  in  the  45 
nm  CMOS  technology  node.  Hafnium-based  dielectrics  with  a  relative  dielectric  constant  of  ^15  to 
26  have  emerged  as  the  materials  of  choice  for  high-k  gate  dielectrics.  Although  radiation  damage 
of  Hf02-based  MOS  devices  has  been  studied  in  recent  years,  the  effects  of  uniaxial  mechanical 
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stress  on  the  ionizing  radiation  response  of  Hf02-based  MOS  devices  have  not  been  reported. 
Hole  trapping  is  observed  to  be  dominant  in  Hf02  [17,  18]  dielectric  layers,  similar  to  Si02.  The 
effects  of  mechanical  stress  on  the  ionizing  radiation  response  of  Si02-based  MOS  devices  have 
been  reported,  but  the  mechanical  stress  (<  4  MPa)  produced  in  these  studies  by  changing  the  gate 
electrode  thickness  is  much  smaller  than  that  used  in  strained-Si  technology  (~  1  GPa). 

In  high-k  transistors,  remote  Coulomb  scattering  (RCS)  is  one  of  the  main  factors  limiting  the 
mobility.  In  commercial  devices,  most  of  the  RCS  effect  comes  from  fixed  charges  generated  by 
the  fabrication  process.  Similarly,  radiation-induced  charges  may  increase  RCS.  Owing  to  the 
ubiquitous  use  of  high-k  based  strained-Si  technology  in  the  future,  it  is  important  to  address  the 
radiation  response  of  these  devices  as  a  function  of  stress.  In  this  work,  we  investigate  the  effects  of 
uniaxial  stress  on  the  radiation-induced  threshold  voltage  shifts  and  mobility  degradation  in  Hf02- 
based  nMOSFETs  using  controlled  external  mechanical  stress. 

Experimental  Details 

Radiation-induced  threshold  voltage  shifts  and  mobility  degradation  in  mechanically  stressed 
Hf02-based  nMOSFETs  are  extracted  from  drain  current-gate  voltage  (ID-VGS)  characteristics. 
The  samples  used  in  this  study  are  TiN/Hf02  gate  stack  nMOSFETs  on  Si  (001)  wafers  with  <110> 
channel  direction.  The  high-k  gate  dielectric  is  7.5  nm  in  physical  thickness.  The  thickness  of  the 
SiOx  interlayer  is  1  nm.  The  effective  oxide  thickness  (EOT)  of  these  transistors  is  2.3  nm.  The 
aspect  ratio  (W/L)  of  the  transistors  used  in  this  investigation  is  10  //m/3  //m.  Strain  in  the  devices 
is  produced  by  applying  uniaxial  mechanical  stress  using  a  four  point  bending  jig,  as  shown  in  Fig. 
49.  The  devices  are  irradiated  in  an  ARACOR  10-keV  X-ray  irradiation  system  under  different 
values  of  stress  (100  MPa  and  200  MPa  tensile,  no  stress,  and  200  MPa  compressive),  at  -2  V  gate 
bias,  with  the  source,  drain,  and  body  grounded. 


Table  1.  Total  ionizing  dose  (in  Si02)  experimental  matrix. 
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Although  the  maximum  applied  stress  (~200  MPa)  in  this  investigation  is  about  20%  of  that 
produced  by  process-induced  stressors  (~  lGPa),  the  experiments  provide  insight  into  the  effects 
of  stress  on  radiation-induced  threshold  voltage  shifts;  this  approach  is  analogous  to  previous  work 
describing  the  effects  of  stress  on  unirradiated  MOS  threshold  voltages.  Samples  are  irradiated  to  a 
cumulative  dose  of  5  Mrad(Si02)  at  a  dose  rate  of  31.5  krad(Si02)/min;  the  experimental  matrix  is 
shown  in  Table  1 .  Post  irradiation  ID-  VGS  curves  are  measured  using  an  Agilent  4156  semiconductor 
parameter  analyzer.  Threshold  voltage  shifts  (AVT)  are  monitored  as  a  function  of  radiation  dose  for 
different  applied  uniaxial  mechanical  stresses.  Threshold  voltages  are  extracted  using  a  constant- 
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current  (CC)  method.  Pre-  and  post-irradiation  subthreshold  behavior  is  also  studied.  -2  V  gate  bias 
stress  experiments  without  irradiation  are  performed  to  separate  the  contributions  of  bias-induced 
charge  trapping/creation  in  these  gate  dielectrics  from  the  radiation-induced  shifts.  Electron 
mobility,  /i  =IdLJWCqx  (Vgs  -  Vr)  VDS  is  extracted  as  a  function  of  gate  over-drive  voltage  (VGS-VT) 
at  low  drain  voltage  ( VDS=  0.1  V),  where  Lg  is  the  channel  length,  W is  the  channel  width,  Cox  is  the 
gate  oxide  capacitance  per  unit  area,  and  VT  is  the  threshold  voltage. 

Threshold  voltage  instability  is  one  of  the  main  issues  in  high-k  devices.  In  this  work,  threshold 
voltage  shifts  are  measured  after  being  stabilized  by  sweeping  the  ID~VGS  curves  before  and 
after  irradiation.  For  the  samples  used  in  this  investigation,  the  threshold  voltage  shift  stabilizes 
generally  after  the  second  ID~VGS  curve  after  irradiation  or  bias  stress.  The  stabilization  occurs 
because  charges  in  shallow  hole  trap  sites  in  the  SiOxor  Hf02  may  detrap,  or  holes  in  deeper  trap 


Fig.  49.  Experiment  setup  for  measuring  total  ionizing  dose  (TID)  effects 
on  a  mechanically  strained  nMOSFET,  and  cross  section  of  gate  stack  of 


sites  near  the  interface  may  recombine  with  injected  electrons  during  the  first  I D-VGS  sweep  under 
a  given  irradiation  or  bias  condition.  The  effect  of  hot  carrier  injection  on  the  threshold  voltage 
shift  during  multiple  I-V  sweeps  is  expected  to  be  negligible  since  VDS  is  very  low  (~100  mV)  and 
there  is  no  abrupt  increase  in  gate  leakage.  Previous  work  by  Dixit  et  al.  focused  on  the  threshold 
voltage  shift  of  the  first  I-V  curve,  which  includes  the  effects  of  transient  shifts.  This  work  focuses 
instead  on  the  changes  in  radiation-induced  charge  trapping  under  uniaxial  mechanical  stress  after 
the  threshold  voltage  stabilizes. 


Fig.  50.  Semilog  and  linear  plots  of  the  ID-VGS 
characteristics  as  a  function  of  the  accumulated  x-ray 
dose  (in  Si02)  under  tensile  stress  of  200  MPa. 
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Results  and  Discussion 

Hole  trapping  in  the  gate  oxide  is  the  dominant  radiation-induced  charge  for  these  devices  under 
the  irradiation  conditions,  as  seen  by  the  decrease  in  the  threshold  voltage  as  shown  in  Fig.  50. 
This  agrees  with  previous  results  on  Hf02-based  nMOSFETs.  Threshold  voltage  shifts  (AVj)  can 
be  caused  by  interface  trapped  charge  (A  V.)  and  oxide  trapped  charges  (A  VJ .  Since  there  is  no 
significant  change  observed  in  the  subthreshold  slope  in  Figure  50,  the  threshold  voltage  shifts 
are  caused  mainly  by  an  increase  in  AVot.  Transistors  irradiated  under  other  stress  conditions 
(compressive  (200  MPa),  no  stress,  and  tensile  (100  MPa))  are  also  dominated  by  positive  charge 
trapping. 

Radiation  Induced  Threshold  Voltage  Shifts  Under  Mechanical  Stress 

The  effects  of  applied  mechanical  stress  on  charge  trapping  are  characterized  by  monitoring 
threshold  voltage  shifts  at  each  radiation  dose.  The  tensile  stress  effect  is  shown  in  Fig.  5 1 .  Increasing 
tensile  stress  results  in  less  threshold  voltage  shift  at  each  dose  level  than  that  measured  for  devices 
irradiated  with  no  stress  applied.  Smaller  threshold- voltage  shifts  are  observed  after  applying  only 
bias  for  a  total  time  of  2.5  h,  which  is  equivalent  to  the  time  required  for  5  Mrad(Si02)  irradiation. 
Tensile  stress  also  reduces  the  threshold  voltage  shift  resulting  only  from  the  bias.  The  data  points 
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Fig.  51.  Threshold  voltage  shifts  (AVT)  due  to  —2  V  gate 
bias  +  radiation  and  —2  V  gate  bias  without  radiation 
under  no  stress,  tensile  stress  of  100  MPa,  and  200  MPa 
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Fig.  52.  Threshold  voltage  shifts  (AVT)  due  to  —2  V  gate  bias 
+  radiation  and  —2  V  gate  bias  without  radiation  under  no 
stress,  compressive  stress  of  200  MPa  conditions. 
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Mechanical  Stress  (MPa) 

Fig.  53.  Threshold  voltage  shifts  (AVT)  vs.  mechanical  stress 
after  5Mrad  (Si02)  and  2.5  h  under  -2  V  gate  bias  (positive 
(+):  tensile,  negative  (-):  compressive) 

are  the  average  threshold  voltage  shifts  at  each  radiation  dose  or  bias  time.  The  error  bars  in  the 
data  points  represent  the  standard  deviation  in  the  data  at  each  dose  and  stress  level.  Fig.  52  shows 
a  similar  trend  for  the  threshold  voltage  shifts  under  200  MPa  of  compressive  stress. 

In  contrast  to  the  reported  stress  dependence  of  the  Si02  nMOSFET  threshold  voltage  shift  under 
irradiation,  both  applied  tensile  and  compressive  uniaxial  stress  reduce  the  threshold  voltage  shifts 
in  devices  with  Hf02  and  SiOx  dielectrics  in  Fig.  53  that  are  either  irradiated  under  bias,  or  subjected 
only  to  bias  stress  without  irradiation,  for  comparison.  A  possible  explanation  of  these  results  is  that 
compressive  and  tensile  uniaxial  mechanical  stress  both  lower  the  hole  trap  energy  level  in  Hf02 
and/or  SiOx,  reducing  hole  trapped  charges.  In  recent  work,  trap-assisted  gate  tunneling  current 
in  high-k  MOS  capacitors  increased  under  both  compressive  and  tensile  stress,  suggesting  that 
the  hole  trap  energy  distribution  may  be  shifted  to  lower  average  values  by  uniaxial  mechanical 
stress.  In  addition,  an  increase  in  trap  assisted  tunneling  in  thin  (2.5  nm)  Si02  results  from  reduced 


(a) 


Fig.  54.  (a)  Charge  detrapping/ 
neutralization  model;  (b)  multiple  trapping¬ 
detrapping  hole  transport  model. 
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trap  activation  energy  due  to  both  compressive  and  tensile  stress.  Uniaxial  mechanical  stress  may 
change  the  trap  energy  levels  by  changing  bond  lengths  and  angles  in  Hf02  and  SiOx.  This  is 
consistent  with  previous  works  on  the  oxygen  vacancy  defect  that  show  that  the  trap  microstructure 
and  energy  levels  can  be  changed  by  stretching  the  Si-Si  bonds  and/or  changing  the  bond  angles. 


We  consider  two  possible  reasons  why  lowering  hole  trap  energy  levels  may  reduce  charge  trapping 
in  Hf02  and/or  SiOx.  First,  uniaxial  mechanical  stress  may  enhance  the  detrapping  of  holes  in 
shallow  trap  sites  or  the  neutralization  in  deep  trap  sites  by  electron  injection,  as  illustrated  in  Fig. 
54(a).  This  reduction  in  hole  trap  energy  increases  the  probability  that  these  defects  can  emit  a 
trapped  hole  or  capture  an  electron  to  compensate  a  nearby  trapped  hole.  Trapped  holes  in  deep  trap 
sites  can  be  neutralized  by  capturing  electrons. 

Second,  the  effective  hole  mobility  in  the  gate  dielectrics  along  the  <00 1>  direction  under  uniaxial 
mechanical  stress  may  be  increased  by  reducing  the  average  trap  energy  level.  Hole  transport  in  thin 
(~10  nm)  high-k  oxides  can  be  described  by  a  multiple  trapping  model.  This  motion  is  illustrated 
schematically  in  Fig .54(b).  Reduced  trap  energy  levels  also  can  increase  the  effective  hole  mobility, 
which  is  proportional  to  exp (-EJkT)  [38],  where  Ea  is  the  hole  trap  activation  energy,  k  is  the 


Fig.  55.  Electron  mobility  vs.  gate  over-drive  voltage  (VGS- 
VT)  for  the  cases  of  no  stress,  tensile  stress  of  200  MPa,  and 
5  Mrad(Si02)  under  tensile  stress  of  200  MPa. 


Fig.  56.  Electron  mobility  enhancement  vs.  mechanical 
stress  before  and  after  5  Mrad(Si02)  irradiation. 
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Boltzman  constant,  and  T  is  temperature.  Hence,  strain-induced  lowering  of  hole  trap  energy  levels 
may  reduce  charge  trapping  in  Hf02  and  SiOx. 

Radiation  Induced  Mobillity  Degradation  under  Mechanical  Stress 

Electron  mobility  is  shown  as  a  function  of  gate  over-drive  voltage  (VGS-VT)  in  Fig.  55.  200  MPa 
of  tensile  stress  enhances  the  electron  mobility  at  all  gate  biases.  Electron  mobility  degradation 
for  devices  irradiated  to  5  Mrad(Si02)  under  200  MPa  of  tensile  stress  is  ~  1%  at  VGS-VT  =  0.55 
V,  compared  to  the  pre-irradiation  200  MPa  case,  but  the  electron  mobility  is  still  higher  than  the 
unstressed  case.  Fig.  56  shows  the  electron  mobility  enhancement  as  a  function  of  stress  before 
and  after  5  Mrad(Si02)  irradiation.  The  mobility  enhancement  compared  to  unstressed  devices  is 
positive  after  a  total  dose  of  5  Mrad(Si02)  at  all  tensile  stress  levels  above  70  MPa,  indicating  that 
the  benefit  of  strained  Si  is  not  lost  after  irradiation. 

B.  Laser-  induced  Current  Transients  in  Strained  Si  Diodes 

Single  event  transients  (SETs)  and  single  event  upsets  (SEUs)  are  related  to  collection  of  radiation¬ 
generated  charge  at  sensitive  circuit  nodes.  Although  strained-Si  technology  is  widely  adopted,  the 
effects  of  mechanical  stress  on  current  transients  generated  by  laser  or  ion  strikes  at  the  source/drain 
regions  have  not  been  reported.  It  is  important  to  understand  how  mechanical  stress  affects  these 
transient  pulses  since  the  transport  of  the  radiation-generated  carriers  in  the  substrate  is  affected 
by  stress.  Laser-induced  current  transients  on  a  uniaxially  stressed  Si  N+/P  junction  diode  are 
reported  in  this  paper.  An  N+/P  diode  is  a  good  representation  of  the  source/drain  junctions  that  are 
responsible  for  charge  collection  in  n-channel  MOSFETs.  P-channel  MOSFETs  are  also  important 
for  considering  SETs  and  SEUs.  However,  stress-induced  electron  mobility  enhancement  is  easier 
to  understand  than  that  of  holes,  so  N+/P  diodes  are  used  in  this  work. 

The  shapes  of  current  transients  and  the  amount  of  collected  charges  are  measured  as  a  function 
of  stress,  because  both  of  them  are  crucial  in  predicting  SETs  and  SEUs  in  circuits.  Controlled 
external  mechanical  stress  is  applied  via  a  four-point  bending  jig  while  the  samples  are  irradiated 
using  a  picosecond  pulsed  laser.  The  characterization  system  is  based  on  the  direct  measurement  of 
the  current  transients.  The  FLOODS  simulation  tool  is  used  to  explain  the  mechanisms  responsible 
for  the  differences  in  charge  collection  between  stressed  and  unstressed  devices.  Additionally,  the 
simulations  provide  insight  into  the  effects  of  high  mechanical  stress  (~1  GPa)  on  laser-induced 
current  transients,  above  the  maximum  stress  that  could  be  applied  using  the  four-point  bending  jig 
(240  MPa  on  these  samples). 

Experimental  Details 

The  laser-induced  current  transients  in  mechanically-stressed  Si  N+/P  diodes  are  captured  using  a 
high-speed  measurement  system  with  an  integrated  four-point  bending  jig,  as  shown  in  Fig.  9.  The 
samples  used  in  this  study  are  N+/P  diodes  fabricated  on  (001)  Si  wafers  using  a  standard  130-nm 
CMOS  technology.  The  active  area  of  the  diodes  is  50  jim  x  100  pim.  Nickel  silicide  (NiSi),  silicon 
oxide  (SiOx),  and  copper  (Cu)  patterns  are  present  on  top  of  the  diodes  as  shown  in  Fig.  10  using 
transmission  electron  microscopy  (TEM)  and  energy-dispersive  X-ray  spectroscopy  (EDS).  The 
thickness  of  the  NiSi,  SiOx,  and  Cu  patterns  is  ~20  nm,  720  nm,  and  280  nm,  respectively.  The 
doping  densities  of  the  n+,  p-well,  and  p-substrate  are  ~1020,  ~1018,  and  ~1016  cm'3,  respectively. 
The  shallow  trench  isolation  (STI)  in  these  devices  is  an  additional  source  of  mechanical  stress. 
However,  STI-induced  stress  in  the  center  of  the  large  diodes  is  negligible.  Uniaxial  mechanical 
stress  along  the  <110>  direction  is  applied  using  a  four-point  bending  jig. 
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Fig.  57.  Laser-induced  current  transient  measurement  system  using  a 
four-  point  bending  jig. 
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Fig.  58.  Schematic  of  N+/P  diode  structure  through  TEM 
and  EDS  analysis  (not  to  scale)  and  TEM  image. 

A  cavity-dumped  dye  laser  with  a  wavelength  of  590  nm,  a  pulse  energy  of  218  pJ,  and  a  pulse  width 
of  1  ps  is  used  to  inject  electron-hole  pairs  in  the  diode.  The  laser  direction  is  normally  incident  to 
the  diode  surface  and  has  a  spot  size  of  12// m  in  diameter.  The  peak  carrier  concentration  produced 
by  the  laser  is  ~1.6  x  10 19  cm'3.  The  pulse  laser  energy  reaching  the  diode  active  area  is  smaller 
than  the  value  measured  at  the  surface  of  the  structure  due  to  the  optical  properties  of  the  layers  on 
top  of  the  diode.  The  transient  measurement  system  uses  a  Tektronix  TDS8200  digital  sampling 
oscilloscope  with  80E03  sampling  module  (20  GHz  bandwidth)  and  is  connected  to  the  device 
using  a  bias  tee  (10  kHz  to  40  GHz)  and  a  ground-signal-ground  (GSG)  probe  tip  (DC  to  40  GHz). 

Current  transients  on  the  N+/P  diode  are  measured  under  different  values  of  stress  (160  MPa  and 
240  MPa  tensile,  no  stress,  and  160  MPa  compressive)  with  a  5  V  reverse  bias.  Although  the 
maximum  applied  stress  (~240  MPa)  in  this  investigation  is  about  16%  of  that  produced  by  process- 
induced  stressors  (~  1.5  GPa),  the  experiments  still  show  the  dominant  mechanisms  in  the  effects 
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Fig.  59.  Laser-induced  current  transients  and  the  ratio 
of  collected  charge  measured  as  a  function  of  <1 10> 
uniaxial  mechanical  stress. 


of  stress  on  SETs;  this  approach  is  analogous  to  previous  works  describing  the  effects  of  stress  on 
unirradiated  MOS  devices. 

Experimental  Results 

The  effect  of  the  applied  mechanical  stress  on  maximum  current  and  charge  collection  is  characterized 
by  monitoring  laser- induced  current  transients  at  each  uniaxial  stress  value.  Increasing  tensile  stress 
results  in  lower  maximum  currents  (/  )  and  collected  charges  ( Q )  than  those  measured  under  no 
stress,  as  shown  in  Fig.  59  for  times  up  to  10  ns  after  the  laser  pulse  strikes  the  device.  Each 
transient  curve  is  measured  using  an  averaging  technique  (100  points)  in  the  sampling  oscilloscope. 
Q  is  obtained  by  integrating  the  measured  transient  as  a  function  of  time.  The  data  points  are 
the  average  Q  at  each  level.  The  error  bars  in  the  data  points  represent  the  standard  deviation  in 
the  data  at  that  stress  level.  Opposite  to  tensile  stress,  compressive  stress  increases  Imax  and  Q.  A 
decrease/increase  in  Imax  and  Q  under  tensile/compressive  stress  can  be  explained  by  a  1-D  transient 
analytical  solution. 

Current  transients  (. I(t ))  are  proportional  to  N]inU in  the  solution,  where  A  is  the  number  of  laser¬ 
generated  electron-hole  pairs  and  pinU  is  electron  mobility  along  the  <00 1>  direction.  ]inU  is  the 
dominant  contribution  for  electron  mobility,  because  electrons  are  mainly  moving  in  the  <00 1> 
direction  due  to  applied  field  along  the  <00 1>  direction  for  the  large  diodes  used  in  the  experiment. 
A  as  a  function  of  depth  in  Si  is  defined  as 

N(z)  -  ?j&exp(-<x(o)z)f  xI0(z,t)dt  (1) 

where  a  is  the  absorption  coefficient  of  Si,  hco  is  the  photon  energy  (2.1  eV),  z  is  depth  in  the  Si,  IQ 
is  the  intensity  of  the  laser  beam,  and  t  is  the  time,  a  depends  on  the  band  gap,  where  a  normalized 
stress  dependent  a  is  defined  as 


A  a(o) 
a 


A£  (a) 
ha)-Et 


«\, 


hco  >  Ex 


(2) 


where  a  is  the  mechanical  stress,  and  E  is  the  Si  band  gap  (1.12  eV).  Based  on  (1)  and  (2),  a  change 
in  A  as  a  function  of  uniaxial  tensile  stress  and  the  depth  into  Si  is  plotted  in  Fig.  12.  Since  stress- 
induced  bandgap  narrowing  in  Si  is  minimal  (~0.01  eV  at  240  MPa  of  both  compressive  and  tensile 
stress),  an  increase  in  a  of  Si  for  this  range  of  mechanical  stress  is  negligible.  As  a  result,  there  is 
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Fig.  60.  The  number  of  laser-generated  electron-hole 
pairs  as  a  function  of  depth  (2)  and  <11 0>  uniaxial 

no  significant  increase  calculated  in  N  at  each  depth  under  mechanical  stress,  less  than  1%  at  240 
MPa  of  tensile  stress,  as  shown  in  Fig.  60. 

Since  the  change  in  N  due  to  stress  is  minimal,  a  ~6.5%  decrease  in  I  at  240  MPa  of  tensile 
stress  is  caused  mainly  by  a  decrease  in  electron  mobility  along  the  <001>  direction.  Likewise, 
for  compressive  stress,  an  increase  in  /  ^results  from  an  increase  in  electron  mobility  along  the 
<00 1>  direction,  because  strain-induced  band  gap  narrowing  is  very  small  (~0.01  eV  at  240  MPa). 
This  suggests  that  a  decrease/increase  in  electron  mobility  along  the  <00 1>  direction  under  tensile/ 
compressive  stress  {Apt  ^  results  in  a  decrease/increase  in  Imay.  The  experimental  results  and 
qualitative  analysis  both  can  be  explained  by  previous  results  on  piezoresistance  (jt)  coefficients 
in  Si. 


n=  /i(a)-/i(Q)  A  n 

0)  '  afi(0) 

The  Jt  coefficient  represents  changes  of  mobility  resulting  from  applied  stress, 

where  pi  (a)  and  pi(0)  are  the  mobility  with  and  without  stress,  respectively,  and  Apt  is  the  change 
in  the  mobility.  Changes  (increase  or  decrease)  of  the  electron  mobility  result  from  changes  of  the 
average  electron  effective  mass  (m*),  due  to  repopulation  of  electrons  under  mechanical  stress. 
For  example,  Fig.  61  shows  that  tensile  stress  splits  the  conduction  bands  into  A^  and  A4.  Electrons 
repopulate  from  the  A4  valley  into  the  valley.  Average  effective  mass  along  the  <110>  direction 
(mz/)  decreases  under  tensile  stress  in  the  same  direction,  but  average  effective  mass  along  the 
<001>  direction  (m^  increases  under  tensile  stress  in  the  <110>  direction.  Thus,  <110>  tensile 
stress  decreases  the  electron  mobility  along  the  <00 1>  direction  (//nD),  because  pi  is  inversely 
proportional  to  m*.  The  opposite  dependence  is  expected  with  compressive  stress.  The  concept  of 
the  Jt  coefficient  is  implemented  in  current-transient  simulations  for  diodes  under  mechanical  stress 
in  the  next  section.  Q  is  also  proportional  to  the  funneling  length,  L  =  (1  +pi nJ]A where pipU is  the 
hole  mobility  along  the  <00 1>  direction,  and  W  is  the  depletion  width  .  A  change  in  hole  mobility 
along  the  <00 1>  direction  (Apt ^  under  uniaxial  mechanical  stress  is  negligible  (~  0.3%  at  250 
MPa).  Therefore,  the  change  of  the  collected  charge  ( AQ )  as  a  function  of  the  applied  mechanical 
stress  is  also  dominated  by  the  change  in  electron  mobility  (Api^.  By  applying  the  same  concepts 
for  analyzing  I(t)  and  L  as  in  the  N+/P  diodes  above,  it  is  possible  to  predict  how  current  transients 
in  P+/N  diodes  would  change  under  mechanical  stress.  Since  I(t)  is  proportional  to  NpipU  and  ApipU 
(~  1%  at  1  GPa)  is  not  significant,  the  change  in  Imax  is  not  expected  to  be  significant  under  stress. 
Q  is  likely  to  increase  with  tensile  and  decrease  with  compressive  stress,  because  L  is  equal  to 
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Fig.  61.  Uniaxial  tensile  stress  effect  on  electron  mobility. 


( 1  +pi ^/pi^W  and  affected  by  Zl///0.  However,  further  experimental  data  are  required  to  verify  the 
analytical  analysis  of  current  transients  in  P+/N  diodes  under  mechanical  stress. 

TABLE  2. 

Values  of  piezoresistance  (n)  coefficients  (10 5  MPa  ')  used  in  FLOODS[57] 


Kn 

Kl2 

K44 

Electron 

-102.2 

53.4 

-13.6 

Hole 

6.6 

-1.1 

138.1 

TCAD  Simulation  and  Results 


Based  on  the  previous  experimental  analysis,  FLOODS  simulations  are  performed  to  understand 
the  mechanisms  of  carrier  transport  under  uniaxial  stress  and  to  predict  how  high  stress  (~lGPa) 
affects  the  current  transients  in  diodes.  The  Masetti  and  Brooks-Herring  mobility  models  are  used 
to  account  for  carrier  transport  in  a  high  injection  case.  Shockley-Read-Hall  and  Auger  band-to- 
band  recombination  models  are  also  considered.  The  number  and  distribution  of  electron-hole  pairs 
generated  by  the  laser  pulse  is  calculated  by  a  single-photon  absorption  (SPA)  equation. 


Before  analyzing  the  effects  of  stress  on  current  transients,  baseline  simulations  under  no  stress  are 
performed.  These  results  are  matched  to  the  measured  current  transient  under  no  stress.  It  is  very 
important  to  understand  the  physics  that  dominates  current  transients  in  an  unstressed  case  in  order 
to  predict  the  results  under  a  stressed  case.  A  2-dimensional  simulation  structure,  shown  in  Fig.  62, 
is  built  based  on  analysis  of  the  structure  and  material  of  the  N+/P  diodes,  as  shown  in  Fig.  58.  The 
width  and  depth  of  the  diodes  are  100  //m  and  10  //m,  respectively.  The  SiOx,  Cu  dummy  patterns, 
and  NiSi  are  not  implemented  in  the  simplified  FLOODS  simulations.  However,  the  omitted  layers 
can  reduce  laser  energy  due  to  the  reflection,  absorption,  and  transmission  properties  of  each 
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Fig.  62.  Schematic  of  laser-induced  current  transients  and 
2-dimensional  simulation  structure  of  an  n+p  diode. 


Fig.  63.  Simulated  energy  dependence  of  laser-induced 
current  transients. 


material.  Fig.  63  shows  that  a  decrease  in  laser  pulse  energy  results  in  a  decrease  in  the  peak  current 
and  charge  collection.  The  simulated  result  for  the  case  of  pulse  energy  of  13.5  pJ  agrees  with  the 
experiment  result,  as  shown  in  Fig.  63.  However,  the  pulse  energy  (13.5  pJ)  is  much  different  from 
the  measured  pulse  energy  (218  pJ).  The  discrepancy  can  be  explained  by  the  absorption,  reflection, 
and  transmission  properties  of  each  layer  over  the  active  region  of  the  diode.  Since  43%  of  the 
spot  area  of  incident  laser  is  occupied  bv  Cu  dummy  patterns  which  block  the  laser,  only  57%  of 
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incident  laser  energy  is  transmitted.  Next,  only  78%  of  the  energy  is  transmitted  through  the  720 
nm  SiOx,  due  to  reflection  losses  at  the  interfaces.  Lastly,  16%  of  the  energy  is  transmitted  through 
NiSi.  Therefore,  the  calculated  pulse  laser  energy  reaching  the  diode  active  area  is  ~15.5  pJ,  ~7.1% 
(=  0.57  X  0.78  x  0.16)  of  the  incident  energy.  If  the  thickness  of  each  layer  varies  by  ~10%  and 
the  composition  of  the  NiSix  also  varies,  the  calculated  laser  energies  range  from  12.5  to  22  pJ. 
The  collected  charge  in  the  simulation  (12.8  pC)  agrees  well  with  the  average  collected  charge 
in  the  experiment  (12.3  pC).  As  a  result,  the  amount  of  energy  used  in  the  simulations  to  produce 
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Fig.  64.  Simulated  laser-induced  current  transients  as  a  function  of  <11 0>  uniaxial 
mechanical  stress. 

agreement  with  the  experiments,  13.5  pJ,  is  reasonable.  A  piezoresistive  mobility  model  based  on 
Smith’s  JT-coefficients  is  used  to  consider  mobility  enhancement  under  mechanical  stress,  as  shown 
in  Table  2.  The  6x6  piezoresistive  model  is  defined  as 

where  jt  ,  o  ,  q  ,  and  fx  are  components  of  the  piezoresistance  coefficient,  mechanical  stress, 
resistivity,  and  carrier  mobility,  respectively,  and  A q  /  q  and  A/x  /  /x  are  fractional  changes  in 
resistivity  and  mobility. 

The  doping  dependence  of  the  jr-coefficients  is  considered.  From  (4),  currents  are  calculated  as  a 
function  of  mechanical  stress.  Current  densities  are  expressed  as 


where  J.(0)  and  J.(o)  are  current  density  components  with  and  without  stress  based  on  a  Cartesian 
coordinate  system,  respectively.  The  simulated  current  transients  in  Fig.  64  show  the  same  trend 
as  the  experimental  data  in  Fig.  59.  Imax  and  Q  in  the  simulations  also  agree  with  the  experiments, 
as  shown  in  Figs.  65  and  66.  The  data  points  in  the  experiments  are  the  average  Imax  and  Q  at 
each  stress  level.  The  error  bars  in  the  data  points  represent  the  standard  deviation  in  the  data  at 
each  stress  level.  The  simulation  results  predict  that  I  and  O  under  1  GPa  of  tensile  stress  will 
decrease  by  ~23%  and  ~21%,  respectively.  Analogous  to  tensile  stress,  1  GPa  of  compressive  stress 
increases  Imax  and  Q  by  17%  and  13%,  respectively.  These  experiment  and  simulation  results  for 
strained  N+/P  diodes  show  that  uniaxial  stress  changes  the  shape  of  current  transients  and  collected 
charges. 

To  apply  the  strain  engineering  concepts  for  mitigating  SETs  and  SEUs  in  deep  submicron  devices, 
we  need  to  understand  the  differences  between  these  large  diodes  and  deep  submicron  devices. 
Due  to  the  scaling  of  devices,  the  size  of  source/drain  junctions  is  smaller  than  that  of  a  radiation¬ 
generated  e-h  pair  cloud.  While  out-of-plane  transport  of  those  carriers  dominates  current  transients 
in  these  large  diodes,  both  out-of-plane  and  in-plane  transport  of  radiation-generated  carriers  under 
mechanical  stress  should  be  considered  for  scaled  devices.  Mechanical  stress  can  either  enhance 
or  degrade  mobility,  depending  on  the  orientation.  For  example,  <110>  uniaxial  tensile  stress 
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Fig.  65.  Peak  current  (Imax)  as  a  function  of  mechanical  stress,  (positive  (+)  :  tensile, 
negative  (-):  compressive) 
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Fig.  66.  Collected  charges  (Q)  until  10  ns.  (positive  (+)  :  tensile,  negative  (-): 
compressive) 

increases  electron  mobility  along  the  <110>  direction,  but  decreases  it  in  the  <00 1>  direction.  For 
commercial  off-the-shelf  (COTs)  chips,  strained-Si  technology  is  implemented  to  increase  carrier 
mobility  along  the  channel  (in-plane)  direction.  It  is  necessary  to  evaluate  whether  strained-Si  have 
beneficial  or  detrimental  trade-offs  between  chip  performance  and  radiation  effects. 

C.  Summary  and  Conclusion 

Positive  charge  trapping  is  the  dominant  radiation-induced  degradation  mechanism  in  both 
unstressed  and  mechanically  stressed  Hf02-based  nMOSFETs.  Uniaxial  tensile  and  compressive 
stresses  in  nMOSFETs  decrease  the  amount  of  net  positive  charge  trapping  and  reduce  the  threshold 
voltage  shift.  This  is  attributed  to  enhanced  detrapping  of  holes  or  compensating  electron  trapping 
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in  Hf02  or  SiOx,  and/or  increasing  effective  hole  mobility  in  the  gate  dielectric.  Electron  mobility 
enhancement  produced  by  the  stress  is  retained  in  the  irradiated  devices.  Uniaxial  strain  engineering 
has  the  potential  to  maintain  drive  current  (mobility)  enhancement  significantly  in  advanced  Hf02- 
based  MOSFETs,  even  after  irradiation  to  the  relatively  high  total  doses  (5  Mrad(Si02))  reported 
here. 

Uniaxial  tensile  stresses  in  Si  N+/P  diodes  decrease  the  maximum  peak  currents  and  collected 
charges  for  laser-induced  current  transients.  Quantitative  analysis  and  FLOODS  simulation  results 
suggest  that  this  can  be  attributed  to  the  degradation  of  electron  mobility  along  the  <00 1>  direction. 
Unlike  uniaxial  tensile  stress,  uniaxial  compressive  stress  shows  the  opposite  trend.  Therefore, 
uniaxial  strain  engineering  has  the  potential  to  control  the  shape  of  single  event  transients  and 
the  amount  of  charges  collected  in  devices.  Furthermore,  these  results  suggest  that  strained-Si 
technology  has  a  significant  impact  on  SETs  and  SEUs  at  the  circuit  level. 
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3.23  Analysis  and  Modeling  of  Total-lonizing-Dose  (TID)  Effects 


This  report  summarizes  technical  results  and  key  findings  from  Arizona  State  University  for 
the  Vanderbilt  University  AFOSR  MURI  program  on  radiation  effects.  Over  the  past  five  years. 
ASU  efforts  have  focused  on  the  analysis  and  modeling  of  total-ionizing-dose  (TID)  effects  in 
semiconductor  materials  and  devices.  The  specific  topics  summarized  in  this  final  report  pertain  to 
ASU’s  research  on: 

•  basic  mechanisms  of  defect  buildup  in  isolation  oxides 

•  total  ionizing  dose  effects  characterization  and  modeling  in  advanced  CMOS  transistors 

The  ASU  effort  was  directed  by  Dr.  Hugh  Barnaby.  Several  graduate  students  were  supported 
during  the  program  including  three  Ph.D.  students:  Xiao-Jie  Chen,  Michael  Mclain,  and  Ivan 
Sanchez  Esqueda.  Both  Mclain  and  Chen  completed  their  doctorates  during  the  program. 

The  results  of  the  modeling  and  experimental  work  performed  by  the  ASU  team  have  lead  to 
significant  improvements  in  the  basic  models  that  describe  ionizing  radiation-induced  defect 
buildup  semiconductor  oxides.  Modeling  refinements  have  helped  identify  key  mechanisms 
related  to  dose  rate  sensitivity,  in-package  hydrogen  contamination,  response  saturation,  and 
defect  annealing.  Physical  models,  primarily  expressed  as  systems  of  differential  equations,  were 
implemented  in  numerical  simulators  and  the  results  compared  to  experimental  data  for  validation. 
The  rigorous  model  validation  and  parameterization  performed  during  the  course  of  the  program 
was  the  first  of  its  kind  for  advanced  CMOS  technologies  (below  the  100  nm  node).  In  addition  to 
the  physical  model  development,  techniques  for  device  level  analytical  modeling  of  TID  effects, 
from  defect  generation  to  radiation-enabled  compact  modeling  were  developed  at  ASU  for  the 
MURI.  The  techniques  were  integrated  into  a  predicative  technology  modeling  (PTM)  package  for 
radiation  effects  that  generates  validated  compact  models  for  radiation  damage  from  user  specified 
input  conditions  (e.g.,  environments  and  system  states).  The  general  flow  of  the  PTM  for  radiation 
methodology  is  summarized  in  Fig.  67.  This  effort  produced  for  the  first  time  surface-potential- 
based  models  for  ionizing  radiation  effects  in  both  bulk  and  SOI  CMOS. 

A.  Basic  Mechanisms  for  Defect  Buildup 


Fig.  67.  Flow  diagram  describing  approach  to  predictive  technology  modeling  for 
radiation  effects. 


Model  Overview 
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In  this  task,  models  describing  the  processes  involved  in  radiation-induced  defect  generation  were 
adapted  and  in  some  case  extended  to  address  effects  in  emerging  technologies.  These  processes 
include  those  which  lead  to  the  buildup  of  trapped  oxide  charge  (Not)  in  isolation  oxides  and  traps 
at  the  interface  between  dielectric  and  semiconductor  materials.  The  processes  for  Not  generation, 
illustrated  in  Fig.  68,  include  electron-hole  pair  generation,  geminate  recombination,  charge 
transport,  and  trapping.  The  key  differential  equations  used  to  model  these  processes  are 

and 


ft 

dx 


(i) 


Eq.  1  captures  the  relationship  between  ionizing  dose  rate  (D )  and  hole  flux  transport  (f )  and  Eq. 


St 


UJt) 


(2) 


2  models  the  buildup  and  annealing  of  trapped  oxide  charge  as  a  function  of  carrier  reactions  with 


Fig.  68.  Processes  for  oxide  trapped  charge  buildup. 


fixed  defect  sites. 

The  interface  trap  (Nit)  model ,  illustrated  in  Fig .  69 ,  shares  similar  processes  with  the  charge  trapping 
model,  e.g.;  carrier  generation,  geminate  recombination,  and  flux-based  carrier  transport.  The  Nit 
model  is  differentiated  by  its  dependence  on  hydrogen  content  in  the  oxide  (or  other  dielectric) 
during  radiation  exposure.  This  was  demonstrated  experimentally  in  this  program.  For  interface 
traps,  key  model  equations  include 

and 


Vu 

ax 


-  Ndh  oDHfp 


air 

dt 


(3) 
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Eq.  3  captures  the  relationship  between  hole  and  hydrogenated  defect  reactions  and  proton  flux 


av, 

dt 


(4) 


transport  (fR)  and  Eq.  4  models  the  buildup  and  annealing  of  interface  traps  as  a  function  of  proton 
reactions  with  passivated  bonds. 


Hydrogen  and  Dose  Rate  Effects 


•  -surviving  hole  (p) 

^  -  hole  trap  (NT) 

Q  -  trapped  hole  (Not) 

f  -  hole  flux 
p 

area  =  0(0 


Fig.  69.  Processes  for  interface  trap  buildup. 


Experiments,  conducted  by  ASU,  revealed  that  exposure  of  a  device  oxide  to  molecular  hydrogen 
prior  to  irradiation  can  have  a  very  significant  impact  on  the  build-up  of  interface  traps.  This  effect 
is  shown  in  Fig.  70.  In  this  plot  the  density  of  interface  traps  after  a  fixed  radiation  dose  (specifically 
30  krad(Si))  is  shown  to  increase  as  the  H2  content  in  the  gas  surrounding  the  sample  is  increased. 
These  results  were  reported  by  Jie  Chen  after  he  performed  combined  irradiations  and  hydrogen 
soaks  on  annular  transistors  fabricated  in  a  commercial  bipolar  process  [1].  Chen  also  developed 
a  model  for  the  observed  effect  based  on  the  diffusion  of  H2  from  the  gas  ambient  into  the  device 
oxide.  Reactions  between  radiation-induced  carriers  in  the  oxide  and  hydrogen  species  were  used 
to  derive  the  closed  form  expression 


which  relates  molecular  hydrogen  concentration  to  the  post-irradiation  interface  trap  density. 

f  I  \ 


AM, 


W,)2 

I  +  ki(N„fj 


(5) 


Parameters  for  the  expressions  ([],  kp  and  k2)  were  shown  to  be  related  to  defect  density  in  the 
oxide,  reaction  cross-sections  for  holes,  protons,  and  interfacial  defect  sites,  oxide  thickness,  total 
dose,  hydrogen  diffusivity,  and  carrier  yield.  Fig.  70  shows  how  the  model,  with  appropriately 
chosen  parameters,  can  give  an  excellent  fit  to  the  experimental  data. 

Chen  and  Barnaby  followed  up  on  this  work  by  adding  another  important  variable,  irradiation 
dose  rate,  into  the  input  matrix.  The  results  and  models  describing  the  impact  of  hydrogen  on 
dose  rate  response  in  oxides  were  first  reported  in  this  MURI  program.  Fig  71  plots  N.t  buildup 
vs.  dose  rate  for  devices  exposed  to  different  levels  of  H2  (100%,  1%  and  <  0.01%)  prior  to  and 
during  radiation  exposure.  For  these  experiments  the  total  dose  level  was  fixed  at  30  krad  (Si).  As 
the  figure  shows,  increased  levels  of  hydrogen  raise  the  high  and  low  dose  rate  saturation  levels 
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and  shift  the  transition  dose  rate  to  higher  levels.  Chen  reproduced  these  effects  by  modeling  hole 
reactions  with  hydrogenated  defects  (adapted  from  Mclean  and  others)  and  annihilating  reactions 
between  electrons  and  holes.  Eqs.  6  and  7  are  the  key  reactions  in  Chen’s  adapted  model: 

In  (6),  DH  is  a  hydrogenated  defect  in  the  oxide.  In  (7)  Chen  adds  electron  and  hole  recombination. 


DH  +  p  ^  D  +  H+, 

CO 

p  +  X. 

(7) 

Fig.  70.  Fit  of  model  described  by  (5)  to  data  obtained  from  experiments. 


Fig.  71.  Measured  N.t  as  a  function  of  dose  rate  for  different 
hydrogen  concentrations  in  thermal  oxides 


Fig.  72  shows  the  simulated  model  results. 
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Fig.  72.  Simulation  of  N.  formation  as  a  function  of  dose  rate  for 
different  hydrogen  concentrations. 
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Fig.  73.  Layout  and  cross-sectional  view  of  the  NW 
FOXFET. 

Models  for  TID  Defect  Buildup  in  CMOS  Isolation  Oxides 

In  this  effort,  CMOS  field  oxide  test  structures,  FOXFETS  and  FOXCAPs,  were  designed  and 
fabricated  in  order  to  provide  data  for  model  validation.  While  several  design  variants  were 
developed  for  the  program,  the  most  effective  characterization  tool  was  an  n-well  (NW)  FOXFET 
devices  fabricated  in  a  90  nm  commercial  bulk  CMOS  low- standby  power  (LSP)  technology.  Fig. 
73  provides  illustrations  of  the  layout  and  cross-sectional  view  of  the  standard  NW  FOXFET  used 
for  the  testing.  CMOS  test  structures  were  irradiated  in  several  60Co  irradiation  sources  (ASU,  BAE 
systems,  and  NRL)  so  that  the  responses  to  different  dose  rates  could  be  determine.  Fig.  74  plots 
oxide  trapped  charge  extracted  from  the  FOXFET  data  as  a  function  of  irradiation  and  anneal  time. 
These  data  indicate  similar  trends  in  the  total  dose  responses  and  anneal  rates  (i.e.  no  strong  dose 
rate  dependence  observed). 

Data  sets,  such  as  the  one  in  Fig.  74,  was  used  to  build  and  verify  an  analytical  model  for  the 
various  processes  involved  in  trapped  charge  formation  and  annealing  in  isolation  oxides.  The  key 
equation,  expressed  formally  in  [7],  is 


Anot=D Atgo ntapfy.pxp  -  notcmfyjnxn-m- Atln(  1 0  ]  -t  * 


(8) 
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Fig.  74.  Not  vs  time  extracted  during  and  after  60Co  exposures 
at  different  dose  rates. 


which  describes  the  dependence  of  trapped  charge  buildup  and  removal  in  oxides  on  the  flux  and 
trapping  of  carriers  in  the  oxide  (i.e.,  both  holes  and  electrons)  and  a  time  dependent  annealing 
function.  As  illustrated  in  Fig.  75,  the  analytical  model  (solid  line  labeled  ‘ANYL’)  accurately 
reproduces  the  amount  of  trapped  charge  buildup  as  a  function  of  radiation  exposure  time  for  dose 
rates  of  1800  rad/s,  344  rad/s,  and  20  rad/s.  The  experimental  data  in  the  figures  are  represented  by 
the  symbols  and  labeled  ‘EXP’. 


Fig.  75.  Not  vs.  time  data  and  analytical  model  fits  for  three  different  dose  rates. 


B.  TID  Effects  Characterization  and  Modeling  in  Advanced  CMOS  Transistors 

Overview 

In  the  MURI  program,  ASU  utilized  the  PSP  compact  modeling  framework  to  build  analytical 
models  for  TID  effects  in  CMOS  transistors.  Radiation-enabled  compact  models  were  constructed 
with  highly  accurate  surface-potential  based  formulations.  Model  development  was  performed 
on  both  bulk  and  SOI  CMOS  technologies.  Extensive  experiments  and  TCAD  simulations  were 
conducted  to  validate  the  models. 

Bulk  CMOS 

Ionizing-radiation  effects  in  bulk  CMOS  devices  were  incorporated  into  the  industry  standard 
surface-potential  based  compact  model  (PSP).  The  inclusion  of  total  ionizing  dose  (TID)  effects 
into  PSP  was  accomplished  by  renormalization  of  the  surface  potential  equation  (SPE)  in  order 
to  make  the  PSP  formulation  applicable  to  irradiated  devices.  Model  verification  was  performed 
via  comparison  to  experimental  data  obtained  from  the  NW  FOXFET  devices  fabricated  in  a  90 
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nm  low-standby  power  commercial  bulk  CMOS  technology.  The  model  accurately  describes 


contributions  to  off-state  currents  from  inter-device  leakage  and  demonstrates  its  capabilities  for 
simulating  radiation-induced  degradation  in  advanced  CMOS  integrated  circuits  (ICs).  Eqs.  9-11 
represent  the  key  functions  for  the  model. 

(9) 

V FB  “  ~  0^  +  Ar $b  )  aU^ 

Vi 

(10) 

4-y  /vA  fata  + 

(11) 

Eq.  9  is  the  implicit  surface  potential  equation  (SPE).  Eq.  10  is  the  flatband  equation,  which 
captures  the  impact  of  the  radiation-induced  defects.  Eq.  11  is  the  drain  current  equation  using  the 
symmetric  linearization  method.  ASU  reported  the  results  in  [8].  Through  renormalization  of  the 
SPE,  ASU  researchers  were  able  to  build  compact  models  that  produced  extremely  accurate  fits 
to  the  transistor  (FOXFET)  data,  both  prior  to  and  after  high  levels  of  radiation  exposure.  As  an 
example,  Fig.  76  shows  the  alignment  between  the  Id  vs.  Vgs  data  (symbols)  and  model  (solid  line). 


Fig.  76.  Id-V  s  characteristics  obtained  analytically  (solid  lines)  and  experimental 
data  (symbols)  at  the  different  levels  of  TID  (Vd  =  0.1  V,  Vg  =  Vfe  =  0  V). 


Silicon-on-Insulator  CMOS 

ASU’s  research  effort  on  SOI  devices  for  the  MURI  program  focused  on  modeling  radiation  effects 
on  single  and  multiple  gate  devices.  For  the  conventional  single  gate  transistor  (shown  in  Fig.  77), 
ASU  developed  a  model  for  the  GIDL  enhancement  of  back  channel  leakage  which  was  reported 
previously. 

In  the  model,  holes  injected  into  the  floating  body  from  a  band-to-band  tunneling  current  (JpBBT) 
charge  the  body  of  an  n-channel  SOI  FET  to  a  potential  VB.  This  effect  is  expressed  mathematically 
as, 
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(12) 


Radiation  damage  to  the  thick  buried  oxide  introduces  a  defect  potential  which  is  proportional  to  the 
sum  of  oxide  trapped  and  interface  trap  defects.  This  defect  potential  is  expressed  mathematically 
as, 


K-Z?-b*-D.h  -!*/)] 


(13) 


c 


The  combined  effect  of  damage  to  the  buried  oxide  and  body  charging  from  GIDL  significantly 
reduced  the  back-gate  threshold  voltage  and  enhanced  leakage  current  in  floating  body  SOI  devices. 
ASU  successfully  modeled  this  effect  with  calibrated  Eqs.  12,13,  and 


Eq.  14  is  the  modified  double  gate  SPE  used  for  SOI  transistors. 

In  addition  to  the  conventional  single  gate  SOI  modeling  effort,  ASU  also  performed  an  extensive 
examination  of  radiation  effects  on  multiple  gate  transistors.  This  study  revealed  a  strong  geometric 
dependence  of  TID  sensitivity  in  SOI  multiple  gate  FETs,  specifically  FinFETs  structures.  This 
dependence  was  related  to  the  FinFET  width.  Shown  in  Fig.  79  are  2D  TCAD  contour  plots  of  two 
FinFETs  used  in  the  analysi.  The  devices  are  identical  except  in  the  width  of  the  silicon  film  (tsi) 
which  defines  the  fin.  For  the  structures  shown  in  Fig.  79,  one  fin  width  is  10  nm  and  the  other  is  50 
nm.  The  figure  illustrates  the  potential  contour  in  the  silicon  body  as  a  result  of  charge  buildup  in 
the  buried  oxide.  As  the  numerical  results  show,  the  wider  devices  is  more  responsive  to  the  defects, 
thus  more  sensitive  to  radiation  damage.  Sanchez  Esqueda  et  al.  derived  an  analytical  model  for 
this  effect  which  they  reported.  Fig.  80  shows  how  the  model  compares  to  the  results  of  radiation 
experiments  on  a  sub  100  nm  FinFET  technology.  In  this  figure,  model  predictions  of  increased 
off-state  leakage  current  vs.  dose  is  compared  to  experimental  data  on  devices  with  80  nm  and  40 
nm  fin  widths.  The  data  was  obtained  in  collaboration  with  researchers  at  Vanderbilt  Univerisity. 
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Fig.  77.  Cross-section  illustrating  three  current  processes  related  to  the  radiation 
response  on  I-V  characteristics  of  FD-SOI  transistor 
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Fig.  78.  Plot  of  ID  vs.  VGF  characteristics  obtained 
experimentally  (symbols)  as  well  as  the  simulated  curves  (solid 
lines)  with  VD  =  1.4  V. 


Fig.  79.  Numerical  simulation  results  showing  greater  impact  of  Not 
buildup  in  wider  FinFET  designs. 
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Fig.  80.  Theshold  voltage  shift  vs.  dose  for  different  fin  widths  (solid 
line  is  the  model  and  symbols  are  experimental  data) 

3.24  Theory  of  Phenomena  Induced  by  Energetic  Ion  Beams 

A.  Stopping  Power  of  Channeled  Ions  in  Si 

When  an  ion  beam  hits  a  crystal,  the  primary  event  is  energy  transfer  from  the  ions  to  the  solid. 
The  energy  transfer  may  induce  atomic  displacements.  In  the  simplest  case,  an  ion  beam  can 
be  “channeled”  through  the  low-electron-density  interstitial  regions  of  Si.  The  energy  loss  to 
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Z1 


Fig.  81.  Stopping  power  of  ions  channeled  through  the  low-electron-density  channles  of 
a  Si  crystal.  Z1  is  the  atomic  number  of  the  channeled  ion.  Experimental  data  from  F.  H. 

Eisen,  Can.  J.  Phys.  46,  561  (1968). 

the  solid  is  usually  cast  as  “stopping  power”  and  is  a  function  of  the  atomic  number  Z  of  the 
ions.  Experiments  are  available  that  show  oscillations  of  stopping  power  as  a  function  of  Z.  We 
used  recently  implemented  codes  of  time-dependent  density  functional  theory,  a  state-of-the-art 
approach  for  calculating  the  dynamical  evolution  of  a  many-atom  system,  to  compute  the  stopping 
power  for  a  series  of  ions  (boron  to  calcium  in  the  periodic  table)  channeling  in  <110>  Si  channels 
The  theory  has  no  adjustable  parameters  and  produced  results  in  excellent  agreement  with  the  data 
of  Eisen  (1968)  (Fig.  81).  All  prior  calculations  included  adjustable  parameters  and  more  limited 
agreement  with  the  data.  The  calculations  demonstrated  the  need  to  include  the  full  dynamics  of 
the  electron  system  in  Si  to  capture  the  observed  effects.  [R.  Hatcher,  M.  Beck,  A.  Tackett,  and  S. 
T.  Pantelides,  “Dynamical  effects  in  the  interaction  of  ion  beams  with  solids”,  Phys.  Rev.  Lett.  100, 
103201  (2008)]. 

B.  Displacement  Damage  Differences  in  n-type  and  p-type  Si 

Radiation-induced  displacement  damage  rates  in  semiconductor  materials  generally  correlate  well 
with  the  calculated  non-ionizing  energy  loss  (NIEL)  portion  of  the  total  radiative  energy  deposition 
rate.  Developing  and  refining  models  for  the  calculation  of  NIEL  has  been  a  continuing  focus  of 
the  radiation  effects  research  community.  While  NIEL  has  enjoyed  wide  success  as  a  first-order 
predictor  of  effective  displacement  damage  rate,  it  is  well  known  that  p-  and  n-type  Si  can  exhibit 
differing  damage  rates,  as  determined  from  observed  reductions  in  effective  minority  carrier 
lifetimes  or  diffusion  lengths.  In  these  cases,  p-type  Si  exhibits  a  lower  rate  of  reduction  in  carrier 
lifetimes  than  n-type  Si  (see  Fig.  82).  This  difference  is  typically  interpreted  to  indicate  a  lower 
concentration  of  radiation-induced  defects-that  is,  less  damage-in  the  p-type  Si  samples,  due  to  an 
atomistic  mechanism  that  must  differ  in  effect  between  p-  and  n-type  Si. 

Using  state-of-the-art  quantum  mechanical  calculations,  we  have  shown  that  the  damage  rate 
differences  in  n-  and  p-type  Si  correlate  directly  to  the  fraction  of  radiation-induced  displacement 
damage  due  to  low-energy  primary  knock-on  atom  (PKA)  recoils.  These  low-energy  damage 
events  produce  isolated  vacancy-interstitial  pairs  (Frenkel  Pairs,  FPs)  that  quickly  evolve  towards 
thermodynamic  equilibrium  following  the  initial  radiation  event.  We  show  that  the  details  of  this 
evolution,  initially  controlled  by  the  fundamental  properties  of  Frenkel  Pairs  themselves,  contribute 
to  the  observed  nversus  p-type  damage  rate  differences  through  previously  reported  dopant- 
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Fig.  82.  Measured  minority  carrier  diffusion  length  damage  coefficients  for  1  Q-cm 
n-  and  p-type  Si  after  irradiation  with  various  incident  particles.  Figure  from  G.  P. 

Summers  et  al.  “Damage  correlations  in  semiconductors  exposed  to  y-radiation, 
electron-radiation  and  proton-radiation,”  IEEE  Trans.  Nucl.  Sci,  vol.  40,  no.  6,  pp. 

1372-1379,  Dec.  1993. 

dependent  stabilities  of  vacancy-impurity  complexes.  Further,  we  show  that,  for  certain  conditions, 
FPs  can  directly  contribute  to  n-  versus  p-type  damage  rate  differences  by  functioning  as  metastable, 
shallow  electron  traps  that  can  retard  carrier  transport  and  thereby  increase  the  inferred  carrier 
lifetimes.  [M.  J.  Beck,  L.  Tsetseris,  M.  Caussanel,  R.  D.  Schrimpf,  D.  M.  Fleetwood,  and  S.  T. 
Pantelides,  “Atomic-scale  mechanisms  for  Low-NIEL  dopant-type  dependent  damage  in  Si”,  IEEE 
Trans.  Nucl.  Sci.  vol.  53,  No.  6  p.  3621-3628  (2006)] 

To  study  the  details  of  the  atomic  mechanisms  that  ultimately  lead  to  differing  damage  rates  in  n- 
and  p-type  Si,  we  performed  new  calculations  of  the  fundamental  properties  of  vacancy-interstitial 
(Frenkel)  pairs  (FPs)  in  Si  using  density  functional  theory.  We  showed  that  isolated  FPs  either 
recombine,  remain  as  interacting  vacancy-interstitial  pairs,  or  dissociate  into  free  vacancies  and 
interstitials,  with  probabilities  that  depend  on  the  position  of  the  Fermi  level  in  the  Si,  i.e.,  the 
doping  type.  In  particular,  we  find  that,  in  p-Si  small  FPs  (vacancy-interstitial  separation  distances 
~8  A)  are  more  likely  to  be  stable  with  respect  to  defect  recombination  than  equivalent  FPs  in  n-Si. 
The  stability  of  FPs  against  dissociation  is  found  to  be  similar  to  that  of  vacancy-B  pairs,  with  a 
FP  binding  energy  (relative  to  the  energy  of  non-interacting,  isolated  vacancies  and  interstitials) 
of  -  0.5  eV.  Thus,  the  formation  dynamics  of  vacancy-impurity  complexes,  shown  here  to  be 
indirectly  influenced  by  the  fundamental  properties  of  FPs,  are  likely  dominant  in  determining 
observed  n-  and  p-type  Si  damage  rate  differences  for  finite  temperatures.  In  addition,  FPs  may 
directly  influence  n-  and  p-type  differences  at  low  temperatures.  [M.  J.  Beck,  L.  Tsetseris,  and  S. 
T.  Pantelides,  “Stability  and  dynamics  of  Frenkel  pairs  in  Si”,  Phys.  Rev.  Lett.  99,  215503  (2007)]. 

C.  Displacement  Damage  Formation  in  Si 

Atomic-scale  processes  during  displacement  damage  formation  have  been  previously  studied 
using  molecular  dynamics  (MD)  calculations  and  empirical  potentials.  Low-energy  displacements 
(  1  keV)  are  characterized  by  a  high  cross-section  for  producing  secondary  knock-on  atoms  and 
damage  clusters,  and  determine  the  threshold  displacement  energy  (an  important  parameter  in 
NIEL  calculations).  We  used  first-principles,  parameter- free  quantum  mechanical  calculations  of 
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the  dynamics  of  low-energy  displacement  damage  events  in  Si.  We  found  that  isolated  defects 
formed  by  direct  displacements  can  result  from  damage  events  of  100  eV.  For  higher  energy  events, 
the  initial  defect  profile,  which  subsequently  undergoes  thermal  annealing  to  give  rise  to  a  final 
stable  defect  profile,  is  the  result  of  the  relaxation  and  recry stallizatioof  an  appreciable  volume 
of  significantly  disordered  and  locally  heated  crystal  surrounding  the  primary  knock-on  atom 
displacement  trajectory.  [M.  J.  Beck,  R.  Hatcher,  R.  D.  Schrimpf,  D.  M.  Fleetwood,  and  S.  T. 
Pantelides,  “Quantum  mechanical  description  of  displacement  damage  formation”,  IEEE  Trans. 
Nucl.  Sci.  vol.  54,  No.  6  pp.  1906-1912  (2007)]. 

Displacement  Damage  Effects  In  Single-Event  Dielectric  Rupture 

Building  on  our  results  in  c-Si,  we  examined  the  effects  of  Si  and  O  recoils  in  a-Si02  systems.  In 
particular  we  probed  the  role  of  low-energy  displacement  damage  in  triggering  single-event  gate 
rupture  (SEGR),  or,  more  generally,  single-event  dielectric  rupture  (SEDR).  Using  SRIM  2008, 
we  showed  that  irradiation  of  thin  a-Si02  gate  oxides  with  heavy  ions  at  energies  associated  with 
SEGR  or  SEDR  produces  low-energy  (-100  eV)  atomic  recoils  within  the  oxide  layer.  While  the 
overall  probability  of  individual  recoils  varies  strongly  with  incident  ion  species  and  energy,  the 
distribution  of  recoil  energies  is  independent  of  both  energy  and  ion  species  for  the  high-LET 
irradiations  we  considered.  Low-energy  recoils  (100  eV)  in  a-Si02  produce  defected  regions  in 
the  amorphous  network  that  result  in  a  high  density  of  defect  states  throughout  the  a-Si02  band 
gap  (see  Fig.  83).  These  electronic  states  are  localized  on  the  recoil-induced  defects  and  are 
therefore  spatially  correlated.  The  combination  of  neighboring  defects  and  associated  defect  states 
lying  within  the  a-Si02  band  gap  results  in  a  “stepping  stone”  path  (similar  to  conducting  paths 
invoked  in  percolation  models  of  stress-  or  radiation-induced  leakage  current)  for  low  resistivity 
charge  transfer  along  the  length  of  the  recoil  damaged  region.  These  low-resistivity  paths  serve  as 
conducting  pipes  for  capacitive  discharge. 


rime  after  recoil,  fs 


Fig.  83.  Top:  Three  snap  shots  of  the  evolution  of  an  Si02  network  following  the  recoil  of  one  Si  atom.  The  colored  atoms 
indicate  the  presence  of  network  defects.  Bottom:  number  of  defects  and  their  energy  level  in  the  band  gap  as  a  function  of 

time. 
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Conducting  pipes  generated  directly  by  SHI-induced  displacement  damage  provide  a  unifying 
physical  mechanism  for  the  direct  formation  of  electrically-active  structural  defects  in  dielectric 
layers.  The  direct  formation  of  defects  as  displacement  damage  obviates  the  need  to  invoke  complex 
models  for  conversion  of  electron  excitation  energy  to  ion  kinetic  energy  to  explain  the  initial 
triggering  mechanism  for  permanent  heavy-ion  damage  to  dielectrics.  Hence  SEDR,  radiation- 
induced  soft  breakdown  (RISB),  and  long-term  reliability  degradation  (LTRD),  latent  track 
formation  and  transient  ion-induced  leakage  currents  are  damage  events  sharing,  at  least  in  part,  the 
same  physical  origin:  ion-induced  displacement  damage  leading  to  the  formation  of  low-resistivity 
conducting  pipes  (see  Fig.  84).  For  sufficiently  high  fields  applied  during  irradiation,  high  currents 
through  these  paths  will  lead  to  thermal  runaway,  and  SEDR  will  be  observed.  At  lower  fields, 
lower  currents  will  not  necessarily  lead  to  thermal  runaway,  and  RISB  will  be  observed.  At  longer 
times,  device  degradation  enhanced  by  previous  ion-induced  damage  (in  the  form  of  residual  defect 
pockets)  will  give  rise  to  LTRD.fM.  J.  Beck,  B.  R.  Tuttle,  R.  D.  Schrimpf,  D.  M.  Fleetwood,  and  S. 
T.  Pantelides,  ’’Atomic-displacement  effects  in  single-event  gate  rupture”,  IEEE  Trans.  Nucl.  Sci. 
vol.  55,  No.  6,  pp.  3025-3021  (2008)] 


Fig.  84.  Band  diagram  of  a  dielectric  stack  immediately  following 
a  heavy-ion  strike.  Filled  and  empty  circles  represent  electrons  and 
holes,  while  black  and  gray  circles  represent  thermal  and  ion-LET 
induced  carriers.  Here  we  schematically  show  the  displacement- 
damage-induced  defect  states  representing  a  physical  conducting 
path.  Bold  arrows  indicate  the  field-injected  current  flow  through 
the  oxide.  This  current  causes  thermal  runaway,  resulting  in 
permanent  structural  damage  in  the  oxide  layer. 


We  subsequently  laid  the  foundations  of  a  multi-scale  first-principles  approach  to  calculate  leakage 
currents  through  thin  dielectric  films  containing  recoil-induced  defects.  We  used  the  “source-and- 
sink”  for  quantum  transport  calculations  developed  recently  by  Varga  and  Pantelides  to  calculate 
current- voltage  characteristics  of  thin  Si02  films  between  A1  electrodes.  To  isolate  the  detailed 
effects  of  individual  atomic-scale  defects,  we  modeled  ion-induced  defects  by  removing  different 
combinations  of  oxygen  atoms  from  the  Si02  cell  and  calculating  the  transport  properties  as  a 
function  of  removed  atoms.  We  then  constructed  a  transport  model  using  three-dimensional 
percolation  theory  and  Mott  defect- to-defect  tunneling.  This  model  is  similar  to  previous  percolation 
models  of  dielectric  leakage,  and  analogous  to  the  tunneling  models  invoked  in  extensive  previous 
studies  of  trap-assisted  tunneling  (TAT).  A  key  difference  here  is  that  our  defect  distributions  are 
obtained  from  quantum  mechanical  calculations  and  evolve  with  time.  We  constructed  a  specific 
model  using  defected  Si02  thin  films  produced  by  dynamical  DFT  calculations  (Figs.  85-86).  The 
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Fig.  85  (a)  Localization  of  electronic  density  projected 
onto  a  defect  state  near  conduction  band-edge  and  (b) 
corresponding  high-electron  occupancies  in  percolation 
calculation  of  leakage  current. 


Fig.  86.  Example  of  a  percolation  path. 

model  has  two  parameters:  a  defect  radius  rQ  that  defines  the  defect’s  localization  and  controls  the 
exponential  hop,  and  a  hop  frequency  factor  vQ.  Ultimately  these  parameters  can  be  obtained  from 
quantum  mechanical  calculations.  In  the  first  implementation  of  the  model  we  assumed  a  value  rQ 
=  2  A  for  the  average  “defect  radius”  and  vQ  =1012  Hz,  as  estimated  using  experimental  data  of  hop 
attempt  rates.  With  these  values,  the  average  current  through  a  defected  Si02  structure  is  shown 
in  Fig.  87.  It  is  in  both  qualitative  and  quantitative  agreement  with  experimental  current- voltage 
characteristics  showing  RISB,  giving  us  an  initial  validation  of  the  approach.  The  results  illustrate 
that  low-energy  recoils  that  are  directly  generated  by  heavy  ion  strikes  on  oxide  layers  can  produce 
low -resistivity  paths  through  the  oxide  layers,  resulting  in  appreciable  leakage  currents.  [M.  J. 
Beck,  Y.  S.  Puzyrev,  N.  Sergueev,  K.  Varga,  R.  D.  Schrimpf,  D.  M.  Fleetwood,  and  S.  T.  Pantelides, 
“The  role  of  atomic  displacements  in  ion-induced  dielectric  breakdown”,  IEEE  Trans.  Nucl.  Sci. 
vol.  56,  No.  6,  pp.  3210-3217  (2009)]. 
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Gate  Voltage  (V) 

Fig.  87.  Calculated  defect-mediated  leakage  current  vs.  voltage  (red  curve),  overlain 

on  measured  post-irradiation  leakage  currents,  labeled  by  the  gate  voltages  during 
irradiation.  The  curves  labeled  -3.2  to  -3.9  V  correspond  to  RISB. 

Physical  Mechanisms  Responsible  for  the  Abruptness  of  the  Si-Si02  Interface 

We  used  extensive  first-principles  density-functional  theory  calculations  to  identify  key  mechanisms 
that  control  the  morphology  of  the  Si-Si02  interface.  The  key  findings  that  support  a  smooth  layer- 
by-layer  propagation  of  oxidation  fronts  are  the  following:  i)  there  is  an  effective  barrier  of  about 
2  eV  for  02  molecules  to  cross  the  near  interfacial  region  and  arrive  at  the  substrate.  This  barrier 
accounts  for  the  observed  2-eV  activation  energy  for  the  “reaction”  step  in  the  Deal-Grove  model 
of  thermal  oxidation,  ii)  the  oxidation  reaction  of  Si-Si  bonds  through  dissociation  of  02  molecules 
is  almost  barrier-less  and  strongly  exothermic.  This  fact  alone  points  to  a  random  deposition  (RD) 
problem  with  large  roughness,  contrary  to  the  established  smoothness  of  the  Si-Si02  interface,  iii) 
diffusion  of  O  atoms  between  Si-Si  bonds  along  the  oxidation  front  is  activated  at  thermal  oxidation 
temperatures  Tth  due  to  relatively  small  barriers  of  2-2.3  eV.  iv)  in  certain  cases,  interfacial  Si-O-Si 
bridges  agglomerate  with  binding  energies  of  about  0.3-0 .4  eV  per  O  atom.  In  combination  with 
finding  iii),  we  thus  uncover  the  basic  condition  that  favors  possible  nucleation  and  island  growth 
mode,  v)  migration  of  O  in  bulk  Si  is  also  activated  at  Tth,  which  raises  the  question  why  O  does 
not  simply  dissolve  in  Si.  As  a  mechanism  that  suppresses  solvation,  we  find  that  errant  O  species 
that  arrive  from  the  substrate  get  trapped  at  the  interface  and  eventually  relax  in  configurations 
consistent  with  a  smooth  front  (an  example  of  the  latter  relaxation  is  shown  in  Fig.  88).  Occasional 
desorption  back  into  the  bulk  Si  decreases  the  long-range  roughness  of  the  interface.  Thus,  the  results 
described  in  iii)-v)  resolve  the  quandary  posed  in  ii)  and  provide  a  description  for  the  evolution 
of  the  oxidation  front  as  a  random  “deposition”  process  with  relaxation.  This  kind  of  process  is 
known  to  result  in  smooth  film  growth  on  surfaces  and  also  provides  for  overall  smoothness  in  the 
case  of  the  Si-Si02  interface,  as  observed.  [L.  Tsetseris  and  S.  T.  Pantelides,  “Oxygen  migration, 
agglomeration,  and  trapping:  key  factors  for  the  morphology  of  the  Si-Si02  interface”,  Phys.  Rev. 
Lett.  97,116101  (2006)]. 
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Fig.  88.  Sticking  and  reorganization  (smoothing)  of  an  (a)  O  dimer,  (barrier  is  2.2  eV),  and  an  (b)  O 
trimer  at  the  Si-Si02  interface. 


D.  Properties  and  Phenomena  in  Strained-Si  Channels 

First-principles  Calculations  of  Electron  Mobilities  in  Strained-Si  channels 

We  implemented  a  first-principles,  parameter- free  approach  based  on  density  functional  theory 
to  calculate  electron  mobilities  in  Si  channels.  Initial  calculations  were  based  on  SOI  models  for 
which  atomic-scale  Si02-Si-Si02  structures  were  constructed.  Atomic-scale  roughness  at  the 
interface,  consisting  of  either  suboxide  bonds  (Si-Si  bonds  on  the  oxide  side  of  the  interface) 
or  oxygen  protrusions  (Si-O-Si  bonds  on  the  Si  side  of  the  interface).  We  applied  the  method  to 
explore  the  role  of  strain  in  enhancing  the  electron  mobility.  The  origin  of  this  enhancement  has 
not  been  clearly  established  on  the  basis  of  the  usual  mechanisms  (phonons,  Coulomb  scattering, 
and  macroscopic  interface  roughness).  We  found  that  scattering  from  the  primitive  atomic-scale 
roughness  elements  at  the  Si-Si02  interface  can  account  for  the  data,  as  shown  in  Fig.  89.  The  small 
open  circles  represent  strain-induced  mobility  enhancement  when  the  only  interface  scattering 
centers  are  suboxide  bonds;  the  open  squares  represent  strain-induced  mobility  enhancement  when 
the  only  interface  scattering  centers  are  oxygen  protrusions.  The  large  open  circles  represent  the 
combined  effect.  Black  symbols  are  data  from  various  references  taken  in  single-gate  devices.  [G. 


Strain  (%) 
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Fig.  89.  Calculated  electron  mobilities  in  a  strained-Si  channel  compared  with 
available  data. 
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Hadjisavvas,  L.  Tsetseris,  and  S.  T.  Pantelides,  “The  origin  of  electron  mobility  enhancement  in 
strained  MOSFETs”,  IEEE  Electron  Device  Letters,  vol.  28,  No.  11,  p.  1018-1020  (2007)]. 

Hydrogen-dopant  Interactions  in  SiGe  and  Strained  Si 

The  appearance  of  carrier  traps  and  the  deactivation  of  dopants  are  typical  hydrogen-related 
phenomena  that  are  of  prime  importance  to  the  reliability  of  traditional  Si-based  devices.  We  used 
density  functional  calculations  to  probe  the  dynamics  of  hydrogen  as  individual  impurities  or  in 
complexes  with  dopants  in  strained  Si  (s-Si)  and  SiGe  systems.  We  found  that  the  charged  state 
determines  the  tendency  of  hydrogen  to  be  released  from  dopant  sites  and  to  shuttle  between  a  SiGe 
substrate  and  a  s-Si  overlayer.  In  this  way,  the  effect  of  hydrogen  differs  between  accumulation  and 
inversion  cycles  of  s-Si  and  SiGe  devices.  [L.  Tsetseris,  D.  M.  Fleetwood,  R.  D.  Schrimpf,  and  S. 
T.  Pantelides,  “Hvdrogen-dopant  interactions  in  SiGe  and  strained  Si”  Appl.  Phys.  Lett.  96, 251905 
(2010)]. 

Vacancies  in  SiGe/Strained-Si  Structures 

In  other,  still  unpublished  work,  we  examined  the  relative  energetic  of  Si  vacancies  in  strained-Si 
and  a  SiGe  substrate.  We  found  that  Si  vacancies  and  Si  self-interstitials  have  lower  formation 
energies  in  SiGe,  which  has  desirable  consequences  for  radiation-induced  defects:  Si  vacancies  and 
Si  self-interstitials  would  tend  to  leave  the  strained-Si  channel  and  migrate  into  the  SiGe  substrate. 
A  secondary  effect  is  that  Si  interstitials  in  SiGe  tend  to  replace  a  network  Ge  and  create  a  Ge 
interstitial.  Another  secondary  effect  is  that  Ge  atoms  tend  to  agglomerate  around  Si  vacancies  and 
Ge  interstitials  in  the  SiGe  substrate. 

E.  Theory  of  Ge-Ge02  Interfaces 

Defect  Passivation  and  the  Morphology  of  the  Ge-Ge02  Interface 

Si  has  been  the  electronic  material  of  choice  because  of,  among  other  reasons,  its  high-quality 
interface  with  its  oxide  compound  Si02.  But  as  high-k  materials  replace  Si02  in  the  gate  dielectric, 
new  materials  emerge  also  to  take  the  place  of  Si  as  the  active  semiconductor  in  devices.  Ge  is  a 
possible  alternative  because  of  its  large  carrier  mobility.  There  is  currently  a  resurgence  in  efforts 
to  grow,  analyze,  and  optimize  Ge-based  devices,  especially  those  that  combine  Ge  with  high-k 
systems,  such  as  Hf02.  Typically,  a  thin  GeOx  interlayer  appears  between  the  Ge  substrate  and 
Hf02  the  properties  of  which  has  attracted  significant  experimental  and  theoretical  interest  in  the 
last  couple  of  years.  One  of  the  most  important  issues  for  this  matter  is  the  existence  and  nature 
of  interfacial  carrier  traps.  In  the  case  of  Si,  Pb  centers  (i.e.  interfacial  Si  dangling  bonds)  are  well 
established  as  the  dominant  interface  traps.  There  is  an  ongoing  debate  whether  Pb  centers  are 
important  also  in  the  case  of  Ge-Ge02  interfaces.  Assuming  that  Ge  dangling  bonds  do  exist  at 
the  Ge-Ge02  interface,  the  most  plausible  passivation  scenario  is  through  hydrogenation.  We  have 
found,  however,  that  an  interfacial  Ge-H  bond  dissociates  through  an  exothermic  reaction  that 
releases  H  to  the  Ge02  network  and  0.4  eV.  The  reaction  barrier  of  the  process  depicted  in  Figs. 
90(a)-(b)  is  1.2  eV.  The  analogous  dissociation  reaction  of  Figs.  90(c)-(d)  with  a  F  passivant  is 
endothermic,  but  the  small  reaction  energy  (0.2  eV)  and  barrier  (0.4  eV)  suggest  that  fluorination  is 
not  a  very  effective  way  of  dangling  bond  passivation  either. 

Another  important  issue  for  Ge-based  devices  is  the  morphology  of  Ge-Ge02  interfaces.  In  the  Si- 
Si02  case,  we  previously  identified  a  set  of  mechanisms,  namely  oxygen  migration,  agglomeration, 
and  trapping,  that  support  a  scenario  of  abrupt  and  smooth  oxidation  front.  We  probed  the  same  set 
of  mechanisms  for  Ge-Ge02  boundaries.  We  found  that  an  02  species  migrates  as  peroxy  bridges 
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Fig.  90.  Dissociation  of  passivated  bonds  at  the  Ge-Ge02 
interface:  (a),  (c)  depict  passivated  Ge-H  and  Ge-F  bonds,  (b),  (d) 
show  Ge  dangling  bonds  and  released  H,  F  atoms. 


towards  the  Ge  substrate  where  it  creates  a  pair  of  vicinal  Ge-O-Ge  structures  releasing  about  3 
eV  of  energy.  We  then  investigated  the  stability  of  various  Ge-O-Ge  protrusions  and  found  that, 
compared  to  the  Si-Si02  system,  there  is  a  stronger  tendency  for  atomic-scale  roughness  at  the  Ge- 
Ge02  interface.  For  example,  the  2-0  protrusion  of  Fig.  91(b)  is  more  stable  (by  0.4  eV)  than  the 
“smoother”  structure  of  Fig.  91(a).  Further  studies  will  clarify  if  such  protrusions  are  in  fact  thermal 
donor  precursors  that  can  transform  to  carrier  traps.  [L.  Tsetseris  and  S.  T.  Pantelides,  “Morphology 
and  defect  properties  of  the  Ge-Ge02  interface”,  Appl.  Phys.  Lett.  95,  262107  (2009)]. 


Fig.  91.  O  protrusions  at  a  Ge-Ge02  interface. 


Ge  Volatilization  Products  in  High-k  Dielectrics  on  Ge 

GeO  molecules  are  often  emitted  by  Ge  substrates  under  high-temperature  annealing  and,  in  the  case 
of  gate  stacks,  they  diffuse  through  high -k  oxides.  We  used  first-principles  quantum-mechanical 
calculations  to  probe  the  stability  of  these  impurities  in  La203  and  Hf02  and  their  effect  on  the 
electronic  properties  of  the  host  systems.  We  found  that  the  GeO  species  introduce  several  different 
levels  inside  the  energy  band  gaps  of  La203  and  Hf02.  As  a  result,  the  impurities  may  act  as  charge 
carrier  traps.  Hydrogenation  of  the  GeO  defects  modifies  the  position  and  numbers  of  gap  states, 
but  does  not  eliminate  the  carrier  trap  levels  completely.  The  results  suggest  a  possible  role  of  Ge 
volatilization  in  enhancing  leakage  currents  and  degradation  in  high -k  gate  stacks  of  Ge-based 
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TABLE  1.1 

SIGE  HBT  DEVICE  LIST  FOR  IBICC  TESTING 


Emitter  Area 

Ab  (fin l) 


Resistivity 

(Gt  ■  cm) 


IBM  5 HP  (IBM  5HP) 

Jazz  SiGe-120  Bulk  (Jazz  bulk) 

Jazz  SiGe-120  SOIt  (Jazz  SOI) 

Jazz  SiGe-120  High  Resistivity  Substrate  f  (Jazz  HRS) 
National  SiGe -Si ED  Bulk  (NSC  bulk) 

National  SiGe-SiED  epif  (NSC  epi) 


0.5  x  1O.0 
0.2  X  10.16 
0,2  x  10,16 
0.2  x  10.16 
0.4  x  2O.0 
0.4  x  2ao 


*- 10 
B- 10 
S40 
1500 
S40 
0.009 


t  Not  standard  eomniercial  product  offerings.  Experimental  hardware  lots, 


devices.  [E.  Golias,  L.  Tsetseris,  A.  Dimoulas,  and  S.  T.  Pantelides,  “Ge  volatilization  products  in 
high-k  dielectrics”,  Microelectron.  Engin.,  in  press.] 

3.25  Substrate  Engineering  Concepts  to  Mitigate  Charge  Collection 
in  Deep  Trench  Isolation  Technologies 

A.  Abstract 

Delayed  charge  collection  from  ionizing  events  outside  the  deep  trench  can  increase  the  SEU 
cross  section  in  deep  trench  isolation  technologies.  Microbeam  test  data  and  device  simulations 
demonstrate  how  this  adverse  effect  can  be  mitigated  through  substrate  engineering  techniques. 
The  addition  of  a  heavily  doped  p-type  charge-blocking  buried  layer  in  the  substrate  can  reduce 
the  delayed  charge  collection  from  events  that  occur  outside  the  deep  trench  isolation  by  almost 
an  order  of  magnitude,  implying  an  approximately  comparable  reduction  in  the  SEU  cross  section. 

B.  Overview 

Some  device  technologies  that  typically  make  use  of  Deep  Trench  Isolation  (DTI),  such  as  Silicon- 
Germanium  HeterojunctionBipolarTransistors(  SiGeHBTs),haveSingleEvent  Upset  (SEU)  cross 
sections  exceeding  the  active  area  of  the  device.  Recent  work  reveals  that  shift  registers  composed 
of  a  particular  type  of  SiGe  HBT  have  an  SEU  cross  section  approximately  ten  times  the  active 
area  of  the  device,  which  is  defined  by  the  area  enclosed  by  the  DTI.  This  work  shows  that  the  SEU 
cross  section  enhancement  at  highly  ionizing  energies  is  the  result  of  delayed  (long  time)  charge 
collection  from  regions  of  charge  generation  many  micrometers  away  from  the  DTI.  The  range  and 
magnitude  of  charge  collection  is  primarily  governed  by  the  doping  concentration  and  physical 
structure  of  the  substrate. 

Heavy  ion  microbeam  irradiation  reveals  the  effects  of  substrate  engineering  in  SiGe  HBTs, 
specifically  how  it  affects  delayed  charge  collection.  Based  on  the  microbeam  charge  collection 
data,  a  substrate  engineering  concept  in  the  form  of  a  p-type  charge-blocking  buried  layer  at  the 
bottom  of  the  DTI  is  proposed.  Three  dimensional  Technology  Computer-Aided  Design  (TCAD) 
simulations  reveal  that  the  buried  layer  concept  reduces  charge  collection  by  more  than  70%  for 
lightly  ionizing  events  that  occur  outside  the  DTI  and  by  more  than  a  factor  of  ten  for  highly 
ionizing  events. 
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c. 


Ion  Microprobe  Testing 


All  of  the  SiGe  HBTs  underwent  four-probe  IBICC  so  that  charge  collection  on  all  the  terminals 
(emitter,  base,  collector,  and  substrate)  could  be  monitored  simultaneously.  The  emitter,  base,  and 
collector  were  held  at  0  V  while  the  substrate  was  biased  to  -  5  V  via  the  charge  collection  amplifier 
biasing.  These  charge  collection  bias  conditions  were  intended  to  replicate  the  “off-state”  of  the 
device-under-test  (DUT),  the  most  sensitive  state.  This  charge  collection  setup  had  fidelity  to 
those  conditions.  Experiments  were  conducted  using  normally-incident  36  MeV  oxygen  ions  with 
approximately  a  1  pm  spot  size.  The  ions  have  a  range  of  25.5  pm  in  silicon,  a  surface  incident  LET 
of  5.2  MeV  cm2  /mg, 

IBICC  data  is  reported  for  six  devices  from  three  different  vendors:  IBM  Corporation,  Jazz 
Semiconductor  and  National  Semiconductor,  with  relevant  parameters  detailed  in  Table  I.  The 
abbreviations  used  in  Table  I  will  be  used  throughout  the  rest  of  the  manuscript.  All  devices  were 
subjected  to  the  same  ionizing  radiation  under  identical  bias  conditions.  The  Jazz  HRS  and  SOI 
SiGe  HBT  devices,  as  well  as  the  NSC  epi  SiGe  HBT  device,  supplied  by  BAE  Systems,  were  from 
experimental  hardware  lots  and  are  not  standard  commercial  product  offerings.  All  devices  are  with 
p-type  substrates.  The  DTI  in  all  six  devices  is  approximately  the  same. 

The  data  sets  shown  in  Fig.  92  (a)-(c)  displays  integrated  charge  collection  on  the  collector  terminal, 
since  ion-induced  collector  current  is  the  major  cause  of  SEU  in  most  standard  SiGe  HBT  circuits. 
All  subsequent  plots  of  charge  collection  will  be  for  the  collector  terminal.  While  the  data  plots 
here  focus  on  the  collector  terminal,  the  majority  of  the  current  sensed  on  the  collector  terminal 
also  appears  on  the  substrate  terminal.  The  currents  have  opposite  polarities  since  one  current  is 
traveling  into  the  device  and  one  is  traveling  out  of  the  device. 

The  data  clearly  show  the  structure  of  the  devices.  The  peak  charge  collection  occurs  for  strikes 
within  the  DTI  and  the  tails  represent  charge  collection  from  events  occurring  outside  the  DTI.  These 
two  regions  of  data  reveal  separate  charge  collection  characteristics,  which  are  position  dependent. 
Within  the  DTI,  drift  transport  dominates  due  to  the  extension  of  the  potential  into  the  substrate 
from  the  bottom  of  the  subcollector  (collector-substrate)  depletion  region.  Normally-incident  36 
MeV  oxygen  ions  deposit  about  26  MeV  in  the  substrate  of  a  typical  etched  device,  which  generates 
roughly  1.1  pC  of  charge.  With  the  noted  exception  for  the  IBM  5HP  devices  applies  here,  since 
the  presence  of  the  thick  overlayers  decreased  the  range  of  the  ion  in  the  substrate  and  thus  the  total 
charge  generated. 

The  peak  charge  collection  in  Fig.  92(a)-(c)  occurs  within  the  DTI  for  ions  that  cross  the  reverse- 
biased  subcollector  junction.  The  microbeam  data  show  this  peak  to  be  about  1  pC,  a  charge 
collection  efficiency  of  approximately  90%.  Outside  the  DTI,  there  is  no  pre-existing  electric  field 
to  move  the  charges  once  they  separate  from  the  ion  track,  so  the  charge  collection  is  slower  and 
less  efficient  since  the  electrons  must  diffuse  to  the  subcollector  junction  to  induce  current  on  the 
collector  terminal.  The  tails  to  either  side  of  the  DTI  seen  in  all  the  figures  are  representative  of  this 
collection  mechanism.  They  peak  at  about  200  fC  and  fall  off  to  less  than  50  fC  in  most  cases.  This 
amount  of  charge  is  sufficient  to  cause  upset  in  typical  unhardened  SiGe  HBT  circuits. 

D.  Charge  Collection  Modeling  and  Mitigation  Strategies 

Present  simulation  and  experimental  results  suggest  that  delayed  charge  collection  from  ionizing 
events  outside  the  DTI  can  be  mitigated— if  not  eliminated— by  introducing  a  charge-  blocking 
buried  layer  of  heavily  doped  p-type  silicon  at  the  level  of  the  bottom  of  the  DTI.  A  cross  section 
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Fig.  92.  IBICC  data  for  the  six  SiGe  HBT  DUTs.  The  date  sets  shown  display 
charge  collection  on  the  collector  terminal.  The  deep  trench  isolation  boundary  is 
indicated  approximatley  by  vertical  dashed  lines. 
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of  the  IBM  5HP  device  with  the  addition  of  the  buried  layer  is  shown  in  Fig.  93.  The  p-type  layer 
shown  at  the  level  of  the  bottom  of  the  DTI  is  approximately  2  mm  thick  .  First  generation  SiGe 
HBT  devices  with  this  buried  layer  could  be  manufactured  with  little  to  no  degradation  in  device 
performance  since  the  retrograde  well  is  more  than  4  mm  below  the  subcollector  junction. 

To  evaluate  the  effectiveness  of  the  proposed  buried  layer,  a  series  of  DESSIS  simulations  was 
conducted.  As  evidenced  by  the  simulation  points  shown  in  Fig.  94,  the  buried  layer  reduces  charge 
collection  from  events  outside  the  DTI  to  the  level  measured  in  a  SiGe  HBT  with  a  much  higher 
substrate  doping  concentration.  The  buried  layer  does  not  affect  the  charge  collection  magnitude 
for  ionizing  events  that  occur  within  the  active  region  between  the  DTI  walls,  but  it  does  result  in 
a  large  decrease  in  the  sensitive  area. 


Subcol  lector 


i 


l 


Charge-blocking  Buried  Layer 


Fig.  93.  Cross  section  of  the  buried  layer  concept  for  the  IBM  5HP  device. 


E.  Conclusion 

Substrate  engineering  is  a  valuable  tool  for  controlling  delayed  charge  collection  from  events  that 
occur  outside  the  deep  trench  in  DTI  technologies,  like  SiGe  HBTs.  Moderate-to-low  substrate 
doping  concentrations  permit  large  numbers  of  electrons  to  diffuse  hundreds  of  micrometers, 
thereby  increasing  the  charge  collection  efficiency  of  devices  fabricated  on  these  substrates.  This 
effect  increases  the  SEU  cross  section  for  highly  ionizing  particles,  perhaps  an  order  of  magnitude 
or  more  beyond  the  active  region  of  the  device.  The  devices  that  showed  the  least  efficient  charge 
collection  had  the  most  heavily  doped  substrate  or  employed  some  other  means  of  sheltering  the 
sub-  collector  junction— as  in  the  case  of  a  heavily  doped  substrate  or  buried  oxide  layer. 

Based  on  present  simulation  results  and  experimental  conclusions  from  the  NSC  epi  and  Jazz  SOI 
experimental  IBICC  data,  a  device  hardened  with  the  proposed  charge-blocking  buried  layer  would 
be  relatively  insensitive  to  heavy  ion  and  proton  events  that  occur  in  the  substrate  outside  the  DTI, 
without  suffering  the  speed  or  complexity  penalties  of  other  hardening  approaches. 
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Fig.  94.  DESSIS  charge  collection  simulations  of  the  IBM  5HP  device  with  three 
different  versions  of  the  p-type  charge-blocking  buried  layer. 

3.26  A  Generalized  SiGe  HBT  Single-Event  Effects  Model  for  On- 
Orbit  Event  Rate  Calculations 

A.  Abstract 

This  work  draws  on  experimental  and  simulation  results  to  derive  a  generalized  SEU  response 
model  for  bulk  SiGe  HBTs.The  model  was  validated  using  published  heavy  ion  and  new  proton 
data  gathered  from  high-speed  HBT  digital  logic  integrated  circuits  fabricated  in  the  IBM  5AM 
SiGe  BiCMOS  process.  Calibrating  to  heavy  ion  data  was  sufficient  to  reproduce  the  proton  data 
without  further  adjustment.  The  validated  model  is  used  to  calculate  upset  event  rates  for  low-earth 
and  geosynchronous  orbits  under  typical  conditions. 

B.  Overview 

The  generalized  model  developed  in  this  work  was  applied  to  two  different  DFF  shift  register 
designs  fabricated  in  the  IBM  5AM  SiGe  BiCMOS  (IBM5AM)  process:  a  baseline  design  and  a 
radiation  hardened  by  design  (RHBD)  variant.  This  process  is  characterized  by  a  0.5  pin  drawn 
emitter  width,  a  unity-gain  cutoff  frequency  of  50  GHz,  and  a  of  3.3  V. 

Of  the  two  IBM  5AM  designs  considered  here,  one  was  a  baseline,  nominal  switching  current,  DFF 
shift  register  design  and  the  other  employed  a  RHBD  dual-interleaving  technique  that  included 
duplicated  pass  and  storage  cells,  which  effectively  decoupled  the  differential  inputs  and  outputs 
in  the  storage  cell.  This  input/output  decoupling  increased  the  critical  charge  of  this  design.  For 
the  sake  of  simplicity,  these  designs  are  referred  to  as  “baseline  design”  and  “RHBD  design” 
throughout.  Both  shift  register  designs  are  127-bits  long  and  were  fabricated  solely  out  of  IBM 
5AM  SiGe  HBTs;  no  complementary  metal-oxide-semiconductor  (CMOS)  transistors  were  used. 

The  primary  broadbeam  heavy  ion  data  set  upon  which  this  work  is  based  is  shown  in  Fig.  95.  It 
is  obvious  that  with  the  low  FET  neon  ion  the  cross  section  of  the  RHDB  design  does  not  increase 
with  increasing  angle,  but  instead  decreases  with  increasing  angle— i.e.,  decreasing  cross  section 
with  increasing  effective  FET.  This  behavior  violates  the  assumptions  of  the  original  rectangular 
parallelepiped  (RPP)  model,  which  generally  assumes  increasing  cross  sections  with  increasing 


MURI  FINAL  REPORT  •  2005  -  2010  •  VANDERBILT  UNIVERSITY 

PAGE  75 


Fig.  95  Broadbeam  heavy  ion  data  for  the  baseline  and  RHBD  127-stage  shift  register  designs.  The 
important  feature  is  the  cross  section  decrease  with  increasing  angle  for  the  RHBD  device  with  a 

higher  critical  charge. 

effective  LET.  Since  these  data  are  not  described  adequately  by  the  default  RPP  model,  they  have 
been  re-plotted  with  the  RPP  cosine  corrections  removed.  The  data  are  plotted  as  a  function  of  angle 
instead  of  effective  LET  and  the  cross  section  was  scaled  by  cos(theta)  to  remove  the  effective 
fluence  correction.  All  subsequent  data  sets  will  be  plotted  in  this  manner  to  avoid  confusion.  For 
the  sake  of  reference,  the  normally-incident  LETs  for  the  ions  in  Fig.  95  are  Ne  2.8  MeV  cm2/mg, 
Ar  8.3  MeV  cm2/mg,  and  Xe  53  MeV  cm2/mg. 

Since  a  large  number  of  space  borne  particles  are  incident  at  oblique  angles,  understanding  the 
angular  response  of  bulk  SiGe  HBTs  is  critical  to  developing  a  representative  rate  prediction  model. 

C.  Energy  Deposition  Response  Model 

This  work  relied  on  the  energy  transport  and  calorimetry  capabilities  of  the  MRED  tool  set.  Using 
this  tool,  it  is  possible  to  compute  the  energy  deposited  in  one  or  more  sensitive  (fiducial)  volumes 
due  to  impinging  ions.  Furthermore,  these  fiducial  volumes  can  have  weights.  The  volumes  and 
their  weights  function  in  an  ensemble  to  form  a  linear  combination  that  approximates  the  total 
collected  charge.  The  approach  is  described  by  equation  2.1.  The  total  collected  charge  is  the  sum 
over  all  fiducial  volumes  of  the  product  of  the  weight  and  total  charge  liberated. 

The  total  charge  liberated  (Q)  is  related  to  the  total  energy  deposited  (ED)  through  the  relationship 
Q=lpC/22.5  MeV  '  ED  .  This  linear  combination  of  weighted  fiducial  volumes  is  the  construct  that 
will  be  used  to  model  the  energy  deposition  response  of  the  SiGe  HBTs  considered  in  this  work. 
Once  calibrated  to  data,  usually  heavy  ion  broadbeam  cross  section  data,  this  modeling  method 
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provides  an  accurate,  high-speed  approximation  to  the  initial  conditions  and  ensuing  temporal 
evolution  of  charge  transport  and  collection. 

A  2-D  projection  of  the  basic  energy  deposition  response  model  is  shown  in  Fig.  96.  The  fiducial 
volumes  have  been  overlaid  on  the  TCAD  cross  section.  The  top-down  area  of  each  volume  is 
determined  by  the  normal-incident  cross  section  of  the  broadbeam  heavy  ion  data,  which  includes 
Ne,  Ar,  and  Xe.  The  weights  and  depths  of  each  of  the  volumes  are  calculated  by  correlating 
microbeam  data  to  TCAD  simulations.  This  model  was  used  in  a  fully  reconstructed  3-D  model  of 
a  shift  register  stage  for  all  subsequent  simulations  and  the  event  rate  calculations. 


Fig.  96.  This  is  the  basic  weighted  fiducial  volume  ensemble  used 
to  model  the  radiation  response  of  the  IBM  5AM  SiGe  HBT 
process  considered  in  this  work. 


D.  Heavy  Ion  Response 

Before  computing  on-orbit  event  rates,  the  model  was  verified  against  data  sets  that  covered 
enough  of  the  possible  response  parameter  space  to  ensure  predictable  behavior  in  a  more  diverse 
environment  such  as  geosynchronous  or  low-earth  orbit.  The  model  was  calibrated  to  the  heavy  ion 
datasets  for  the  baseline  and  RHBD  designs. 

The  calibrated  heavy  ion  results  for  both  the  baseline  and  RHBD  circuit  designs  are  shown  in  Fig. 
97(a)  and  (b).  The  critical  charge  for  each  of  the  calibrations  is  displayed  on  the  individual  figures. 

E.  Proton  Calibration  Results 

Proton  data  were  collected  during  experiments  on  the  CREST  chip  at  the  Crocker  Nuclear 
Laboratory  (CNL)  at  the  University  of  California  at  Davis  using  63  MeV  protons.  These  data  were 
taken  on  the  baseline  and  RHBD  designs  at  several  different  data  rates  at  normal  incidence  and  a 
grazing  angle.  The  full  data  set  is  plotted  in  Fig.  98(a).  This  data  set  can  now  be  used  to  check  the 
proton  response  of  the  model  developed  in  the  previous  section.  These  modeling  results  are  shown 
in  Fig.  98(b).  The  heavy  ion  model  was  validated  against  proton  data  by  only  a  change  of  particle 
and  energy  in  the  simulation  environment.  The  strong  match  between  simulation  and  data  validates 
a  larger  portion  of  the  model’s  acceptable  parameter  space,  making  it  usable  for  environments  with 
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Fig  97.  Plots  (a)  and  (b)  show  the  calibrated  results  of  the  computer  simulations  for  the 
entire  127-stage  shift  register.  In  each  case  the  open  symbols  are  the  date  from  Fig.  95 
and  the  closed  symbols  are  the  derived  simulation  resits,  (a)  Baseline  design  heavy  ion 
calibration;  (b)  RHBD  heavy  ion  calibration. 


Fig.  98.  (a)  63  MeV  proton  data  from  CNL  taken  on  the  baseline  and  RHBD  CREST  chip 
designs,  (b)  A  comparison  between  the  proton  data  from  (a)  and  the  simulation  results 
based  on  the  model,  which  was  calibrated  to  heavy  ion  data.  The  strong  agreement  shown 
above  was  obtained  by  simply  changing  the  particle  and  energy  int  he  simulation;  no  further 
adjustments  were  made.  The  open  symbols  are  data  based  ont  he  average  across  data  rate 
from  (a)  and  the  solid  symbols  are  simulation  resultes  using  the  response  model. 

large  proton  fluxes.  These  simulation  results  were  obtained  from  the  model  calibrated  with  heavy 
ion  data  only;  no  further  adjustment  was  required. 

F.  Event  Rate  Calculation 

Now  that  we  have  a  calibrated  model,  we  can  use  MRED  to  predict  the  failure  rate  in  any  space 
environment .  An  example  of  this  is  shown  I  Fig  99 ,  were  flux  spectra  for  the  geosynchronous  (GEO) 
and  low-earth  orbit  (LEO)  environments.  Both  of  the  environments  were  solar  minimum/quiet 
conditions,  included  all  available  ion  species,  and  assumed  100  mil  of  aluminum  shielding.  The 
LEO  spectra  were  for  the  space  station  orbit,  which,  according  to  CREME96,  is  at  an  inclination  of 
51 .6  degrees  and  an  orbital  radius  of  500  km. 

G.  Conclusions 

The  shape  and  relationship  between  the  fiducial  volumes  represents  a  critical  aspect  of  this  study. 
The  fiducial  volumes  explain  in  a  quantitative  and  qualitative  way  the  once-anomalous  angular 
response  of  the  technology,  the  low  SEU  threshold,  and  the  large  saturated  cross  section  observed  for 
the  most  highly  ionizing  particles.  Putting  this  model  together  unifies  many  years  of  experimental 
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Fig.  99.  These  figures  show  the  simulated  event  rates  for  the  two  shift  reiter  designs.  The 
data  markers  have  been  thinned  to  aid  viewing.  The  large  marekrs  show  where  the  rate  was 
evaluated  based  on  the  derived  critical  charge,  (a)  Geosynchronous  orbit  event  rate;  (b)  low- 

earth  orit  event  rate. 

and  theoretical  work  and  provides  intuition  to  designers  considering  SiGe  HBT  projects.  It  also 
opens  the  door  for  more  complex  types  of  modeling  that  will  begin  to  look  at  time-domain  effects 
and  other  aspects  of  extreme,  high-speed  digital  technologies  like  this. 

3.27  Laser-Induced  Current  Transients  in  Silicon-Germanium  HBTs 

A.  Abstract 

Device-level  current  transients  are  induced  by  injecting  carriers  using  two-photon  absorption  from 
a  subbandgap  pulsed  laser  and  recorded  using  wideband  transmission  and  measurement  equipment. 
These  transients  exhibit  three  distinct  temporal  trends  that  depend  on  laser  pulse  energy  as  well  as  the 
transverse  and  vertical  charge  generation  location.  The  nature  of  the  current  transient  is  controlled 
by  both  the  behavior  of  the  subcollector-substrate  junction  and  isolation  biasing.  However,  substrate 
potential  modulation,  due  to  deformation  of  the  subcollector-  substrate  depletion  region,  is  the 
dominant  mechanism  affecting  transient  characteristics. 

B.  Introduction 

In  this  work  we  examine  two-photon  absorption  current  transients  from  a  passively-  loaded  single 
SiGe  HBT  with  the  intent  of  providing  accurate  feedback  for  analogous  TCAD  simulations.  The 
experimental  results  presented  here  were  obtained  from  a  SiGe  HBT  test  structure  in  the  IBM  5  AM 
BiCMOS  process  technology.  There  are  three  substrate  taps  for  this  device,  located  approximately  9 
mm  from  the  outside  of  the  deep  trench.  All  four  of  the  device  terminals  were  connected  to  individual 
pads  on  the  die.  The  device  under  test  was  mounted  in  a  custom  broadband  40  GHz  package 
designed  for  topside  and  backside  laser  irradiation  as  well  as  heavy  ion  microbeam  irradiation.  All 
four  device  terminals  were  connected  to  40  GHz  bias  tees  with  the  AC+DC  connection  facing  the 
device  and  the  AC-only  connection  facing  the  oscilloscope.  The  experiments  used  a  Tektronix  12 
GHz  TDS6124C  real-time  digital  storage  oscilloscope.  All  transmission  lines  were  50,  2.9  mm  (40 
GHz)  coaxial  assemblies.  The  die  pads  were  bonded  to  the  microstrip  transmission  lines  with  1  mil 
gold  bond  wires  less  than  1  mm  in  length.  Time  domain  reflectometry  measurements  indicated  that 
the  four  bond  wire  inductances  and  bond  pad  capacitances  were  each  approximately  0.2  nH  and 
60  fF. 

The  device  was  characterized  by  injecting  carriers  using  two-  photon  absorption  from  a  subbandgap 
pulsed  laser.  This  technique  generates  significant  densities  of  electron-hole  pairs  in  proximity  of  the 
beam  focal  point,  which  can  be  deterministically  positioned  in  three  dimensions. 
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Experiments  were  centered  around  three  bias  conditions:  Case  1  had  a  bias  of  on  the  substrate  with 
all  other  terminals  grounded,  Case  2  had  a  bias  of  3  V  on  the  collector  with  all  other  terminals 
grounded,  and  finally,  Case  3  had  a  bias  of  on  the  substrate  with  all  other  terminals  grounded. 
Previous  broadbeam  and  single-photon  absorption  pulsed-laser  experiments  on  pure  bipolar 
applications  were  conducted  with  a  substrate  bias.  However,  BiCMOS  applications  call  for  the 
substrate  to  be  grounded,  which  means  the  device  isolation  bias  must  be  applied  through  a  positive 
voltage  on  the  collector.  The  bias  condition  of  on  the  substrate  was  included  to  highlight  any 
differences  observed  when  placing  3  V  on  the  collector  terminal. 

C.  Experimental  Results 

Pulsed  lasers  have  been  used  for  many  years  to  study  single-event  effects  (SEE)  in  microelectronics. 
With  the  application  of  intense  femtosecond  light  pulses,  carrier  generation  by  non-linear  absorption 
becomes  a  useful  way  to  study  SEE. 

The  data  shown  in  Fig.  100  are  from  experimental  Case  1 ,  with  a  bias  of  on  the  substrate.  The  plot 
shows  the  maximum  measured  total  integrated  charge  vs.  peak  pulse  irradiance  on  a  log-log  scale. 
This  is  a  quadratic  relationship,  and  since  the  slopes  are  approximately  equal  to  two,  two-photon 
absorption  is  confirmed  for  the  collector  and  for  the  base.  The  base  is  a  more  narrow  region,  so 
charge  collection  is  less  efficient  and  some  charge  will  be  lost,  making  the  trend  subquadratic. 


Fig.  100.  Verification  of  the  two-photon  absorption  process 
in  these  experiments.  The  data  are  from  experimental 
Case  1 ,  with  -4  V  ont  he  substrate  nd  all  other  terminals 
grounded.  Int  he  ideal  case,  each  line  would  have  a  slope 
of  two  since  the  number  generated  carriers  depends  on 
the  square  of  the  pulse  irradance.  There  is  a  6%  fractional 
error  int  he  pulse  power  and  a  6%  fractional  error  int  he 
focal  spot  size,  which  produces  are  irradiance  error  of  6% 
and  a  charge  collection  errof  of  6%.  The  errors  in  charge 
collection  are  smaller  than  the  data  point. 


The  data  shown  in  Fig.  101  (a)-(d)  show  maps  of  the  peak  transient  magnitude  for  the  emitter, 
collector,  base,  and  substrate  for  the  bias  conditions  of  Case  1 .  The  emitter  and  collector  transients 
are  negative;  the  base  and  substrate  transients  are  positive.  A  schematic  of  the  deep  trench  isolation 
layout  is  overlaid  on  the  data.  The  substrate  taps  for  these  devices  are  located  9  mm  in  the  negative 
-direction  from  the  bottom  of  the  deep  trench  isolation.  The  laser  spot  size  is  approximately  1.6 
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Fig.  101.  These  maps  show  the  maximum  current-transient  magnitude  for  a  5  x.  1012 
W / cm2  peak  irradiance  laser  pise  on  an  IBM  5AM  siGE  HBT  with  bias  conditions  of 

Case  1. 


in  diameter  and  focused  at  the  device  surface  in  the  -dimension.  The  emitter  data  in  Fig.  101(a) 
support  the  two-photon  absorption  mechanism  justified  in  Fig.  100.  The  nonzero  data  points  are 
well  confined  to  the  emitter  stripe,  which  could  not  occur  if  a  significant  amount  of  linear  absorption 
was  present.  If  linear  absorption  were  present,  the  spatial  distribution  would  be  smeared  out  after 
the  beam  passed  through  almost  300  mm  of  silicon  from  the  backside  of  the  wafer  to  the  surface  of 
the  device.  The  peak  response  of  the  base  and  collector  terminals  is  coincident  with  the  emitter-base 
and  base-collector  junctions.  The  peak  response  is  not  located  at  the  base  and  collector  contacts. 


The  data  in  Fig.  102  (a)  through  (d),  show  maps  of  the  peak  transient  magnitude  for  the  base  and 
collector  for  the  bias  conditions  of  Case  2  and  Case  3.  The  collector  transients  produced  by  the  bias 
conditions  in  Case  2  are  larger  than  the  transients  in  Case  3  for  hits  to  the  active  device  region, 
defined  by  the  emitter-base  and  base-collector  junctions.  While  the  base  transients  in  Case  2  are 
slightly  larger  than  those  produced  by  Case  3,  by  up  to  25%  for  the  peak  response,  the  difference  is 
much  less  than  the  possible  2  []  increase  in  peak  response  for  the  collector  in  Case  2  at  a  peak  pulse 
irradiance  of  9  []  1012  W/cm2  .  The  profile  associated  with  the  base-collector  junction  in  Case  2  is 
also  better  resolved  than  in  Case  3  where  the  response  is  smoother  across  the  active  region  of  the 
device.  As  in  the  case  for  the  data  shown  in  Fig.  102  (a)-(d),  the  laser  was  focused  at  the  surface  of 
the  device  to  a  diameter  of  approximately  1 .6  pm. 


The  data  shown  in  Fig.  103  (a)  and  (b)  display  current  transient  waveforms  for  Case  2  and  Case 
3  at  different  pulse  irradiance.  As  with  the  previous  data  sets,  the  laser  was  focused  at  the  surface 
of  the  device  to  a  diameter  of  1 .6  pm  .  As  shown  in  the  spatially  correlated  maps,  there  is  a  small 
difference  between  the  base  transients,  but  the  Case  2  collector  transients  are  up  to  2  times  larger 
for  the  5  0  1012  W/cm2  and  9  []  1012  W/cm2  traces  for  hits  to  the  base-collector  junction.  In  all  cases, 
the  pulse  lengths  are  approximately  the  same. 


Finally,  the  data  presented  in  Fig.  1 04  show  collector  transients  when  the  focus  of  the  laser  is  moved 
20  pm  below  the  surface  and  scanned  in  the  usual  15  pm  by  15  pm  square.  The  only  difference  in 
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Fig.  102.  These  maps  show  the  maximum  current- transient  magnitude  for  a  5  x.  1012W/ 
cm2  peak  irradiance  laser  pise  on  an  IBM  5AM  siGE  HBT  with  bias  conditions  of  Case  2 

and  Case  3.. 
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Fig.  103.  Base  and  collector  current  transients  for  (a)  Case  2  and  (b)  Case  3. 


the  data  are  that  the  response  inside  the  deep  trench  isolation  has  to  be  measured  on  a  50  mV/div 
vertical  scale  and  the  other  two  traces  were  measured  on  a  5  mV/div  scale.  This  difference  explains 
the  noise  floor  discrepancy  in  the  data.  The  most  important  things  about  these  data  are  the  three 
distinct  responses  depending  on  trans-verse  location  relative  to  the  deep  trench  isolation.  When 
far  outside  the  deep  trench  isolation,  a  low-magnitude,  long-duration  diffusion  pulse  is  observed, 
lasting  about  25  ns.  When  inside  the  deep  trench  isolation,  the  response  is  similar  in  magnitude 
and  length  to  pulses  measured  at  the  surface  despite  moving  into  the  substrate  by  20  pm.  The  pulse 
measured  at  approximately  1  outside  the  deep  trench  isolation  shows  two  response  modes,  a  prompt 
collapse  of  the  substrate  potential  (rapid  current  induction)  and  a  delayed  diffusion  component, 
in  other  words  a  superposition  of  the  response  from  inside  the  deep  trench  isolation  and  from  far 
outside  the  deep  trench  isolation. 
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Fig.  104.  Collector  current  transients  for  a  5  x  1012  w/cm2.. 


D.  Discussion 

There  are  three  important  results  based  on  the  data  in  the  previous  section:  the  effect  of  placing  the 
device  isolation  bias  on  the  collector  instead  of  the  substrate,  the  measurement  of  significant  device 
response  for  laser  pulses  outside  the  deep  trench  isolation,  and  the  measurement  of  a  bimodal 
collector  response  for  laser  pulses  focused  deep  in  the  substrate  and  outside  the  deep  trench  isolation. 

The  difference  in  the  collector  isolation  bias  case,  and  the  reason  the  collector  current  transients  are 
up  to  2  times  larger,  is  because  of  base  width  modulation  related  to  the  Early  effect  and  the  onset 
of  avalanche  multiplication. 

Previous  broadbeam  heavy  ion  data  sets  from  irradiating  shift  register  chains  implied  that  sensitive 
device  cross  sections  were  approximately  an  order  of  magnitude  larger  than  the  active  device  area. 
This  pointed  to  the  fact  that  upsets  occurred  for  events  originating  outside  the  deep  trench  isolation, 
even  for  normally  incident  ions.  However,  whether  or  not  events  outside  the  deep  trench  isolation 
could  be  responsible  for  upsets  has  been  debated.  These  data  should  help  settle  that  discussion 
concerning  specific  circuit  applications  considering  the  abundance  of  events  measured  in  excess  of 
1  mA  for  pulses  located  outside  the  deep  trench  isolation. 

Thirdly,  and  perhaps  most  important  for  device  physics  and  upset  mechanism  perspectives,  is  the 
issue  of  diffusion-induced  potential  modulation.  Charge  funneling,  as  conceived  in  the  early  1980s, 
successfully  described  localized  potential  deformation,  where  the  transverse  area  disturbed  by  the 
incident  particle  was  much  smaller  than  either  the  junction  crossed  or  the  affected  device.  However, 
most  bulk  SiGe  HBTs,  like  the  IBM  5AM  process,  are  manufactured  on  lightly  doped  p-type 
substrates,  typically  around  1  ,  which  means  that  when  potentials  deform  into  the  substrate  of  these 
devices,  the  potential  modulation  covers  a  very  large  area  in  order  to  maintain  the  potential  drop 
fixed  by  the  applied  voltage  at  either  the  substrate  or  collector  contact.  The  potential  modulation 
is  no  longer  a  localized  phenomenon,  but  a  delocalized  effect  that  does  not  require  the  ion  to  cross 
the  active  junction  and  can  easily  span  lateral  and  vertical  distances  in  excess  of  10  um.  For  this 
reason,  potential  modulation  effects  in  SiGe  HBTs  should  not  be  confused  with  classical  funneling 
processes.  This  substrate  potential  modulation  effect  can  be  observed  in  SiGe  HBTs  with  a  pulsed 
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laser  when  the  focal  point  is  below  the  surface  of  the  device.  When  the  carriers  are  generated  in  the 
field- free  region  deep  in  the  substrate,  outside  the  deep  trench  isolation,  they  diffuse  isotropically. 
If  a  sufficient  density  of  carriers  can  reach  the  subcollector-substrate  junction,  they  will  compensate 
the  ionized  dopants  and  force  the  potential  contour  surfaces  to  expand  into  the  substrate,  initially 
limited  by  the  deep  trench  isolation,  until  the  density  of  free  carriers  becomes  low  enough  that  the 
equilibrium  position  of  the  potential  contours  can  support  the  applied  bias.  The  number  of  carriers 
that  reach  the  junction  is  determined  by  the  amount  of  solid  angle  the  junction  occupies.  This  also 
means  that  the  effect  is  both  location  and  pulse  power  dependent.  In  effect,  the  potential  pushes  out 
from  its  equilibrium  position,  creates  a  large  volume  of  near  equipotential  moving  through  charges 
in  the  substrate,  which  induces  a  prompt  current  on  the  collector  terminal.  The  bimodal  collector 
transient  response  in  Fig.  105  is  the  first  direct  evidence  of  diffusion-initiated  potential  modulation 
in  SiGe  HBTs.  Since  the  carriers  are  generated  a  long  distance  from  the  critical  junction,  the 
potential  expands  and  collapses  before  all  the  charge  is  collected,  so  there  is  a  secondary  diffusion 
current  that  is  realized  after  the  prompt  potential  expansion  and  collapse,  which  yields  the  resulting 
superposition  of  potential  collapse  and  diffusion  responses. 

E.  Conclusion 

The  results  presented  here  confirm  the  presence  of  nonlocal  potential  modulation  in  IBM  5AM 
SiGe  HBTs.  While  this  effect  has  demonstrable  consequences  for  bulk  SiGe  HBTs  in  general,  it 
is  also  expected  to  apply  to  many  critical  situations  in  CMOS,  such  as  well  potential  collapse  and 
recovery.  This  work  has  also  established  that  the  isolation  bias  terminal  could  be  a  critical  aspect 
when  considering  the  peak  magnitude  of  collector  current  transients  and  that  certain  applications 
may  suffer  larger  collector  current  transients  for  events  that  strike  the  active  device  regions  directly. 


2.50 


Fig.  105.  Forced  ouput  characteristics  for  a  IBM  5AM  SSiGe  HBT  under  high  injection. 
The  experimental  bias  condition  of  Case  2  is  marked. 
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3.28  Heavy  Ion  Microbeam-and  Broadbeam-lnduced  Transients  in 
SiGe  HBTs 

A.  Abstract 

Silicon-germanium  heterojunction  bipolar  transistor  (SiGe  HBT)  heavy  ion-induced  current 
transients  are  measured  using  Sandia  National  Laboratories’  microbeam  and  high-and  low-energy 
broadbeam  sources  at  the  Grand  Accelerateur  National  d’lons  Lourds,  Caen,  France,  and  the 
University  of  Jyvaskyla,  Finland.  The  data  were  captured  using  a  custom  broadband  IC  package 
and  real-time  digital  phosphor  oscilloscopes  with  at  least  16  GHz  of  analog  bandwidth.  These  data 
provide  detailed  insight  into  the  effects  of  ion  strike  location,  range,  and  LET. 

B.  Introduction 

This  work  presents  heavy  ion  microbeam  position-correlated  data,  coupled  with  a  range  of 
broadbeam  energies  and  linear  energy  transfers  (LETs),  providing  detailed  device-level  data  on 
the  temporal  profile  of  ion-induced  current  transients  in  this  important  semiconductor  technology. 
The  different  LETs  and  particle  energies  show  the  consequences  of  heavy  ion  charge  generation 
and  collection  in  devices  with  lightly-doped  substrates.  These  data  capture  essential  information 
required  for  accurate  device  physics  modeling. 

We  collected  single-event  current  transients  at  three  different  facilities:  Sandia  National  Laboratories  ’ 
(SNL,  in  Albuquerque,  NM)  Ion  Beams  Materials  Research  Lab  using  a  6  MV  EN  tandem  Van  de 
Graaff  microbeam,  the  Department  of  Physics  at  the  University  of  Jyvaskyla  (JYFL,  in  Finland) 
using  a  K-130  cyclotron,  and  the  Grand  Accelerateur  National  d’lons  Lourds  (GANIL,  in  Caen, 
France)  high-energy  beamline.  The  microbeam  data  gathered  at  SNL  are  based  on  36  MeV  O  and 
include  relative  -coordinates  of  each  ion  strike  based  on  the  microbeam  coordinate  system.  The 
data  collected  at  JYFL  and  GANIL  are  broadbeam  data,  gathered  without  knowledge  of  ion  strike 
location,  but  possess  higher  energy,  the  possibility  of  angled  irradiation,  and  a  wide  selection  of 
LETs.  JYFL  heavy  ion  exposures  include  9.3  MeV/u  Ne,  Ar,  Kr,  and  Xe.  The  Ar  irradiations  at 
JYFL  were  performed  at  a  tilt  of  60  in  addition  to  normal  incidence.  We  performed  irradiations  at 
GANIL  using  45.5  MeV/u  Xe.  We  conducted  all  exposures  at  normal  incidence  unless  otherwise 
noted. 

The  device  under  test  (DUT)  is  an  IBM  5AM  SiGe  HBT  At  JYFL  and  GANIL,  the  transients  on  the 
base  and  collector  were  measured  and  recorded  with  a  Tektronix  DPO71604A  16-GHz  (40  GS/s), 
real-time  digital  phosphor  oscilloscope  (DPO).  At  SNL,  substrate,  collector,  and  base,  transients 
were  measured  and  recorded  with  a  Tektronix  DP072004  20-GHz  (50  GS/s),  real-time  DPO.  The 
oscilloscope  triggered  on  the  collector  channel  for  all  experiments..  The  experiments  focused  on 
three  bias  conditions  for  the  DUT:  (Case  1)  Vsub  =  -4,  (Case  2)  Vc  =  3  V,  and  (Case  3)  Vc  =  -  3  V  If 
the  terminal  is  not  listed,  it  is  grounded. 

C.  Results 

Fig.  106  (a)  and  (b)  plot  the  peak  base  and  collector  currents  as  a  function  of  position  obtained 
from  a  36  MeV  O  time-resolved  ion  beam-induced  charge  (TRIBIC)  scan  on  an  IBM  5AM  SiGe 
HBT  under  the  bias  conditions  of  Case  1 .  The  scan  area  is  20  pm  []  20  pm  with  200  nm  steps  and 
a  spatial  resolution  of  <  1  pm.  The  scans  produced  approximately  400  data  points  based  on  a  mV 
trigger  on  the  collector.  The  peak  collector  responses  are  confined  to  the  base-collector  junction. 
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Fig.  106.  MeV  O  TRIBIC  scan  on  a  IBM  5AM  SiGe  HBT  Bias  for  Case  1 
and  all  other  terminals  grounded. 

The  oscilloscope  triggers  on  the  collector  channel  only  when  the  ion  strikes  in  or  near  the  area 
enclosed  by  the  deep  trench  isolation.  Coupling  this  microbeam  information  with  the  peak  transient 
current  and  integrated  transient  charge  enables  position  correlation  of  the  broadbeam  strikes  in 
reference  to  the  DUT’s  physical  structures. 


In  Fig.  107  (a),  with  a  3  V  bias  on  the  collector  (Case  2),  instead  of  a  -3  V  bias  on  the  substrate  (Case 
3)  as  shown  in  Fig.  107  (b),  the  collector  current  transients  within  the  base-collector  junction  are 
magnified  by  more  than  a  factor  of  two.  The  magnification  of  the  base-collector  junction  transients 
in  Case  2  is  presumably  due  to  a  combination  of  the  Early  effect  and  avalanche  multiplication.  The 
data  in  each  of  these  figures,  106  (a)  and  106  (b),  represents  the  same  voltage  dropped  across  the 
subcollector  junction,  producing  an  equivalent  depletion  layer.  The  nominal  transient  peak  current 
of  0.5  mA  remains  the  same  for  strikes  that  do  not  cross  the  base-collector  junction. 


Fig.  108  shows  two  of  the  larger  base  and  collector  current  transients  obtained  from  the  JYFL 
broadbeam  heavy  ion  results,  demonstrating  the  significance  of  ion  LET  on  the  production  of 
current  transients.  Based  on  knowledge  of  the  microbeam  data  already  presented,  the  transients 
shown  in  Fig.  6  are  the  result  of  direct  hits  to  the  active  region  of  the  device;  each  pair  shown 
are  correlated  events  from  a  single  ion.  As  expected,  the  xenon  transient  in  Fig.  6  produces  more 
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Fig.  107.  36  MeV  O  TRIBIC  scan  on  an  IBM  5AM  SiGe  HBT  with  Case  2  in  (a)  and  with 

Case  3  in  (b). 

Charge  resulting  in  large  transients  on  both  the  collector  and  base  terminals.  The  plateau  in  the 
xenon  collector  transient,  and  large  amount  of  collected  charge,  is  due  to  the  fact  that  the  device 
terminals  are  tied  to  external  voltage  sources  and  capacitors.  Connecting  the  DUT  to  a  circuit  that 
allowed  the  voltage  to  collapse  under  high  current  draw  would  modify  the  plateau  and  perhaps 
shorten  the  transient.  However,  the  plateau  suggests  some  form  of  saturation,  which  is  likely  due 
to  systematic  effects,  such  as  the  bias  tee  capacitor.  The  neon  transients  are  similar  to  the  SNL 
microbeam  transients  in  Fig.  4  and  compare  well  to  previous  pulsed  laser  testing,  indicating  data 
consistency  and  LET  proportionality. 

The  JYFL  argon  results,  shown  in  Fig.  108  (a)  and  (b),  at  normal  incidence  and  a  60  degrees  tilt 
confirm  that  increasing  the  angle  of  incidence  relative  to  the  device  surface  normal  produces  fewer 
transients,  based  on  an  oscilloscope  trigger  value  of  -15  mV  on  the  collector.  At  normal  incidence 
with  -4  V  on  the  substrate,  the  oscilloscope  captured  50  transients  after  a  fluence  of  3.85 " 107  cm- 
2  .  However,  at  a  tilt  of  60  degrees  with  the  same  bias  conditions,  only  16  events  were  measured 
after  a  fluence  of  1.94 ' 108  cm-2  ,  a  16"  decrease  in  cross  section.  This  result  confirms,  at  the 
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device-level,  the  effect  of  cross-section  decrease  with  increasing  ion  angle  for  low  LET  particles, 
observed  in  many  previous  broadbeam  tests  of  SiGe  HBT  circuit  applications  and  described  via 
TCAD  simulation.  It  is  critical  to  understand  this  effect  in  order  to  calculate  event  rates  for  space- 
based  applications. 


Fig.  108.  JYFL  —  ion  LET  comparison  for  Case  1.  Neon  and  xenon  transients 
caputred  at  JYFL  demonstrating  device  response  to  two  extreme  LETs,  3.7  MeV 
cm2/ mg  for  the  neon  ions  and  60  MeV  cm2/ mg  for  the  xenon  ions;  the  integrated 
charge  is  labeled  next  to  each  curve. 


D.  CONCLUSION 

These  heavy  ion  current  transient  data  represent  a  significant  improvement  to  the  state-of-the-art 
understanding  of  heavy  ion-induced  charge  and  transients  in  SiGe  HBTs  by  using  microbeam  data 
to  position-correlate  broadbeam  single-event  current  transients.  These  results  complete  heavy  ion 
microbeam  charge  collection  and  pulsed  laser  data  sets.  Previous  microbeam  data  only  measured 
collected  charge  and  previous  two-photon  pulsed-laser  transient  measurements  cannot  be  easily 
correlated  to  ion  LET  and  ion-specific  effects.  Taken  as  a  whole  and  coupled  with  device  simulation, 
a  complete  picture  of  charge  generation,  transport,  and  collection  in  bulk  SiGe  HBTs  is  possible. 

3.29  Simulation  of  Single  Event  Effects  in  New  Technologies  and 
Ultra-small  Devices 

A.  Overview 

The  thrust  of  this  activity  has  been  to  develop  the  capability  to  predict  the  single  event  behavior  of 
ultra-small  and  non-silicon  technologies  from  information  about  their  composition  and  structure 
using  first-principles  physics  and  a  knowledge  of  the  radiation  environment.  Central  to  this  effort 
has  been  the  development  of  a  Monte  Carlo  simulation  system  for  analysis  of  energy  deposition  in 
microelectronic  structures.  In  this  section,  we  describe  the  computational  system  that  has  emerged 
and  some  of  the  protocols  that  have  been  developed  to  simulate  ion-induced  single  events.  Through 
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Fig.  109.  Architecture  of  the  Geant4  application  MRED,  showing  both  the  Geant4  core  and  other  components.  The  actual 
program  “mred”  is  a  unix  shell  script  that  invokes  the  Python  executable  run_mred.py.  The  run_mred.py  script  loads  several 
Python  modules  including  G4Core.py,  which  contains  most  of  the  computing  machinery,  which  in  turn  resides  in  the 
_G4Core.so  shared  library.  The  _G4Core.so  library  is  built  by  the  Python  distutils  routine  from  the  components  to  the  left  in 

this  illustration  using  a  custom  setup.py  script  [2] . 


this  research  and  related  programs  a  comprehensive  and  unifying  mathematical  framework  for 
single  event  rate  calculation  has  been  developed.  This  framework,  which  encompasses  not  only 
the  Monte  Carlo  methods  that  are  the  focus  of  this  program,  but  also  the  analytical  methods  that 
preceded  it,  and  the  methods  of  applying  it  are  described  in  detail  in  a  recent  review.  Here,  we  give 
an  overview  of  software  and  methods.  Applications  are  described  elsewhere  in  this  report. 

3.30  MRED  -  A  System  for  Simulating  Single  Event  Effects  in 
Microelectronics 

The  key  component  of  our  single  event  simulation  system  is  a  Monte  Carlo  program  that  describes 
the  evolution  of  events  initiated  by  single  radiation  quanta.  The  program  is  called  MRED  (Monte 
Carlo  Radiative  Energy  Deposition)  and  its  structure  is  shown  in  the  following  figure. 

The  core  of  the  MRED  program  is  Geant4,  a  collection  of  C++  libraries  created  by  the  high 
energy  physics  community  for  describing  detectors  [3]  [4].  We  have  adapted  Geant4  for  use  with 
semiconductor  devices  through  a  unique  structure  based  on  the  computer  language  Python,  and 
the  addition  of  several  critical  Fortran  programs  for  specific  functionality.  Three  of  these  modules 
deal  with  nuclear-reaction  physics,  while  the  third  adds  state-state-of-the-art  electron,  positron,  and 
gamma  transport  down  to  <  100  eV  along  while  giving  plausible  results  for  low-energy  electron 
trajectories  in  deep  sub-micron  structures. 

The  C++  classes  that  were  available  in  Geant4  at  the  beginning  of  this  program  to  describe  ion-ion 
nuclear  processes  were  inadequate  to  support  research  on  single-event-effects  in  semiconductors. 
Leveraging  synergies  with  a  DTRA  Basic  Research  program  to  assess  nuclear  reaction  models, 
we  added  three  critical  Fortran  modules  to  MRED  and  have  used  them  extensively  in  single  event 
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Fig.  110.  A  collection  of  100  trajectories  of  250  eV  electrons  incident  on  a 
10  nm  cube.  Computed  with  MRED  using  the  interface  to  PENELOPE2008. 

This  image  demonstrates  the  capacity  of  MRED  to  establish  the  details  of 
energy  deposition  in  deep  sub-micron  structures,  one  of  the  objectives  of  this 
program.  Previously  unpublished  figure. 


analyses.  The  modules  are  JQMD,  a  component  of  the  PHITS  program  from  Japan  [5],  and  two 
programs  LAQGSM  and  CEM03  from  Los  Alamos  [6] . 

JQMD  was  initially  important  because  the  Geant4  libraries  did  not  include  physics  modules  for 
the  analysis  of  ion-ion  collisions  in  which  the  projectile  had  atomic  mass  greater  than  12  u.  Since 
a  substantial  number  of  cosmic  rays  are  in  this  category,  including  Fe,  significant  uncertainly 
was  introduced.  JQMD  is  a  quantum  molecular  dynamical  code,  which  treats  colliding  nuclei 
as  collections  of  nucleons  and  attempts  to  trace  their  motion  in  collision  and  to  identify  specific 
surviving  nuclear  species  following  the  collision.  Unfortunately,  the  upper  limit  for  JQMD  is 
approximately  3  GeV/u,  so  a  small  but  significant  fraction  of  high  energy  cosmic  rays  could  not  be 
described. 

To  broaden  the  energy  range  of  cosmic  rays  that  could  be  plausibly  described  we  turned  to  another 
Fortran  program,  LAQGSM,  created  at  Los  Alamos  National  Laboratory  from  an  earlier  code 
QGSM  that  originated  at  Dubna,  in  the  Soviet  Union.  After  a  protracted  negotiation,  we  became 
the  first  site  outside  Los  Alamos  to  obtain  a  license  for  use  of  this  software.  As  of  this  writing 
LAQGSM  and  a  related  program  for  protons  CEM03  are  the  default  modules  for  ion-ion  nuclear 
reactions  in  MRED. 

Although  the  general  quality  of  electron  transport  in  Geant4  is  high,  this  program  and  others  have 
specific  requirements  for  analysis  in  small  volumes.  Our  research  suggested  that  while  the  estimates 
of  energy  deposition  in  ultra-small  volumes  using  Geant4  alone  were  probably  relatively  accurate, 
the  details  of  trajectories  were  not.  Near  the  end  of  the  program  we  were  fortunate  to  obtain  access 
to  a  state-of-the-art  electron,  positron  and  photon  transport  program,  PENELOPE2008,  through  a 
collaborator,  Makoto  Asai,  at  the  Stanford  Linear  Accelerator  Laboratory  [7] .  The  following  figure 
shows  the  result  of  bombarding  a  10  nm  silicon  cube  with  100  250  eV  electrons.  The  detail  in  these 
trajectories  evident.  In  fact,  this  figure  includes  several  instances  where  secondary  electrons  with 
energies  <  100  eV  have  been  created. 
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A.  Emergence  of  a  Hierarchy  of  Rate  Prediction  Methods 


One  of  the  motivating  factors  for  this  program  was  the  suggestion  from  observations  on  hardened 
SRAMs  that  materials  such  as  W  surrounding  the  active  devices  could  play  a  role  in  establishing 
the  rate  of  observed  single  event  effects.  Using  capabilities  developed  under  this  program  and 
exploiting  synergies  with  NASA  and  DTRA  programs  we  were  able  to  add  addition  evidence  for 
the  basic  conclusion  that  non- silicon  materials  in  integrated  circuits  contribute  to  single  event  effect 
rates,  and  subsequently  to  demonstrate  experimentally  the  validity  of  the  hypothesis.  In  each  case, 
observation,  supporting  observation,  and  confirmation,  the  software  and  techniques  described  here 
proved  to  be  key  to  establishing  the  finding. 


The  following  figure  shows  schematically  the  relationship  between  various  components  of  the 
MRED  simulation  system.  The  Geant4  portion  of  the  program  corresponds  roughly  with  the  Virtual 
Irradiator  of  this  figure.  Energy  deposition  computed  using  the  virtual  irradiator  is  converted  to 
thermalized  charge  linearly  at  a  rate  of  3.6  eV/e-h-pair  and  one  of  several  methods  is  used  to 
determine  whether  a  single  event  has  occurred.  When  this  program  began,  we  believed  that  it  would 
be  possible  to  achieve  a  substantial  improvement  in  the  time  it  takes  for  a  Poisson  solver  to  make  a 
determination  that  a  given  quantum  has  or  has  not  caused  an  effect.  In  spite  of  significant  progress 
made  to  the  FLOODS  simulator  through  the  course  of  this  program,  it  has  still  not  been  possible 


Fig.  111.  Monte  Carlo  simulation  of  single  event  effects.  Energy  deposition  by  individual  radiation  quanta  is  generated  using 
the  MRED  program  and  analyzed  by  any  of  a  hierarchy  of  methods  to  determine  the  probability  of  a  single  event  effect.  The 
frequency  of  such  events  is  a  direct  measure  of  the  single  event  cross  section. 


to  introduce  detailed  physical  simulations  into  the  standard  single-event-analysis  workflow.  As  a 
result,  a  number  of  techniques  have  been  devised  through  this  and  related  programs  to  determine 
the  electrical  effects  of  a  specific  pattern  of  energy  deposition.  One  of  the  most  elaborate  and 
successful  has  been  the  introduction  of  SPICE  simulations  to  predict  the  circuit  effects  in  events  in 
which  several  transistors  simultaneously  receive  charge. 

The  use  of  SPICE  as  a  method  for  the  analysis  of  circuit  effects  is  an  example  of  one  constituent  of  a 
hierarchy  of  approximation  that  has  emerged  in  the  Monte  Carlo  the  analysis  of  single  event  effects. 
This  hierarchy  is  illustrated  in  the  following  figure.  Supplementing  the  conventional  rectangular 
parallelepiped  (RPP)  and  integral  RPP  (IRPP)  methods  are  a  number  of  additional  ones  that  involve 
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Monte  Carlo  solution.  The  simplest  of  these  is  a  calorimetry  in  a  single  sensitive  volume.  This  differs 
from  the  RPP  method  in  the  possible  shape  of  the  sensitive  volume,  the  materials  surrounding  the 
sensitive  volume  and  most  importantly  in  the  number  and  detail  of  the  physics  processes  that  are 
simulated. 

To  begin  to  approximate  collected  charge  rather  than  deposited  energy  a  multiple-sensitive-volume 
scheme  has  been  developed  that  involves  the  linear  combination  of  deposited  energy  values  in 
nested  volumes.  This  weighted  sum  can  be  calibrated  to  approximate  collected  charge  in  a  surprising 
number  of  situations. 

The  SPICE  method  referenced  above  uses  multiple  nested-sensitive-volume  groups  to  approximate 
individual  transistor  behavior,  but  determines  the  overall  effect  with  circuit  simulation.  This  scheme 
makes  it  possible  to  analyze  events  in  which  multiple  transistors  much  be  involved,  with  each  being 
perturbed  in  specific  ways. 


RPP 


IRPP 


MC  One  SV 


MC  multi-SVs  MC  SVs  +  SPICE  MC  SVs  +  TCAD 


Fig.  112.  Hierarchy  of  approximation  in  the  simulation  of  single  event  effects.  The  conventional  RPP  and  IRPP  methods  are 
supplemented  by  methods  involving  a  single  sensitive  volume,  a  group  of  sensitive  volumes  with  varying  efficiencies,  sensitive 
volume  groups  and  SPICE  analysis,  and  ultimately  the  use  of  a  Poisson  solver  to  track  charge  motion 


Finally,  at  the  most  detailed  a  full  Poisson  solution  can  be  applied  to  an  initial  condition  based  on 
the  actual  structure  of  the  radiation  event.  Presently,  this  method  is  simply  too  slow  to  be  practical, 
but  it  has  been  demonstrated  and  may  one  day  be  the  preferred  approach,  particularly  if  advanced 
methods  for  computing  the  evolution  of  the  event  can  be  developed  that  are  both  physically  accurate 
and  computationally  efficient. 
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B. 


Mathematical  Formulation  of  the  Rate  Prediction  Hierarchy 


The  capstone  achievement  of  the  single  event  simulation  portion  of  this  work  has  been  the 
establishment  of  a  mathematical  framework  for  describing  rate  prediction  that  is  sufficiently 
comprehensive  that  both  Monte  Carlo  and  analytical  approximations  can  be  handled  in  a  unified 
way.  Thus,  the  same  formal  analysis  can  be  applied  to  all  of  the  distinct  methods  shown  in  the 
above  figure. 

The  mathematical  basis  for  single-event  rate  prediction  has  been  comprehensively  reviewed  by 
Petersen.  However,  all  of  the  methods  that  he  describes  begin  from  starting  points  that  are  already 
relatively  specialized  to  the  assumptions  that  underlie  the  RPP  approach.  For  example,  all  begin 
with  a  direct  and  deterministic  relationship  between  energy  deposition  and  the  properties  of  the 
incident  particles.  However,  beginning  at  this  level  of  abstraction  excludes  approaches  in  which  the 
energy  deposition  is  handled  implicitly  or  indirectly.  Similarly,  it  is  not  easy  to  untangle  the  precise 
way  in  which  energy  enters  the  computation. 

We  have  analyzed  single  event  effects  (SEE)  rate  prediction  using  a  general  approach  that  allows 
the  rate  of  any  effect  initiated  by  a  single  radiation  quantum  to  be  computed,  given  the  addition 
of  sufficiently  detailed  assumptions.  These  assumptions,  formulated  in  concise  mathematical 
language,  reduce  the  problem  to  obtaining  the  solution  of  the  general  SEE-rate  integral. 


R, (£,  {}  -  J  )  ■  u )  H (-n (x )  df  EJKxj r) Fr (r,  E.  u,  j£ ,  /' ; g,  t} 

•:  At VE  -® 


There  is  no  universally  accepted  mathematical  notation  for  SEE  rate  computations,  and  even  the 
definitions  of  terms  are  sometimes  ambiguous.  In  the  following  discussion,  the  word  “event”  is 
used  to  describe  the  interaction  of  one  quantum  of  radiation  with  the  system  being  analyzed.  A 
process  of  interest  initiated  by  a  single  quantum  of  radiation  is  called  an  “effect.”  For  each  event, 
there  is  a  probability  of  an  effect.  The  basis  of  rate  prediction  is  the  calculus  of  these  probabilities. 

To  be  definite,  it  is  assumed  that  there  is  a  region  of  space  called  the  “world”  within  or  surrounding 
the  device,  circuit,  or  system  of  interest.  The  laboratory  and  sample  are  assumed  to  be  at  rest  in 
the  world  reference  frame.  The  size  of  the  world  is  unspecified  both  absolutely  and  relative  to  the 
device,  but  it  is  assumed  that  the  occurrence  of  an  effect  is  completely  controlled  by  what  happens 

within  the  world  and  that  the  world  surface  is  convex  with  surface  normal  unit  vector  m)  at 
position  -  ' . 

Radiation  is  assumed  to  have  two  independent  sources,  particles  crossing  the  surface  into  the  world 
volume  from  outside  it,  as  e.g.  cosmic  rays,  and  particles  originating  in  the  world  volume,  as  e.g. 
by  radioactive  decay.  The  external  radiation  source  is  completely  specified  at  every  point  on  the 

world  surface  by  the  flux  distributions  it  0(£n  E,  U  7  X,  .  The  internal  radiation  source  is  specified 

by  a  local  generation  rate  F(z,  E,  .V,  f"1)  .  Here  z  is  the  type  of  particle,  such  as  y  ray,  proton,  a 
particle,  56Fe,  etc.,  E,  is  the  particle  energy,  it  is  the  vector  specifying  particle  direction,  I  is  a  point 
or  origin  on  the  world  surface  or  within  the  world  volume,  and  r  ‘  is  the  time  at  which  the  particle 
is  at  location  .V  .  The  units  of  the  flux  are,  for  example,  particles  per  second,  per  square  centimeter, 
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per  MeV,  per  steradian.  The  units  of  the  volume  generation  rate  are  particles  per  second,  per  cubic 
centimeter,  per  MeV,  per  steradian.  A  typical  space  environment  contains  many  species,  z- 

The  physical  behavior  of  the  system  is  determined  by  the  probability  per  unit  time 

| J)  that  a  quantum  of  type  z,  which  has  initial  energy  E ,  direction  it  ,  and  position 

(Jt?  f^)  ,  will  result  in  an  effect  during  the  interval  tfl  at  time  t.  The  argument  |  represents  a  set  of 
configuration  parameters  on  which  this  probability  may  depend.  For  example,  in  the  RPP  model  the 
size  of  the  RPP  and  the  critical  energy  E  are  such  parameters.  Given  these  three  definitions  —  flux, 

generation  rate,  and  effect  probability  —  the  rate  /?.(!,  1}  of  single  event  effects  of  type  e  at  time 
t  for  configuration  parameters  0  is  given  by 

Here  H  is  the  Heaviside  unit  step  function.  Reading  from  the  inside  out,  the  integrals  in  the  first 
term  are  over  all  times  F'  prior  to  t ,  over  the  whole  of  the  world  surface  area,  over  all  particle 
directions,  and  over  all  particle  energies.  In  the  second  term,  the  surface  integral  is  replaced  by  a 
volume  integral  over  all  of  the  volume  enclosed  by  the  world  surface.  In  both  terms  the  sum  is  over 
all  species  in  the  respective  radiation  environment.  In  the  first  term,  the  step  function  constrains 

the  surface  and  angular  integrals  to  contributions  for  which  u  ^0  (that  is,  inward).  The 

minus  sign  appears  because  the  world  surface  is  described  by  an  outward  normal  vector,  while 
only  particles  on  the  world  surface  that  are  directed  inward  are  capable  of  producing  an  effect. 
The  second  term  contains  no  such  constraints  because  spontaneous  emission  of  radiation  in  any 
direction  within  the  volume  defined  by  the  world  surface  is  potentially  significant. 


Equation  (1)  is  directly  relevant  to  Monte  Carlo  simulation  of  single  event  effects,  since  the  process 
of  simulation,  particularly  involving  the  first  term,  in  essence  builds  up  the  rate  by  determining 
the  probability  P e  by  repetitive  sampling.  In  addition,  however,  with  appropriate  mathematical 
assumptions  equation  (1)  can  be  shown  to  reduce  to  the  expected  formulae  of  the  RPP  model, 
including  its  use  of  the  LET  distribution  in  place  of  energies. 
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M.L.  McLain  (graduate  student,  ASU)  received  the  2008  -  2009  Achievement  Award  for  College 
Scientists,  ARCS  Phoenix  Chapter. 
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